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Section  1 
INTRODUCTION 


Regions  of  shock  wave/turbulent  boundary  layer  interaction,  in  which  large 
pressure  gradients,  flow  separ  avion  and  turbulent  nonequilibrium  flows  are  generated, 
have  traditionally  been  the  testing  ground  of  prediction  techniques,  particularly  based 
on  the  solution  of  the  full  Navier-Stokes  equations  where  the  recirculating  flow  can 
be  handled  "more  exactly."  However,  despite  the  significant  advances  in  computational 
techniques  during  the  past  decade,  there  .emain  significant  gaps  in  our  understanding 
and  ability  to  predict  regions  of  shock  wave/tut  bulent  boundaty  interactions  in  hypersonic 
flows.  At  the  heart  of  the  problem  lies  the  difficulty  of  describing  the  generation  and 
development  of  turbulence  in  the  extremely  large  pressure  gradients  developed  across 
hypersonic  interaction  regions.  However,  many  of  the  problems  that  have  been  blamed 
on  poor  turbulence  modeling  may  in  fact  sterr.  from  the  grid  selection  or  nature  of 
the  numerical  scheme.  In  the  rush  to  demonstrate  that  the  Navier-Stokes  codes  can 
be  successfully  applied  to  describe  a  variety  of  interaction  problems,  very  little  emphasis 
has  been  placed  on  the  demonstration  that  the  numerical  schemes  are  indeed  an  accurate 
representation  of  the  equations  upon  which  they  are  based.  At  a  minimum,  the  sensitivity 
to  grid  size  should  be  examined,  and,  particularly  for  turbulent  interacting  flows,  an 
analysis  of  the  characteristic  scale  lengths,  like  that  employed  in  triple  deck  theory, 
should  be  performed  to  aid  in  grid  positioning. 

The  complexity  of  turbulent  interaction  regions  makes  it  essential  that  detailed 
information  from  experiments  be  used  to  construct  realistic  models  of  the  turbulent 
transport  mechanisms.  The  harsh  aerothermal  environment  in  hypersonic  turbulent 
interacting  flows  makes  delicate  probing  of  these  flows  difficult,  and  because  the  typical 
natural  frequencies  in  these  flows  are  of  the  order  of  500  kHz,  the  instrumentation 
and  recording  requirements  for  capturing  fluctuating  features  of  the  flow  are  far  from 
simple.  While  skin  friction  and  heat  transfer  measurements  are  useful  in  the  evaluation 
of  turbulence  modeling  techniques,  employing  comparisons  with  surface  pressure  data 
to  support  the  models  used  in  a  Navier-Stokes  code  is  a  weak  verification. 

The  complexity  of  the  flow  field  in  regions  of  shock  wave/turbulent  boundary 
layer  interaction  is  such  that  it  is  unrealistic  to  expect  to  describe  such  regions  in  any 
detail  within  the  framework  of  the  boundary  layer  equations.  Indeed,  there  are  some 
who  would  question  whether  the  time-  or  mass-avetaged  Navier-Stokes  equations  capture 
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the  basic  fluid  mechanics  associated  with  the  intrinsically  unsteady  nature  of  separated 
regions.  In  hypersonic  flows,  the  effects  of  compressibility  on  the  structure  and 
development  of  turbulence  must  also  be  considered. 

While  there  have  been  strenuous  efforts  to  obtain  predictions  of  two-dimensional 
and  three-dimensional  turbulent  interaction  regions,  it  is  currently  recognized  that 
"successes"  with  the  "Navier-Stokes"  code  in  describing  some  three-dimensional  turbulent 
interaction  regions  are  a  result  of  the  dominance  of  the  pressure  and  inertial  terms  in 
these  flows.  In  two  recent  studies, ’.2  j+  y,^  shown  that  turbulence  modeling  could  be 
changed  signficantly  without  significantly  changing  the  numerical  solution  for  a  three- 
dimensional  interaction  region.  For  two-dimensional  interactions,  it  appears  that 
turbdlence  modeling  is  more  critical.  To  obtain  good  agreement  for  thtse  latter  flows, 
some  very  gross  assumptions  must  be  made  in  the  turbulence  model.  Shang  and  Hankey^, 
for  example,  chose  to  apply  an  empirical  relationship  (selected  by  matching  the  length 
oi  the  separated  region)  to  rapidly  decrease  the  turbulent  scale  size  through  the 
interaction  region..  Horstmann^,  however,  found  the  best  agreement  with  Settles^ 
measurements  in  wedge-induced  separated  regions  using  a  two-equation  model  for 
turbulence  scale  size  and  vorticity.  Working  with  this  same  turbulence  model,  however, 
Horstmann  was  unable  to  predict  the  occurrence  of  separation  on  two  incident 
shock/turbulent  boundary  layer  configurations  studied  by  Holden**  at  Mach  11.2.  In  the 
latter  studies,  for  flows  which  were  clearly  separated,  tire  numerical  solution  failed  to 
predict  the  characteristic  plateaus  in  either  the  heat  transfer  or  pressure  distributions. 
The  modeling  of  turbulence  in  separated  interaction  regions  at  hypersor.ic  Mach  numbers 
should  account  for  the  effects  of  compressibility  and  the  generation  of  turbulence  by 
the  unsteady  movement  of  the  incident  and  induced  shocks  as  they  traverse  and  interact 
with  a  major  region  of  the  turbulent  boundary  layer  and  the  unsteady  movement  of 
the  separated  region.  Clearly,  further  detailed  experimental  work  on  insightful 
theoretical  modeling  is  required  to  develop  numerical  prediction  techniques  which  are 
capable  of  describing  separated  turbulent  interaction  regions  in  hypersonic  flow. 

In  this  report,  Section  2  describes  the  use  of  "parabolized"  (or  "thin-layer") 
Navier-Stokes  equations  to  solve  turbulent  flow  equations  and  existing  differing  schemes 
to  solve  Navier-Stokes  equations.  Section  3  describes  how  models  and  instruments  were 
developed  and  used  to  obtain  detailed  measurements  of  the  profile  characteristics  of 
a  turbulent  boundary  layer  ahead  of  and  through  regions  of  flow  separation  induced  by 
shock  wave/boundary  layer  interaction  over  a  cone/flare  configuration.  Section  4 
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presents  the  results  of  the  application  of  laser  holography  to  measure  viscous  flow 
interactions  in  hypersonic  shock  tunnel  flows.  Section  5  summarizes  the  research 
conducted  under  the  current  contract  to  investigate  fundamental  problems  associated 
with  flight  at  hypersonic  speeds  with  emphasis  placed  on  aerothermal  effects  of 
viscous/inviscid  interactions  and  boundary  layer  transition.  Section  6  provides  a  list  of 
references.  Two  appendices  provide  copies  of  papers  on  aerothermal  problems  associated 
with  hypersonic  flight  that  were  presented  by  Dr.  M.S.  Holden  at  the  AIAA  24th 
Aerospace  Science  Meetings. 


3 


Section  2 

SUMMARY  OF  EXISTING  SOLUTION  TECHNIQUES 
2.1  INTRODUCTION 

While  the  Navier-Stokes  equations  were  formulated  nearly  140  years  ago,  only 
recently  has  computational  speed  increased  enough  to  handle  the  entire  set  of 
simultaneous  differential  equations.  Even  with  the  advent  of  the  new  generation  "super" 
computers  (e.g.,  CRAY1,  CRAYXMP,  CDC7600,  new  generation  array  processors,  etc.), 
economics  remains  a  major  consideration.  Chapman?,  however,  notes  that  the  cost  of 
these  types  of  calculations  has  been  decreasing  by  a  factor  of  ten  every  eight  years. 
It  may  be  just  a  matter  of  time  before  these  types  of  large  numerical  simulations  ate 
of  some  consequential  value  to  practical  engineering  problems. 

In  most  turbulent  flow  applications,  the  Navier-Stokes  equations  are  simplified 
by  Reynolds  time  averaging.  We  are  then  faced  with  the  difficult  task  of  developing 
a  model  for  turbulence.  There  have  been  many  approaches  to  the  turbulent  closure 
problem,  from  simple  algebraic  closure  to  various  two-equation  models  (e.g.,  k-e,  k-w, 
see  Anderson  and  Tannehill^,  Tennekes  and  Lumley^,  Jones  and  Launder1®),  and  the 
Reynolds  stress  equations.  The  lack  of  progress  is  threefold.  One,  we  lack  the  physical 
understanding  to  model  these  complex  turbulent  flows.  Two,  we  lack  the  computational 
resources  to  deal  with  so-called  higher  order  turbulence  models,  which  add  nine  more 
unknowns  (in  an  incompressible  flow)  and  two  or  three  more  equations.  Three,  we  still 
lack  the  ability  to  solve  the  laminar  portion  of  the  Navier-Stokes  equations  in  an 
efficient  and  accurate  manner.  Most  existing  differencing  schemes  (which  will  be 
discussed  here)  are  second-order  accurate  with  various  restrictions  applied  by  the 
numerical  stability. 

There  are  several  levels  of  complexity  of  Navier-Stokes  equations  which  fall 
between  a  full  Navier-Stokes  equation  (Figure  1)  and  a  boundary  layer  equation.  Common 
names  for  these  types  of  equations  are  "parabolized"  Navier-Stokes  equations  and  "thin- 
layer"  Navier-Stokes  equations.  These  sets  of  equations  are  characterized  by  their 
applicability  in  both  viscous  and  inviscid  flow  regions;  they  all  contain  the  normal 
pressure  te*ms  which  are  usually  dropped  in  any  boundary  layer  formulation. 
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To  arrive  at  the  thin  layer  approximation  shown  in  Figure  1,  a  scaling  argument 
is  used  (see  Anderson  and  TannehillS),  and  terms  on  the  order  of  l/Re^  and  smaller 
are  neglected.  All  viscous  terms  containing  derivatives  parallel  t«.  the  wall  ate  dropped 
because  they  are  substantially  smaller  than  those  viscous  terms  normal  to  the  wall.  If 
this  same  argument  were  to  be  used  for  the  normal  (y)  direction,  the  equation  could 
be  reduced  to  dp/dy=0.  For  this  approximation  to  handle  a  possible  flow  reversal,  it 
is  necessary  to  retain  all  terms  in  the  momentum  equations  except  for  the  viscous 
terms  oarallel  to  the  flow  (the  upstream  diffusion  terms).  The  final  set  of  equations  is 
as  shown  in  Figure  1,  in  the  Cartesian  coordinate  system.  This  mixed  set  of  hyperbolic- 
parabolic  partial  differential  equations  is  solved  instead  of  the  elliptic-hyperbolic  partial 
differential  equations,  using  a  "time-dependent"  approach.  The  overall  complexity  of 
the  equations  is  reduced,  but  the  difficulty  of  solving  a  hyperbolic  system  still  remains. 

The  parabolized  Navier-Stokes  equations  are  derived  using  an  argument  similar 
to  that  employed  to  derive  the  thin-layer  equations.  An  asymptotic  expansion  (see 
Karvorkin  and  Cole**)  is  performed  (Rudman  and  Rubin^)  with  the  Navier-Stokes 
equations,  and  all  of  the  terms  exhibiting  characteristic  order  of  magnitude  of  (1/6*2) 
ur  higher  are  retained.  The  norma,'  pressure  gradient  terms  are  deleted  from  the  x- 
momentum  equation,  which  restricts  the  solutions  of  the  flows  without  upstream 
influences.  The  most  commonly  used  form  of  the  parabolized  Navier-Stokes  equations 
(formulated  by  Cheng^)  contains  the  streamwise  pressure  gradient  term.  Thus,  the 
mosx  common  form  of  this  type  of  equation  is  obtained  by  assuming  that  only  the 
streamwise  viscous  derivative  terms  and  the  heat  flux  terms  are  negligible  compared 
to  the  normal  transverse  viscous  derivative  terms.  The  parabolized  Navier-Stokes 
equations  are  derived  by  dropping  all  viscous  and  heat  flux  terms  containing  partial 
derivatives  with  respect  to  the  streamwise  direction,  from  a  steady  Navier-Stokes 
equation.  A  set  of  parabolized  Navier-Stokes  equations  is  shown  in  Figure  1. 

2.2  COMPARISONS  WITH  MEASUREMENTS  IN  HYPERSONIC  FLOW 

The  majority  of  the  efforts  made  in  the  numerical  simulations  at  CUBRC  were 
geared  toward  understanding  and  utilizing  existing  differencing  schemes  for  solving  the 
Navier-Stokes  equations.  The  initial  efforts  included  obtaining  a  code  that  utilized  an 
explicit  type  of  differencing  scheme  and  running  it  to  simulate  a  set  of  laminar 
shock/boundary  »ayer  interaction  problems.  The  MacCormack-Shang-Hankey  code  was 
obtained  from  J,  Shang  of  Wright-Patterson  AFB.  It  is  a  straightforward  application 
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-  ASSUMPTION  MADE  WAS  THAT 
ALL  VISCOUS  TERMS  CONTAINING 
DERIVATIVES  PARALLEL  TO  THE 
WALL  ARE  DROPPED. 


-  PARABOLIC  TYPE  OF  EQUATIONS  WILL 
NORMALLY  DROP  ($5)  TO  ACHIEVE 
PARABOLIC  NATURE,  BUT  IT  WAS  LATER 
ADDED  TO  ACHIEVE  MORE  GENERAL 
FORMULATION.** 

-  NORMALLY  WRITTEN  IN  A  FORM  WHERE  THE 
VI SCOUS  AND  INVISCID  TERMS  ARE 
SEPARATED.** 


•ALL  EQUATIONS  ARE  WRITTEN  IN  TWO-DIMENSIONAL  FORM. 
•SEE  TANNEHILL  (19841  FOR  FURTHER  DETAIL. 


Figure  1  CLASSIFICATION  OF  NAVIER-STOKES  EQUATIONS 
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SCHEME  PROPOSED  TO  DATE 


•ONLY  KNOWN  APPLICATION  TO  DATE  WAS  FOR  THE  PNS. 

••APPLIED  AS  EULF.R  SOLVER  BY  JAMESON  (19861;  NOT  DISCUSSED  IN  PRESENT  REPORT 

•••ALL  METHODS  ARE  APPLICABLE  TO  BOTH  FULL  COMPRESSIBLE  N  S  EQUATIONS  AND  ANY 
OTHER  SIMPLER  VERSIONS  OF  N*. 


Figure  1  CLASSIFICATION  OF  NAVIER-STOKES  EQUATIONS  (Cont.) 
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of  the  MacCormack  explicit  differencing  scheme,  which  was  vectorized  to  use  the 
CRAZ's  high-speed  vector  processing  (Shang^).  There  were  some  modifications  made 
on  the  pressure  damping  terms  to  handle  the  sharp  leading  edge  problem. 


The  governing  equations  are  in  vector  form: 
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Then  the  governing  equations  are  transformed  into  an  orthogonal  coordinate 
system  based  on  generalized  transformation,  as  proposed  by  Viviand*^  which  made  it 
easier  to  handle  arbitrary  geometry. 
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The  numerical  scheme  applied  is,  as  noted  previously,  MacCormack's^  explicit 


scheme. 
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The  numerical  stability  condition  proposed  by  MacCormack^  is  also  utilized  by 


this  code. 
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The  fourth-order  pressure  damping  was  used  to  add  a  dissipative  control  to  the 
numerical  instabilities  as  proposed  by  ShangA 
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Thus  far.  we  have  made  a  considerable  effort  in  trying  to  obtain  a  set  of  solutions 
for  Holden's^  laminar  shock/boundary  layer  interaction  problem  for  which  a  similar 
type  of  numerical  study  already  exists  (Hung  and  MacCormack^O).  After  a  considerable 
amount  of  "extra"  effort,  we  were  able  to  obtain  a  solution  for  the  flat  plate  section 
of  the  problem.  These  "extra"  efforts  include  applications  of  two  new  mesh  systems: 
or.e  mesh  system  near  the  wall  to  resolve  the  length  scale  problem  within  the  viscous 
layer,  and  the  second  mesh  system  near  the  leading  edge  of  the  flat  plate  to  resolve 
length  scales  near  a  sharp  leading  edge.  The  leading  edge  mesh  system  was  on  the 
order  of  10"^  feet  in  size,  and  the  wall  mesh  system  was  10~^  feet,  which  proved  to 
be  small  enough  to  resolve  all  required  length  scales  within  the  problem.  The  sharp 
leading  edge  problem  is  especially  difficult  for  explicit  type  schemes,  but  unless  it  is 
resolved,  the  explicit  calculations  will  not  be  able  to  proceed.  The  explicit  schemes 
require  that  the  time  marching  is  done  at  the  smallest  length  scale  of  the  problem, 
as  stated  previously,  and  if  the  smallest  length  scale  is  not  resolved,  time  marching 
will  cause  the  numerics  to  overshoot  the  physics,  resulting  in  a  numerical  instability. 
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Once  these  length  scales  were  resolved,  we  were  able  to  obtain,  with  considerable 
computational  time,  a  set  of  solutions  for  the  flat  plate  portion  of  the  shock  boundary 
interaction  problem.  The  calculated  skin  friction  and  heat  transfer  data  are  shown  in 
Figures  2  and  3.  These  figures  show  good  agreement  with  Holden's  early  experimental 

data.^OA 
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Section  3 

EXPERIMENTAL  PROGRAM 


3.1  PROGRAM  OBJECTIVE  AND  DESIGN  OF  THE  EXPERIMENTAL  STUDY 

The  major  objective  of  the  current  study  was  to  develop  and  use  models  and 
instrumentation  to  obtain  detailed  measurements  of  the  profile  characteristics  of  a 
turbulent  boundary  layer  ahead  of  and  through  regions  of  flow  separation  induced  by 
shock  wave/boundary  layer  interaction  over  a  cone/flare  configuration.  Such 
measurements  are  of  key  importance  in  the  evaluation  of  the  theoretical  modeling  of 
the  turbulent  separation  process  in  hypersonic  flows.  As  discussed  in  the  introduction, 
current  turbulence  models  apparently  are  incapable  of  describing  the  d  elopment  of 
turbuience  in  regions  of  strong  pressure  gradients  and  boundary  layer  separation  in 
hypersonic  flow,  possibly  because  of  compressibility,  shock/turbulence  interaction,  or 
unsteady  effects  under  hy*>et  sonic  highly  cooled  wall  conditions.  In  hypersonic  high- 
Reynolds-number  flows  over  highly  cooled  walls,  the  "wall  layer",  in  which  our  earlier 
studies  have  suggested  separation  first  takes  place,  and  which  contains  the  principal 
information  on  the  character  of  the  boundary  layer,  is  very  thin.  Consequently, 
boundary  layer  thicknesses  of  ovet  one  inch  are  required  to  enable  this  layer  to  be 
probed  with  the  required  resolution.  While  the  turbulent  boundary  layers  on  the  walls 
of  hypersonic  nozzles  have  been  used  as  the  source  of  thick  turbulent  boundary  layers 
in  experimental  studies,  it  has  been  shown  that  significant  turbulent  nonequilibrium 
effects  can  exist  in  these  nozzle  flows^.  The  distortion  of  the  structure  and  turbulent 
characteristics  of  the  boundary  layer  generated  through  the  strong  expansion  in  the 
nozzle  can  persist  well  downstream  of  the  nozzle  exit  plane  and  can  significantly 
influence  the  characteristics  o  1  a  separating  turbulent  boundary  layer.  For  this  reason, 
we  elected  to  perform  studies  to  examine  the  characteristics  of  the  turbulent  boundary 
ahead  of  and  in  regions  of  shock  wave/boundary  layer  interaction  on  a  large  slender 
cone/flare  configuration  in  the  large  contoured  "D"  nozzle  in  the  96-inch  Shock  Tunnel. 
The  technique  employed  in  the  design  of  the  "D"  nozzle,  and  indeed,  most  contoured 
nozzles,  is  such  that  the  test  core  is  a  cone-shaped  region  of  uniform  flow  which 
originates  well  upstream  of  the  exit  plane.  Thus,  by  designing  a  conical  model  so  that 
it  can  be  fit  within  this  uniform  conical  region,  it  is  possible  to  develop  a  constant 
pressure  boundary  layer  over  a  large  conical  model  that  extends  well  into  the  contoured 
nozzle.  The  ultimate  objective  of  this  study  was  to  obtain  both  mean  and  fluctuation 
measurements  on  the  surface  and  across  the  turbulent  layer.  However,  during  this 
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phase  of  the  study,  we  concentrated  on  obtaining  measurements  of  the  mean  properties 
across  the  viscous  layer,  more  specifically  to  obtain  measurements  of  the  pitot  pressure, 
total  temperature,  total  heat  transfer  rate  and,  using  holographic  interferometry,  the 
mean  density  distribution. 

3.2  MODELS  AND  INSTRUMENTATION 

The  experimental  studies  were  conducted  at  Mach  numbers  of  11  and  13  in  the 
large  contoured  "D"  nozzle  in  the  Calspan  96-inch  Shock  Tunnel.  As  discussed  in 
subsection  3.1,  the  large  conical  region  of  uniform  flow  that  extends  well  up  into  the 
contoured  "D"  nozzle  allows  us  to  generate  a  constant  pressure  boundary  layer  on  a 
conical  model  which  extends  into  the  nozzle. 

For  these  studies,  we  selected  the  large  6-degree  cone  with  flares  of  30  and  36 
degrees  attached  at  its  base.  The  cone/flare  configuration  is  shown  in  Figure  4.  The 
cone  angle  and  length  were  selected  on  the  basis  of  calculations  to  achieve  the  maximum 
length  over  which  uniform  constant  pressure  flow  could  be  established  within  the  further 
constraints  of  tunnel  blockage  and  sting  loading.  A  diagram  of  the  cone/flare  model 
and  its  positioning  within  the  "D"  nozzle  is  shown  in  Figure  5.  Previously, 21  we 
demonstrated  that  this  large  model  could  be  used  to  produce  the  required  flow,  we 
obtained  pressure  and  heat  transfer  measurements  for  this  model  equipped  with  both 
sharp  and  blunt  nosetips.  The  good,  agreement  between  the  measured  pressure  and  heat 
transfer  distribution  and  theory  for  these  configurations,  shown  in  Figures  6  through  9, 
demonstrates  lhat  the  design  and  positioning  of  the  model  produced  the  required  testing 
environment.  Schlieren  photographs  of  the  flow  field  for  these  cases  are  shown  in 
Figures  10  and  1J. 

3,2. 1  Heat  Transfer  Instrumentation 

Platinum  thin-filrr.  instrumentation  was  used  to  obtain  heat  transfer  measurements 
on  the  surface  of  the  flat  plate/cone  rncdel  and  as  the  sensing  element  of  the  0.0,5-Inch- 
diameter  stagnation  heating  probes,  because  this  gage  has  a  megahertz  frequency 
response,  it  can  be  used  to  examine  the  unsteady  characteristics  of  the  turbulent 
boundary  layer  and  separated  region.  The  large  gradients  that  are  generated  along  the 
walls  and  in  the  flow  in  the  separation  and  reattachment  region  of  shock  wave/boundary 
layer  interactions  make  it  essential  that  distortion  of  the  heat  transfer  distribution 
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Figure  4  SHARP  6°  CONE/300  FLARE  MODEL  INSTALLED  IN  CALSPAN'S  96 
SHOCK  TUNNEL 
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resulting  from  lateral  heat  conduction  be  minimized  by  employing  models  constructed 
with  low  conductivity  materials.  The  pyrex-backed,  thin-film  gages  with  their  high 
resolution,  sensitivity  and  frequency  response  are  almost  ideal  for  this  type  of  study. 
The  platinum  thin-film  probes  (Figure  12)  were  used,  in  conjunction  with  the  total 
temperature  instrumentation,  to  examine  the  structure  of  the  turbulent  boundary  layer 
and  shear  layer. 

3.2.2  Pitot  and  Static  Pressure  Instrumentation 

We  used  Calspan  piezoelectric  pressure  transducers  mounted  in  pitot  pressure 
rakes,  beneath  orifices  in  the  model  surface,  to  obtain  the  mean  pitot  pressure  through 
the  boundary  layer  and  along  the  surface.  The  pitot  pressure  gages  (Figure  12)  had 
G.030-inch  orifices  and  were  staggered  to  achieve  a  transverse  spacing  of  0.010  inch 
at  the  base  of  the  boundary  layer. 

3.2.3  Total  Temperature  Instrumentation 

A  significant  effort  was  devoted  to  the  design  and  development  of  a  total 
temperature  gage  which  responded  within  3  milliseconds,  withstood  the  large  static  and 
dynamic  pressures  generated  in  regions  of  shock/boundary  layer  interaction  in  the  shock 
tunnel  flows,  and  was  small  enough  to  resolve  the  total  temperature  in  the  wall  layer. 
The  result  of  this  development  was  a  gage  0.030  inch  in  diameter  which  used  a 
0.0005-inch,  butt-welded  iron/constantin  thermocouple  in  the  arrangement  shown 
schematically  in  Figure  13.  The  typical  response  of  one  of  these  gages  (Figure  14) 
clearly  shows  that  we  have  adequate  time  to  obtain  accurate  measurements.  A  small 
amount  (~2%)  of  radiation  is  applied  and  this  factor  is  checked  for  measurements  in 
the  freestream  where  the  total  temperature  is  known  accurately.  The  gages  are 
deployed  in  a  staggered  array  (Figure  12)  similar  to  that  employed  for  the  pitot  pressure 
gages. 


3.2,4  Test  Conditions  and  Model  Configurations 

The  experimental  studies  were  conducted  at  Mach  11,  13  and  16  for  Reynolds 
numbers  from  30  x  10^  to  80  x  106-  Under  these  conditions,  the  boundary  layer  is 
fully  turbulent  well  upstream  of  the  cone/flare  junction  and,  as  di  cussed  in  the  following 
sections,  the  measurements  of  heat  transfer  were  in  good  agreement  with  prediction 
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techniques  based  on  a  large  amount  of  measurements  on  highly  cooled  walls  in  high- 
Reynolds-number  hypersonic  flows.  The  test  conditions  at  which  the  studies  were 
conducted  are  listed  in  Table  1. 

3.3  RESULTS  AND  DISCUSSION 

The  experimental  program  had  two  objectives.  First,  we  sought  to  generate  a 
model  and  environment  in  which  we  could  establish  a  thick,  well-developed  turbulent 
boundary  layer  that  had  developed  under  constant  pressure  conditions.  We  then  sought 
to  design,  develop  and  use  instrumentation  to  obtain  profile  measurements— first,  in  the 
constant  pressure  boundary  layer,  and  then,  in  regions  of  strong  adverse  pressure  gradient 
in  regions  of  shock  wave/boundary  layer  interaction  generated  at  a  cone/flare  junction. 
We  planned  to  use  two  cone/flare  junctions:  one  (a  30-degree  flare)  which  promoted  a 
flow  close  to  incipient  separation,  and  a  second  (a  36-degree  flare)  which  promoted  a 
well  separated  flow. 

The  measurements  of  the  heat  transfer  and  pressure  distributions  over  the  two 
cone/flare  configurations  are  shown  in  Figures  6  through  9.  The  corresponding  schlieren 
photographs  are  presented  in  Figures  10  and  li.  The  measurements  of  the  pressure 
along  the  entire  length  of  the  cone  were  in  good  agreement  with  predictions  based  on 
Sim's  solutions  for  a  sharp  cone,  and  the  pressures  at  the  back  of  the  flare  were  in 
good  agreement  with  calculations  based  on  an  inviseid  shock  compression  from  the  cone 
to  the  flare  (Figures  7  and  8), 

Flow  fielo  surveys  were  made  to  determine  the  distribution  of  pitot  pressure, 
total  temperature  and  total  heat  at  a  number  of  stations  through  the  interaction  region 
at  each  of  the  flow  conditions  described  above.  Figures  15  and  16  show  the  pitot  and 
total  temperature  measurements  for  the  Mach  11  condition  with  a  30-degree  flare. 
The  profiles  that  were  obtained  at  2,  1.2,  0.8,  and  0  inches  ahead  of  the  cone/flare 
junction  indicate  that  there  is  very  little  upstream  influence  at  this  condition.  (A 
similar  set  of  measurements  for  the  36-degree  flare  is  shown  in  Figures  17  and  18.) 
The  Mach  number  and  velocity  profiles  shown  in  Figures  19  and  20  were  determined 
from  the  inverse  relationships: 
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Table  1 

TEST  CONDITIONS,  LARGE  6°  CONE 


RUNS 

r 

n2 

3,4,8 

6,7 

9 

6-15 

Mi 

3.345E+00 

3.633E+00 

4.200E-fO0 

2.63SE+00 

PO  PSIA 

7.216E-KJ3 

1.760E+04 

1.70SE+04 

R.430E+03 

HO  FT2/SEC2 

1.825E+07 

2.147E+07 

2.795E-HJ7 

1.287E+07 

TO  °R 
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Figure  15  DISTRIBUTION  OF  PITOT  PRESSURE  ACROSS  BOUNDARY  LAYER  UPSTREAM 
OF  CONE/FLARE  JUNCTION,  30°  FLARE 
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Figure  16  DISTRIBUTION  OF  TOTAL  TEMPERATURE  ACROSS  BOUNDARY  LAYER 
UPSTREAM  OF  CONE/FLARE  JUNCTION,  30°  FLARE 
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Fiqure  17  DISTRIBUTION  OF  PITOT  PRESSURE  ACROSS  BOUNDARY  LAYER 
UPSTREAM  OF  CONE/FLARE  INTERACTION,  36°  FLARE 
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Figure  18  DISTRIBUTION  OF  TOTAL  TEMPERATURE  ACROSS  BOUNDARY  LAYER 
UPSTREAM  OF  CONE/FLARE  JUNCTION,  36°  FLARE 
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Figure  19  DISTRIBUTION  OF  MACH  NUMBER  ACROSS  BOUNDARY  LAYER 
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These  measurements  are  compared  with  the  Crocco  relationship  between  enthalpy 
and  velocity  in  Figure  21.  It  is  clear  that  our  measurements  follow  a  parabolic 
relationship: 
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rather  than  Crocco's  linear  relationship: 


where: 
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In  the  past,  it  has  been  assumed  that  the  "fuller  than  Crocco"  velocity  profile  obtained 
in  studies  over  tunnel  walls  was  associated  with  turbulent  nonequilibrium  effects 
associated  with  the  strong  favorable  pressure  gradient  upstream  on  the  nozzle  wall;21 
however,  no  such  explanation  can  be  advanced  to  explain  our  results. 

Of  the  transformation  techniques  that  have  been  postulated  to  cast  the 
compressible  flow  measurements  into  an  equivalent  incompressible  form,  the  relationships 
derived  by  Van  Driest  have  received  the  greatest  support.  Van  Driest  starts  with  the 
assumption  that  the  Crocco  lelationship  is  valid  and  uses  the  mixing  length  theory  to 
calculate  the  Reynolds  stresses  in  the  flow.  By  inspection,  the  transformation,  is 
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Figure  21  TOTAL  TEMPERATURE  AND  VELOCITY  MEASUREMENTS  PRESENTED 
IN  CROCCO  FRAMEWORK 


where: 
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The  appropriate  form  of  the  incompressible  law  of  the  wall  relationship  is 


However,  because  the  wake  region  of  the  flow  is  so  extensive,  Maise  and  McDonald-22 
have  suggested  that  Coles2^  wake  function  be  included: 


They,  in  fact,  suggest  a  defect  form  of  this  relationship: 


•*  * 


Comparisons  between  our  measurements  in  each  of  the  incompressible  formats  are 
shown  in  Figures  22  and  23.  It  is  clear  that  the  transform  is  not  as  successful  at 
these  high  Mach  numbers  as  it  has  been  below  Mach  5.  Comparisons  between  the 
measurements  and  these  prediction  techniques  in  the  compressible  plane  are  shown  in 
Figure  2k.  It  can  be  seen  that  this  is  a  relatively  insensitive  way  of  examining  the 
measurements. 


Wi.  are  concerned  with  the  differences  between  the  current  set  of  measurements 
and  the  prediction  techniques.  We  will  attempt  to  check  our  measurements  against 
more  sophisticated  prediction  techniques  and  density  measurements  obtained  with 
holographic  interferometry.  We  also  plan  further  refinement  of  the  instrumentation  to 
more  effectively  probe  the  wall  layer. 
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Section  4 


PRELIMINARY  APPLICATIONS  OF  HOLOGRAPHIC  INTERFEROMETRY 
TO  STUDY  HYPERSONIC  REGIONS  OF  SHOCK  WAVE/ 

BOUNDARY  LAYER  INTERACTION 

4.1  INTRODUCTION 

This  section  documents  the  findings  of  an  initial  application  of  pulse  laser 
holography  to  study  hypersonic  turbulent  boundary  layer  interactions,  especially  in  regions 
of  strong  pressure  gradients.  One  important  objective  is  to  provide  detailed  information 
on  the  mechanism  of  flow  separation  in  regions  of  shock  wave/  boundary  layer  interaction 
in  hypersonic  flows.  In  previous  studies,  miniature  probes  were  developed  and  used  to 
examine  the  structure  of  turbulent  boundary  layers  in  high  Mach  number,  high  Reynolds 
number  hypersonic  flows.  The  extremely  large  heat  transfer  rates  and  stresses  developed 
on  the  tiny  probes  in  these  flows  make  this  task  difficult,  and  wall  interference  problems 
near  the  model  surface  present  important  concern  for  the  accuracy  of  the  measurements. 
Hence,  exploring  these  flows  with  non-intrusive,  optical  techniques  is  very  desirable. 

Holographic  interferometry  is  a  potentially  worthwhile  measuring  technique, 
because  it  is  nonintrusive  and  offers  the  potential  to  capture  information  of  the  entire 
flow  field  in  a  single  instantaneous  exposure.  To  the  uninitiated,  the  term  "holographic" 
in  holographic  interferometry  might  be  thought  to  imply  that  this  technique  is  capable 
of  measuring  densities  in  streamwise  slices  of  the  flow.  However,  this  Is  not  the  case, 
because  the  fringe  shifts  result  from  an  integration  of  the  density  changes  along  the 
path  traversed  by  the  laser  light  waves.  Thus,  as  in  conventional  Mach-Zehnder 
interferometry,  distortions  resulting  from  two-dimensional  effects,  surface  refraction, 
diffraction  and  more  importantly,  the  intrinsic  unsteady  three-dimensional  structure  of 
a  turbulent  flow  can  result  in  uncertainty  in  the  interpretation  of  the  optical 
measurements.  Edge  effects  are  of  special  concern  when  this  technique  is  used  to 
examine  separated  flows  over  two-dimensional  models,  because  three-dimensional  effects 
are  induced  by  the  separated  region  at  the  edges  of  the  model. 


4.2  LASER  HOLOGRAPHIC  SYSTEM 


The  holographic  recording  system  2^  (Figure  25)  was  loaned  to  CUBRC  (Calspan- 
UB  Research  Center,  Buffalo,  NY)  by  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL), 
Wright-Patterson  Air  Force  Base,  Ohio. 


The  holograms  obtained  during  this  study  are  made  with  a  pulsed  ruby  laser  that 
is  Q-switched  passively  with  a  dye  cell  to  produce  single  25-nanosecond  (nominal)  light 
pulses.  Both  single  plate  and  dual  plate  techniques  are  required  to  record  holograms, 
which  are  subsequently  used  in  the  playback  step  to  obtain  shadowgrams,  schlieren 
photographs  and  interferograms  of  the  shock  tunnel  tests.  The  optical  reference 
conditions  are  taken  for  both  atmospheric  and  tunnel  "pump  down"  conditions;  the  latter 
condition  is  needed  to  assess  the  optical  effect  of  loading  the  test  chamber  windows. 
During  the  actual  recording  of  the  holograms,  the  shock  tunnel  test  section  is  nearly 
at  vacuum  condition,  which  means  that  the  windows  are  stressed  by  atmospheric  pressure. 


4.2.2  Shock  Tunnel  Installation 

The  holographic  system  is  positioned  and  aligned  with  the  Toepler  schlieren 
system  of  the  shock  tunnel  facility.  A  20-power  microscope  objective  is  used  to  expand 
the  ...ser  oeam  to  fill  a  15-inch  diameter  parabolic  mirror.  A  pin  hole  is  not  required 
to  spatially  filter  the  laser  beam.  A  collimated  beam  reflected  from  the  parabolic 
mirror  is  directed  throjgh  the  windows  of  the  shock  tunnel,  captured  by  a  second 
parabolic  reflector  on  the  opposite  side  of  the  tunnel,  and  directed  toward  the  hologram 
recording  plane.  The  holographic  reference  beam  is  split  from  the  output  beam  of  the 
ruby  laser  by  a  3Q  to  70  percent  beam  splitter.  Thirty  percent  of  the  incident  radiation 
is  reflected  from  the  front  surface  of  the  beam  splitter;  this  light  is  used  to  establish 
the  scene  beam,  and  the  remaining  70  percent  is  directed  over  the  shock  tunnel  by 
specially  built  beam  elevation  towers.  The  holographic  reference  beam  is  collimated  by 
a  beam  expander-collimator  positioned  approximately  7  inches  ahead  of  the  hologram 
recording  plane.  Natural  divergence  of  the  ruby  laser  beam  causes  an  approximate  50 
percent  reduction  in  the  intensity  of  the  light  entering  the  collimator,  so  the  net 
intensities  of  the  scene  and  reference  beams  in  the  hologram  plane  are  nearly  equal. 
The  AFFDL  dual  plate  hologram  holder^  is  used  to  hold  and  align  the  holograms,  and  it 
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is  positioned  so  that  the  converging  light  waves  of  the  scene  beam  illuminate  a  1-inch 
circular  area  in  the  hologram  plane;  the  collimated  reference  waves  illuminate  a  2- 
inch  circular  area,  and  both  beams  are  aligned  relative  to  the  hologram  rac^rding  plane 
so  that  the  light  waves  of  the  scene  beam  are  centered  upon  the  light  waves  of  the 
reference  beam.  Further,  the  holder  is  positioned  so  that  the  incident  angles  of  both 
beams  are  20  degrees  relative  to  normal  to  the  holograms.  The  path  length  of  each 
beam  is  approximately  53  feet,  and  both  path  lengths  are  adjusted  to  be  equal  to  within 
2  inches. 
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The  ruby  laser  in  the  recording  system  is  a  Korad,  single-oscillator  system  that 
is  passively  Q-switched  with  a  Kryptocine  dye  cell.  The  bank  voltage  applied  to  the 
flash  lamp  and  the  dye  concentration  in  the  dye  cell  are  adjusted  so  that  the  laser  emits 
a  single  pulse  energy  of  25  millijoules  in  25  to  50  nanoseconds,  nominally.  The  output 
beam  is  0.06  inch  in  diameter  at  the  front  reflector  of  the  ruby  laser  cavity.  The  laser 
is  remotely  charged,  and  once  armed,  the  laser  firing  circuit  is  triggered  from  a  voltage 
spike  induced  by  the  rupture  pressure  at  the  end  of  the  driver  section  of  the  shock 
tunnel.  The  duration  of  usable  flow  in  the  shock  tunnel  for  these  tests  is  approximately 
5  milliseconds,  and  the  triggering  circuit  is  adjusted  to  cause  the  laser  to  fire  near 
the  center  of  this  period  (2.7  to  2.8  ms,  typically).  A  15-milliwatt,  continuous-wave 
(CW),  helium  neon  laser  is  mounted  behind  the  rear  reflector  of  the  ruby  laser  cavity 
and  is  aligned  so  that  both  laser  beams  are  coaxial.  The  CW  laser  is  used  to  align  all 
the  components  in  the  recording  system,  and  to  preview  the  holograms  in  situ  prior  to 
analyzing  them  in  detail  using  a  separate  playback  system. 

Figure  26  shows  the  hologram  playback  system. 2^  A  duplication  of  the  holographic 
reference  beam  is  established  by  directing  the  output  of  an  argon  ion,  CW  laser  into 
a  beam  exoander-collimator  similar  to  the  one  used  in  the  recording  system.  The  same 
dual  hologram  holder  is  used  to  position  the  holograms,  and  a  positive,  f/8  lens  is  used 
to  image  the  reconstructed  object  light  waves  directly  into  the  field  of  a  20X  telescope. 
The  20X  telescope  is  the  imaging  component  of  an  x-y  microcomparator.  This  system 
is  used  to  measure  the  interferometric  fringe  shifts  for  both  single  plate  and  dual  plate 
applications,  as  well  as  to  obtain  photographic  enlargements  of  selected  apertures  for 
the  flow  fields.  This  latter  feature  is  especially  valuable,  because  it  allows  details  of 
the  flow,  which  are  otherwise  invisible,  to  be  observed  and  studied.  Alternatively,  a 
Graflex  camera  with  a  4  x  5  Polaroid  film  holder  is  used  in  place  of  the  microcomparator 
to  photograph  the  full-aperture  views  of  the  flow  fields.  During  hologram  playback, 
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the  reconstructed  light  waves  converge  from  the  holograms  to  form  a  duplication  of 
the  true  focus  of  the  second  parabolic  reflector  in  the  recording  system,  permitting 
schlieren  studies  to  be  accomplished  at  leisure. 

With  the  exception  of  one  test  for  the  2-D  compression  ramp  configuration,  dual 
hologram  interferometry25  is  the  primary  optical  process  used  in  this  study. 
Interferograms  of  the  flows  are  obtained  by  interfering  light  waves  produced 
simultaneously  from  a  test  hologram  and  a  reference  state  hologram;  shadowgrams  and 
schlieren  results  are  obtained  directly  from  the  test  holograms.  A  test  hologram  is  a 
recording  of  a  specific  shock  tunnel  test,  while  a  reference  state  hologram  is  a  recording 
of  a  no-flow,  model-out  condition  in  the  test  section  for  either  atmospheric  or  pump- 
down  states.  The  test  section  of  the  shock  tunnel  is  pumped  to  near  vacuum  levels 
for  each  test,  and  the  pump-down  reference  state  hologram  is  used  to  negate  the 
effects  of  stress  in  the  test  section  windows.  A  one-atmosphere  pressure  loading  is 
supported  by  each  window  during  a  test.  For  these  tests,  window  stress  appears  to  have 
a  negligible  effect  on  the  optical  measurements,  because  the  reference  waves  of  both 
types  of  reference  state  holograms  produce  identical  interferograms.  For  this  study, 
the  reference  state  holograms  are  the  lead  holograms  in  the  dual  plate  alignment  so 
that  during  the  hologram  playback  step,  the  reconstructed  reference  waves  are  generated 
first  and  pass  through  the  test  hologram.  This  alignment  allows  the  reconstructed  test 
waves  to  propagate  freely  from  the  hologram  into  the  data  recording/analyzing  system. 
Glass  compensators,  which  would  correct  the  lateral  shift  of  the  reference  light  waves, 
are  not  placed  in  front  of  the  test  holograms  during  the  recording  step,  because  this 
ray  shifting  is  small  enough  to  be  corrected  by  the  final  alignment  of  the  holograms 
during  the  playback  step. 

4.3  MODELS,  INSTRUMENTATION  AND  TEST  CONDITIONS 

The  models  and  instrumentation  used  in  this  investigation,  and  the  conditions  at 
which  the  experiments  are  conducted,  are  selected  on  the  basis  of  the  knowledge 
accumulated  from  a  large  number  of  studies  over  the  past  15  years.  The  flow  fields 
over  the  configurations  selected  are  known  to  exhibit  definitive  features  in  the  separation 
and  reattachment  regions  as  well  as  definitively  formed  separated  regions. 

Studies  of  the  characteristics  of  incident  shock-induced  turbulent  boundary  layer 
separation  have  been  made  at  Mach  numbers  from  6  to  13  (Ref.  26)  with  the  model 
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shown  in  Figure  27.  The  strength  of  the  interaction  is  controlled  by  the  angle  of  the 
shock  generator,  and  the  boundary  thickness  is  varied  by  changing  the  length  of  the 
flat  plate.  The  flat  plate,  which  is  18  inches  wide,  is  extensively  instrumented  with 
heat  transfer  and  pressure  gages  along  the  centerline  and  also  along  selected  spanwise 
rays.  A  shock  generator  angle  of  15  degrees  and  a  flat  plate  length  of  38  inches  ahead 
of  the  interaction  is  the  basic  configuration  used  in  these  tests.  This  configuration  gives 
a  well  separated  region  and  well  defined  heat  transfer  and  pressure  distributions. 

Studies  of  the  flow  over  the  compression  surface  are  conducted  using  the  flat 
plate/wedge  model  shown  in  Figure  28.  This  model,  which  has  been  used  in  a  number 
of  studies  of  laminar  and  turbulent  interacting  flows27>  uses  a  major  part  of  the 
instrumentation  from  the  shock  generator  model  described  above.  A  wedge  angle  of 
36  degrees  is  selected  for  the  current  measurements,  again  to  obtain  a  well  separated 
flow. 


The  large  cone  flare  model  shown  in  Figure  4  is  used  for  measurements  over  an 
axisymmetric  configuration  for  which  the  transverse  curvature  effects  are  small.  A 
9-foot  long,  6°  half-angle  cone  with  a  36°  flare  is  used  in  these  studies.  This 
configuration  is  selected  so  that  the  flow  field  surveys  made  through  the  boundary 
layer  and  separated  region  just  ahead  of  the  flare  can  be  compared  to  the  results  in 
a  complementary  study. 

These  experimental  studies  are  conducted  in  Calspan's  96"  Shock  Tunnel  at  Mach 
numbers  11  and  13  and  at  the  respective  Reynolds  numbers  listed  in  Table  2.  The 
boundary  layers  are  completely  turbulent  well  upstream  of  the  interactions. 
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Table  2 

TEST  CONDITIONS  FOR  EXPERIMENTAL  STUDIES 


Flat  Plate/ 
Wedge 

Incident 

Shock 

Cone/ 

Flare 

Mach  Number 

11.34 

11.34 

13.06 

Reynolds  Number 

1.7017E7 

1.039E7 

5.071E6 

Reservoir  Temperature  (°R) 

2.564E3 

2.585E3 

3.005E3 

Reservoir  Pressure  (psia) 

1.739E4 

1.721E4 

1.836E4 

Freestream  Velocity  (ft/s) 

5.079E3 

5.5735E2 

6.340E3 

Freestream  Density  (slugs/ft^) 

1.662E-4 

1.621E4-4 

6.593E-5 

Freestream  Temperature  (°R) 

1.053E2 

1.064E2 

9.797E1 

Freestream  Pressure  (psia) 

2.085E-1 

2.056E-1 

7.696E-2 

Pitot  Pressure  (psia) 

3.481E1 

3.428E1 

1.709E1 

45 


4.4 


RESULTS  AND  DISCUSSION 


4.4.1  Comparison  to  Flat  Plate  Theory 


To  evaluate  the  accuracy  of  deducing  density  from  the  interferograms,  detailed 
fringe  shift  measurements  are  made  for  the  turbulent  boundary  layer  approximately  4 
inches  upstream  of  the  corner  interaction  (Figure  29).  These  results  are  compared  with 
simple  prediction  methods  for  hypersonic  turbulent  boundary  layers  over  a  flat  plate 
with  a  sharp  leading  edge.  To  maintain  the  greatest  accuracy,  the  measurements  are 
made  by  directly  viewing  the  holographically  reconstructed  image  with  an  x-y 
microcomparator.  These  measurements  are  compared  with  predictions  based  on  the 
Van  Driest  compressibility  transformation  of  the  Clauser/Cole  incompressible  velocity 
profiles  and  a  modified  Crocco  relationship.^ 

4.4.2  Reduction  of  Interferometer  Measurements 


The  reduction  of  interferometric  measurements  to  density  for  a  2-D  flow  field  is 
a  linear  relationship  between  density  and  fringe  shift  which  is  written  as, 

s  .  KgG/  (/’-/’ref)  (1) 


where:  S  =  measured  fringe  shift, 

dimensionless 

kDG  =  Dale  Gladstone  constant 

A  =  wavelength  of  light 

£  -  extent  of  2-D  field  traversed 
by  light  waves 

/?-/?  ref  =  relative  change  in  density  between 
the  reference  point  and  the 
point  where  S  is  measured 
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FLOW  STRUCTURE  IN  A  TWO  DIMENSIONAL  CORNER 


Normally,  -t  is  assumed  to  be  the  geometric  length  of  the  field  in  the  direction  of  the 
light  waves,  but  rarely  is  this  the  case  for  wind  tunnel  applications.  When  the  model 
is  bounded  by  vertical  sidewalls,  the  extent  of  the  2-D  field  is  less  than  the  geometric 
width  of  the  model,  and  the  numerical  value  for,/  in  Equation  (1)  is  reduced  usually  by 
an  estimation  for  the  boundary  layer  displacement  thickness  of  the  sidewall  boundary 
layers.  The  opposite  is  true  for  models  placed  in  the  core  of  a  free  jet,  as  is  the 
case  here.  Because  the  flow  is  unconstrained  in  the  lateral  direction  (no  sidewalls), 
the  density  has  to  be  continuous  around  the  surfaces  of  the  model.  The  model  develops 
boundary  layers  along  its  sides  similar  to  the  boundary  layer  developed  on  the  testing 
surface.  In  the  absence  of  shock  waves,  all  fluid  properties,  as  well  as  the  streamlines 
in  these  boundary  layers,  are  continuous  around  the  edges  of  the  model,  because  the 
boundary  layers  are  connected.  This  condition  makes  the  extent  of  the  2-D  field  wider 
than  the  geometric  length  of  the  model,  and  the  precise  increase  in  £  is  unknown.  In 
tne  present  case,  a  theoretical  maximum  fringe  shift  across  the  boundary  layer  (Smax) 
of  10.2  is  computed  from  Equation  (1)  using  the as  prescribed  by  shock  tunnel 
flow  field  conditions,  but  the  maximum  S  measured  from  the  interferogram  is  10.59. 
The  density  of  the  freestream  and  the  flow  at  the  wall  are  known  from  other  data,  and 
the  optical  constants  are  known  from  physics.  The  additional  0.39  fringe  shift  is 
believed  to  be  due  to  the  influence  of  the  edge  conditions.  To  have  agreement  between 
the  theoretical  and  measured  Smax,  1  has  to  be  0.56  inch  longer  than  the  18  inch  width 
of  the  flat  plate.  If  the  0.56  inch  increase  in  £  is  ignored,  and  the  fringe  data  are 
reduced  using  the  freestream  conditions,  the  interferometric  measurements  predict  Tw 
=  643°R,  which  is  impossible!  At  the  start  of  a  test,  the  model  temperature  is  nearly 
the  same  as  the  ambient  temperature— 531  to  540°R,  typically— and  even  for  extreme 
heat  transfer  rates,  the  model  temperature  remains  constant  for  the  5-ms  duration  of 
the  test.  Hence,  the  additional  0.39  fringe  shifts  are  likely  due  to  edge  effects. 

As  an  initial  treatment,  $  in  Equation  (1)  is  assumed  to  be  0.56-inch  longer  than 
the  plate  width  so  that  there  is  agreement  in  the  S  measurements  from  the  freestream 
to  the  model  surface.  Although  incorrect,  this  simple  choice  for  1  provides  a  starting 
point  to  show  the  general  trend  of  the  density  distribution,  as  well  as  a  starting  point 
from  which  the  influence  of  the  edge  effects  might  be  determined;  more  realistic 
corrective  methods  are  being  studied  as  reduction  and  analyses  of  these  data  continues. 
For  now,  with  both  Smax  and^^y^specified,  the  quantity  (Kj3Q//\)/in  Equation  (1) 
is  known,  and,  hence,  the  density  at  all  measured  S  locations  is 
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EquaJon  (2)  is  used  to  compute  the  experimental  density  distribution  shown  in  Figure  30. 
For  a  steady  flor,,  the  distribution  would  be  a  reasonably  smooth,  monotonic  function. 
The  mall  anomalies  in  the  actual  distribution  are  believed  to  result  from  the  instanteous 
measurements  of  the  unsteadiness  in  the  flow. 


The  comparisons  between  theoretical  and  experimental  data  are  shown  in 
Figure  31.  The  difficulty  in  the  accurate  reduction  of  the  measurements  close  to  the 
wall  is  reflected  in  densities  that  are  larger  than  the  wall  value  just  above  the  flat 
plate.  For  a  constant  pressure  boundary  layer,  this  would  result  in  temperatures  at 
the  base  of  the  boundary  layer  which  are  lower  than  the  wall  temperature,  which  would 
result  in  heat  transfer  from  the  wall  to  the  flow.  This  result  is  impossible  for  this 
case,  and  the  heat  transfer  measurements  show  this  is  clearly  not  true. 


Using  the  Van  Driest  compressibility  transformation 
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together  with  Coles  modification  of  the  Clauser  incompressible  relationship 
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where  Ct  =  0  iGr  a  linear  relationship  and  Ct  =  0.5  for  a  parabolic  form.  Close  to 
the  wall,  the  measurements  exhibit  the  anticipated  lower  values  associated  with  a  local 
peak  in  temperature  close  to  the  wall. 

The  distribution  of  velocity  through  the  boundary  layer  can  be  calculated  from 
the  measured  density  distribution  by  using  the  modified  Crocco  relationship  (Equation 
4).  These  calculations  for  Ct  =  0  and  0.5  together  with  the  Van  Driest  transformed 
Clauser/Cole  relationship  are  shown  in  Figure  31.  The  velocities  deduced  from  the 
measurements  are  significantly  larger  than  the  predicted  values  and  those  deduced  from 
pitot  and  total  temperature  measurements  obtained  in  an  earlier  Calspan  study.^8  This 
result  is  consistent  with  the  density  discrepancies. 

Clearly,  a  closer  examination  of  refraction  and  edge  effects  as  well  as  the  data 
reduction  technique,  are  required.  In  retrospect,  a  better  approach  may  be  to  examine 
an  axisymmetric  laminar  flow  to  first  evaluate  the  problems  with  the  basic  system. 
Following  this,  the  problems  associated  with  the  intrinsic  three-dimensional  structure 
of  turbulent  flow  and  the  intricacies  of  the  flow  at  the  edges  of  the  two-dimensional 
model  may  be  easier  to  resolve.  Also,  making  independent  measurements  of  the  density 
in  both  laminar  and  turbulent  flows  using  an  electron  beam  are  desirable,  because  this 
would  provide  redundant  and  independent  measurements.  We  intend  to  pursue  this 
approach  to  further  study  these  flows. 

4.4.3  Discussion  of  Qualitative  Features 


A  single  plate  holographic  interferogram  typical  of  those  taken  during  the  current 
study  is  shown  in  Figure  29.  This  particular  hologram  is  made  by  double-exposing  a 
single  plate,  first  with  the  reference  beam  just  before  the  test,  and  then  approximately 
15  minutes  later  with  the  scene  beam  of  the  test.  Here,  the  interferogram  is  fixed  in 
the  hologram.  In  this  interferogram,  the  model  wall  appears  as  a  solid  black  surface, 
while  in  the  interferogram  reconstructed  using  the  dual  plate  technique,  the  model 
appears  as  a  gray  silhouette.  The  photographs  shown  in  this  oaper  are  enlargements 
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of  35mm  photographs  of  the  reconstructed  object  waves  as  viewed  through  the  eyepiece 
of  a  20X  telescope.  Because  this  direct  viewing  approach  is  used  in  this  study,  details 
in  the  hologram  can  be  magnified  to  achieve  pictures  of  greater  clarity  than  would 
otherwise  i.e  obtained.  Restating  the  fact  that  all  optical  data  shown  here  result  from 
the  integrated  effects  of  density  changes  from  the  beam  splitter  to  the  recording  plane, 
these  photographs  must  be  interpreted  with  care,  especially  in  and  close  to  regions  of 
separated  flow. 

4.4.4  Turbulent  Separation  in  a  Two-Dimensional  Compression  Corner 

The  scale  and  properties  of  separated  regions  formed  over  flat  plate/wedge 
compression  surfaces  are  controlled  by  a  "free"  interaction  between  the  viscous  layer 
and  the  overlaying  inviscid  flow.  Because  of  the  intrinsic  differences  between  the 
structure  and  properties  of  laminar  and  turbulent  boundary  layers,  the  mechanism  of 
flow  separation  and  the  stiucture  of  the  separated  flows  are  significantly  different.  In 
laminar  flows,  the  major  mechanism  controlling  the  flow  in  regions  of  separation  is 
the  mutual  interaction  between  the  growth  of  the  displacement  layer  and  the  outer 
inviscid  flows  that  can  be  described  in  terms  of  second  order  boundary  layer  theory. 

Turbulent  separation  takes  place  in  the  sublayer  from  which  the  strong  separation 
shock  traverses  the  boundary  resulting  in  strong  normal  pressure  gradients.  Such  flows 
cannot  be  described  by  boundary  layer  theory  and  possibly  not  by  the  mass  averaged 
Navier-Stokes  equations.  A  schlieren  photograph  and  distributions  of  heat  transfer  and 
pressure  (Figure  32)  are  shown  for  separated  flow  over  the  flat  plate/wedge 
configuration. The  corresponding  holographic  interferograms,  schlieren  photographs 
and  shadowgraphs  that  are  obtained  in  the  current  study  are  shown  in  Figures  33  and 
34.  The  heat  transfer  and  pressure  distributions  clearly  indicate  that  the  separation 
region  occupies  a  fraction  of  a  boundary  layer  thickness.  The  separation  shocks  that 
are  created  in  this  region  are  seen  clearly  in  the  shadowgraph  and  schlieren  photographs, 
and  they  appear  as  rapid  distortions  of  the  fringe  pattern  in  the  interferogram.  These 
shocks  are  highly  distorted  as  they  traverse  the  major  part  of  the  boundary  layer.  A 
matter  of  continuing  debate  is  whether  these  distortions  are  responses  to  the  unsteady 
movement  of  the  separation  point  or  to  responses  in  the  incoming  turbulent  flow. 
However,  it  is  clear  that  significiant  pressure,  temperature,  and,  hence,  density 
fluctuations  are  generated  by  the  inflections  of  the  separation  shock  in  this  region. 
The  multiple  shocks  that  are  typically  observed  in  the  separation  region  reflect  the  time- 
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dependent  position  of  separation  (since  they  ate  also  observed  on  axisymmetric  models) 
and  the  three-dimensional  nature  of  the  flow  at  the  edges  of  the  model. 

A  prominent  feature  of  the  interferogram  is  a  relatively  uniform  high  temperature 
structure  of  the  constant  pressure  recirculation  region  (Figure  29).  The  straight  shear 
layer  can  be  seen  to  be  lifted  above  the  recirculating  flow,  which  appears  to  be 
constructed  from  a  single  vortex  sheet  with  some  smaller  cells  imbedded  close  to  the 
reattachment  region.  This  feature  is  significant  from  the  viewpoint  of  the  turbulence 
modeling  of  these  flows.  The  rapid  recompression  that  is  always  observed  in 
reattachment  regions  is  accompanied  by  the  formation  of  a  reattachment  shock.  This 
shock  relaxes  to  the  wedge  shock  as  the  influence  of  the  more  efficient  compression 
process  close  to  the  separated  region  disappears.  Meaningful  interpretation  of  the 
interferometer  data  obtained  is  difficult,  if  not  impossible,  because  significant  refraction 
effects  are  clearly  observed  in  the  interferogram.  As  observed  and  discussed  in  earlier 
papers,  the  schlieren  photographs  reveal  an  almost  sinusoidal  shock  structure  in  the 
reattachment  region,  possibly  reflecting  a  movement  of  the  reattachment  point.  Such 
a  movement  indicates  changes  in  the  flow  reversal  in  the  separated  region.  This  aspect 
of  flow  unsteadiness  could  significantly  influence  the  size  of  the  separated  region,  and 
clearly  complicates  a  numerical  description  of  this  flow.  The  development  of  turbulence 
through  the  large  pressure  gradients  in  the  reattachment  region  and  the  subsequent 
relaxation  in  the  constant  pressure  region  downstream  of  reattachment  are  important 
features  of  this  flow.  If  refraction  and  edge  effects  can  be  eliminated  from  the  fringe 
shift  measurement,  holographic  interferometry  will  provide  extremely  valuable  quantative 
measurements  of  this  separated  flow  field. 

The  reattachment  region  presents  additional  difficulties,  because  the  strong  flow 
curvature  and  rapid  thinning  of  the  boundary  layer  here  may  be  generating  three- 
dimensional  flow  structures  (Goertler  vortices)  that  clearly  prevent  a  meaningful 
interpretation  of  these  interferograms.  Such  three-dimensional  effects  also  may  be 
induced  in  the  strong  flow  curvature  close  to  the  separation  point.  Basically,  the 
unsteady  and  three-dimensional  characteristics  of  these  flows  must  be  understood  to 
develop  accurate  numerical  descriptions  of  these  flows,  but  an  understanding  of  the 
unsteady  and  three-dimensional  effects  is  also  needed  to  help  interpret  the  optical 
measurements.  This  point  further  suggests  a  need  to  study  lesc  complicated  flows  first 
(e.g.,  axisymmetric  laminar  cases  previously  mentioned). 
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4.4.5  Incident  Shock  Wave/Turbulent  Boundary  Layer  Interaction 


As  in  the  flow  over  the  flat  plate/wedge  compression  corner,  the  size  and  flow 
structure  in  a  region  of  incident  shock  impingement  is  controlled  by  the  mutual 
interaction  between  the  viscous  and  inviscid  flow.  However,  the  disturbances  are 
generated  outside  the  boundary  layer,  and  the  incident  shock  strength  is  modified  by 
its  interaction  with  the  separation  shock.  The  measurements  of  heat  transfer  and 
pressure  on  this  configuration  are  shown  in  Figure  35,  together  with  schlieren  photographs 
made  in  earlier  studies. 27  The  measurements  made  in  the  current  studies  are  shown 
in  Figures  37  through  39.  The  structure  of  the  upstream  disturbance  generated  by  the 
incident  shock  is  very  similar  to  that  generated  over  the  flat  plate/wedge  model.  The 
basic  viscous/inviscid  interaction  resulting  in  flow  separation  takes  place  at  the  base 
of  the  boundary  layer.  Separation  shock(s)  traverse  through  a  major  segment  of  the 
boundary  layer,  interact  with  the  incident  shock,  and  are  then  turned  almost  parallel 
to  the  plate  surface;  an  expansion  fan  is  generated  as  the  incident  shock  interacts  with 
the  constant  pressure  free  shear  layer.  The  boundary  layer  flow  is  again  turned  parallel 
to  the  flat  plate  by  the  formation  of  the  reattachment  shock.  These  features  are 
clearly  visible  in  the  schlieren  photographs  and  int  ograms  shown  in  Figures  36 
through  38. 

The  structure  of  the  recirculating  region  is  shown  in  detail  in  Figure  36,  which 
again  suggests  a  basic  single  cell  structure  with  a  relatively  uniform  high  temperature 
core.  Again,  however,  the  very  laige  density  gradients  generated  in  the  reattachment 
region  have  caused  proportionately  large  refraction  effects  making  the  interpretation 
of  the  interferograms  virtually  impossible  for  this  region.  Similar  to  the  reattachment 
region  developed  on  the  compression  ramps,  the  strong  flow  curvature  would  likely 
cause  the  formation  of  three  dimensional  disturbances  (Goertler  vortices)  in  this  region 
downstream  of  the  incident  shock.  Again,  the  highly  curved  inflective  shocks  in  the 
separation  and  reattachment  regions  are  indicative  of  an  unsteadiness  which  needs  to 
be  investigated  with  high  frequency  flow  measurements  to  gain  an  understanding  of 
these  unsteady  effects. 

4.4.6  Separated  Flows  at  the  Cone/Flare  Junction 

Interferograms  and  shadowgrams  for  the  axisymmetric  flow  fields  obtained  in 
the  current  study  are  shown  in  Figures  40  and  41.  Measurements  of  the  distribution 
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Figure  35  STREAMWISE  DU  HBUTIONS  OF  HEAT  TRANSFER  AND  PRESSURE 

THROUGH  SKEWED-OBLiQUE-SHOCK  BOUNDARY  LAYER  INTERACTION 
(  0  =  150  \J/  =  QO) 
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Figure  36  HOLOGRAPHIC  DEPICTION  OF  SHOCK  INDUCED  BOUNDARY  LAYER 
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of  heat  transfer  and  pressure  over  this  configuration  made  in  previous  studies^  are 
shown  in  Figure  42  along  with  a  conventional  schlieren  photograph.  The  36°  flare  angle 
generates  a  large  separated  region  which  extends  2  inches  ahead  of  the  cone-flare 
junction.  As  expected,  the  flow  in  the  separated  region  is  very  similar  to  the 

corresponding  regions  on  the  compression  ramp  and  incident  shock  configurations. 
However,  interpretation  of  the  phenomena  causing  the  shock  waves  in  the  neighborhood 
of  the  constant  pressure  separated  region  is  not  understood.  This  multiplicity  of  shocks 
is  also  clearly  discernible  in  the  reattachment  region  of  the  flow,  again  leading  to 
speculation  on  the  intrinsic  unsteady  or  three-dimensional  character  of  the  flow  in  this 

region.  Even  though  the  optical  path  lengths  of  these  reattachment  regions  are 

significantly  shorter  than  those  of  the  two-dimensional  models,  strong  refraction  appears 
to  distort  these  interferograms  as  well.  Finding  the  best  way  to  reduce  the 

interferometer  measurements  of  these  axisymmetric  separated  regions  remains  an 

important  task.  However,  resolving  this  problem  may  be  easier  and  more  productive 
than  attempting  to  evaluate  edge  effects  of  the  two-dimensional  models. 

4.5  CONCLUSIONS 

A  preliminary  investigation  has  been  conducted  into  the  use  of  holographic 
interferometry  to  study  regions  of  shock  wave/boundary  layer  interaction  in  a  hypersonic 
shock  tunnel.  The  interferograms  of  the  complex  flow  fields  generated  in  this 

investigation  provide  a  good  qualitative  basis  for  evaluating  some  of  the  important 
phenomena  that  control  the  characteristics  of  th  ,e  flows.  However,  the  quantitative 
evaluations  of  the  interferograms  are  made  difficult  by  edge  effects  and  three- 

dimensional  flow  effects  on  two-dimensional  models  and  by  refraction  effects  in  regions 
of  high  density  gradients  close  to  the  wall.  The  disagreement  between  density 

measurements  and  simple  prediction  techniques  for  turbulent  flat  plate  boundary  layers 
indicates  that  further  analyses  and  experiments  are  required  to  develop  an  understanding 
of  both  the  refraction  and  three-dimensional  effects.  Future  studies  should  be  directed 
towards  laminar  flows  over  axisymmetric  and  two-dimensional  models  after  which  the 
more  difficult  conditions  of  these  turbulent  flows  can  be  understood. 
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Figure  42  DISTRIBUTION  OF  PRESSURE  AND  HEAT  TRANSFER  IN 

SEPARATED  FLOW  OVER  THE  LARGE  6°  CONE/360  FLARE 
CONFIGURATION 
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Section  5 
CONCLUSIONS 


This  report  summarizes  the  research  conducted  under  the  current  contract  to 
investigate  fundamental  problems  associated  with  flight  at  hypersonic  speeds  with 
particular  emphasis  on  those  related  to  the  aerothermal  effects  of  viscous/inviscid 
interactions  and  boundary  layer  transition. 

Solutions  have  been  obtained  to  Navier-Stokes  equations  for  the  laminar  flow 
over  the  leading  edge  of  a  sharp  flat  plate  in  Mach  16  flow  for  highly  cooled  wall 
conditions  using  a  modified  MacCormack/Shang  fully  explicit  formulation.  To  obtain  a 
stable  converged  solution,  it  was  necessary  to  reduce  the  grid  size  close  to  the  leading 
edge  to  the  order  of  the  mean  free  path,  and,  typically,  20,000  time  steps  were  required 
to  achieve  convergence.  However,  once  obtained,  the  solution  was  in  good  agreement 
with  the  experiment. 

The  experimental  study  conducted  under  the  current  contract  was  directed  toward 
obtaining  detailed  flow  field  measurements  in  a  separating  boundary  layer  over  a  large 
cone/flare  model.  In  this  study,  the  structure  of  the  separating  boundary  layer  was 
examined  with  pitot,  total  temperature  and  laser  holography  measurements.  These 
measurements  suggest  that  the  total  temperature/velocity  relationship  in  the  cone 
boundary  layer  is  quadratic  rather  than  the  usually  assumed  linear  form  suggested  by 
Crocco.  The  Van  Driest  transformation,  which  has  been  used  successfully  in  supersonic 
flows  ovei  adiaoatic  walls  to  relate  the  measured  velocity  to  similar  measurements  in 
subsonic  flow,  is  apparently  not  as  effective  in  hypersonic  flows  over  highly  cooled 
walls.  The  measurements  in  the  separation  region  demonstrate  a  rapid  change  in  the 
structure  of  the  sublayer  as  separation  takes  place  and  the  formation  of  strong  shock 
waves  in  the  turbulent  shear  layer. 

In  addition  to  experimental  studies  of  turbulent  boundary  layer  separation  in 
hypersonic  flow,  during  the  past  year,  further  analysis  has  been  performed  of 
measurements  made  during  earlier  studies  for  AFOSR.  This  work  has  resulted  in  two 
AIAA  publications  during  the  course  of  this  current  contract.  We  investigated  the 
aerothermal  characteristics  of  nosetips  and  leading  edges  as  transition  moves  onto  the 
nosetip.  In  these  studies,  it  was  shown  that  heat  transfer  rates  significantly  larger 
than  the  stagnation  value  can  be  generated  near  the  sonic  line  close  to  the  end  of  the 
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transition  region,  and  such  large  heat  transfer  rates  can  cause  serious  indentations  of 
an  ablative  nosetip.  The  results  of  these  studies  for  transitional  flows  over  hemispherical 
and  indented  noseshapes  were  reported  at  the  annual  AIAA  meeting  in  January  1986 
(see  Appendix  A).  Finally,  an  invited  paper  reviewing  the  status  of  research  in 
aerothermal  problems  in  hypersonic  flight  with  emphasis  on  those  associated  with 
boundary  layer  transition,  and  viscous/interaction  and  flow  separation,  was  presented 
at  the  Annual  AIAA  Meeting  in  January  1986  (Appendix  B)- 
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ABSTRACT 

Experimental  studies  have  been  conducted  to 
examine  the  flow  structure  and  distribution  of  heat 
transfer  and  pressure  over  indented  noseshapes. 
Measurements  were  made  to  examine  the  influence  of 
boundary  layer  trips;  surface  roughness;  and  Mach  number 
and  Reynolds  number  of  the  freestream,  on  the  size  and 
structure  of  the  shock/turbulent  boundary  interaction  over 
the  nosetips.  It  was  demonstrated  that  Mach  number  and 
surface  roughness  are  of  key  importance  in  the  correct 
simulation  of  these  flows,  and  it  is  essential  in  the  correct 
simulation  of  these  flows,  and  it  is  essential  that  the 
boundary  be  fully  turbulent  upstream  of  the  body 
shock/boundary  layer  interaction.  In  similar  studies  of 
the  flow  over  the  recovered  NRV  noseshape  at  conditions 
which  duplicated  the  flight  Mach  number  and  Reynolds 
number,  it  was  demonstrated  the  flow  was  intrinsically 
three-dimens.onal  and  that  the  distribution  of  heat  transfer 
and  pressure  were  also  controlled  by  boundary  layer 
transition  and  shock/boundary  layer  interaction.  In  a  study 
performed  specifically  for  code  validation,  detailed 
measurements  of  heat  transfer  and  pressure  were  made 
in  laminar  flows  over  two  indented  noseshapes.  Navier- 
Stokes  solutions  for  one  of  these  case  obtained  by  Widhopf 
and  Victoria  exhibited  unsteady  flows.  Experiment  showed 
the  flow  to  be  stable;  however,  it  was  observed  that  the 
heat  transfer  measurements  were  in  agreement  with  the 
time  averaged  predictions. 

I.  INTRODUCTION 

Wind  tunnel  studies  using  ablating  models  have 
demonstrated  that  "slender"  concave  nose  shapes,  over 
which  separated  regions  and  gross  flow  instabilities  are 
observed,  can  be  formed  under  certain  free  stream 
conditions  and  nosetip  properties.  While  an  explanation 
of  this  behavior  using  sophisticated  shape  change  codes 
awaits  the  development  of  accurate  models  of  regions  of 
transition  and  shock  wave/turbulent  boundary  layer 
interaction  the  simple  equilibrium  solutions  obtained  by 
Welsh*  provide  an  important  insight  into  the  development 
of  ablated  nose  shapes.  Welsh's  analysis  predicts  the 
occurrence  of  the  well  known  "blunt"  self-preserving  nose 
shape  for  fully  laminar  or  fully  turbulent  boundary  layers; 
however,  for  flows  where  transition  occurs  ahead  of  the 
sonic  point  on  the  nosetip  (and  under  certain  fully  turbulent 
conditions)  this  approach  predicts  the  occurrence  of 
slender  nose  shapes  similar  to  those  shown  in  Figure  1. 
This  result  has  been  verified  in  wind  tunnel  tests  of 
camphor  models  where  shapes  similar  to  those  shown  on 
Figure  2  can  be  developed  under  transitional  conditions. 
The  existence  of  slender  ablated  nose  shapes  have  been 
demonstrated  in  many  experimental  studies2’3’1*1'  in 
general,  however,  these  shapes  exhibit  marked  asymmetric 
characteristics.  Asymmetric  nose  shapes  can  give  rise  to 
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large  de-stabilizing  moments  and  side  forces.  Highly 
indented  nose  shapes,  under  conditions  similar  to  that 
shown  in  Figure  3  were  generated  on  camphor  models, 
and  over  such  shapes  the  viscous  and  inviscid  flow  fields 
interact  in  such  a  way  that  the  flow  became  oscillatory. 

In  studies  of  ablated  nose  shapes,  it  has  become 
customary  to  describe  the  nosetip  geometry  in  terms  of 
the  BI-CONIC  or  TRI-RADIUS  configuration.  Schematic 
representations  of  these  shapes  are  shown  in  Figure  4a. 
The  form  of  a  highly  ablated  nose  shape  which  is  to  be 
characterized  by  these  models  is  shown  in  Figure  5.  The 
increased  heating  resulting  from  boundary  layer  transition 
close  to  the  sonic  region  is  believed  responsible  for  the 
initial  development  of  the  identation  in  ablating  nose 
shapes.  Because  transition  in  such  a  situation  is  almost 
certainly  three-dimensional  in  character,  the  initial 
indenting  of  the  nosetip  will  also  reflect  the  asymmetric 
characteristics  of  the  transition  process.  Subsequently 
the  large  heating  rates  which  are  developed  in  the 
recompression  region  resulting  from  shock-boundary  layer 
interaction  at  the  base  of  the  indentation  can  cause  grossly 
indented  nose  i  apes.  While  it  is  possible  to  predict  .he 
initial  ...J».,ting  using  relatively  simple  shape  change  code, 
when  significant  ablation  occurs  the  computation  of 
heating  over  the  rough  surface  for  flows  with  embedded 
regions  of  shock  wave-boundary  layer  interaction  are 
subject  to  considerable  uncertainties.  To  calculate  the 
effects  of  surface  roughness  on  heating  and  surface 
recession,  the  magnitude  of  the  surface  roughness  and  the 
roughness  augmented  heating  in  the  presence  of  blowing 
must  be  predicted.  The  prediction  of  the  heating  to  the 
embedded  regions  of  shock  wave  turbulent  boundary  layer 
interaction,  some  involving  separated  flows,  over  a  rough 
surface  represents  one  of  the  most  difficult  configurations 
to  calculate  with  precision.  Although  there  have  been  a 
number  of  experimental  studies^i^  to  examine  the 
distribution  of  heating  to  non-indented  nose  shapes  and 
investigations  with  a  major  emphasis  on  examining  the 
gross  stability  of  flows  over  highly  indented  nose  shapes, 
little  detailed  information  is  available  for  indented 
nosetips  where  attached  or  slightly  separated  regions  of 
shock  wave-turbulent  interactions  occur.  Such  studies  are 
required  to  provide  a  method  of  evaluating  and  developing 
the  prediction  methods,  and  obtaining  better  insights  into 
the  basic  phenomena  leading  to  nosetip  indentation. 

This  paper  is  composed  of  three  experimental  studies 
designed  to  provide  fundamental  information  on  the 
characteristics  of  the  flow  over  indented  nose  shapes.  In 
the  first  study  ws  examined  the  parameters  which  control 
the  correct  simulation  of  the  flow  over  a  non-ablating 
axisymmetric  indented  nosetip.  Detailed  measurement  of 
heat  transfer,  pressure  and  flow  field  structure  were  made 
on  rough  and  smooth  configurations  for  a  range  of  Mach 
numctrs  and  Reynolds  numbers.  We  describe  a  study  using 
a  4X  replica  of  the  recovered  NRV  nosetip,  in  which 


A-5 


detailed  heat  transfer,  pressure  and  flow  field 
measurements  were  made  for  Mach  numbers  and  Reynolds 
number  conditions  which  duplicated  those  obtained  in 
flight.  Finally  the  results  of  a  study  of  lai....w  flows 
over  two  indented  nose  shapes  are  reported  to  provide  a 
basis  for  evaluating  solutions  obtained  using  the  full  Navier 
Stokes  equations. 

II.  STUDIES  TO  INVESTIGATE  THE  SIMULATION  OF 

FLOWS  OVER  INDENTED  NOSE  SHAPES 

Introduction 

Because  the  initiation  and  development  of  nosetip 
indentation  is  so  strongly  influenced  by  the  occurrence 
and  position  of  boundary  layer  transition  on  the  nosetip, 
experimental  studies  to  simulate  these  flows  must  be  very 
carefully  constructed.  To  correctly  simulate  the  flow 
over  a  mildly  indented  nose  shape,  it  is  of  fundamental 
importance  that  the  Reynolds  number  be  of  a  sufficient 
magnitude  to  move  transition  relatively  close  to  the 
stagnation  region  and  preserve  turbulent  flow  through  the 
strong  expansion  which  precedes  the  model  concavity.  To 
correctly  simulate  the  flow  field  and  the  distribution  of 
heating  beneath  the  embedded  region  of  shock  wave/ 
turbulent  boundary  layer  interaction,  the  Mach  number, 
Reynolds  number  and  wall-to-free  stream  stagnation 
temperature  ratio  must  be  correctly  simulated  as  well  as 
the  surface  roughness  of  the  model. 

In  this  segment  of  the  study  we  examined  the  flow 
field  and  distribution  of  heating  to  the  indented  nosetip 
at  high  Mach  numbers  for  unit  Reynolds  numbers  large 
enough  to  establish  turbulent  boundary  layers  over  the 
model  without  tripping.  We  wished  to  investigate  the 
effects  of  trips  and  free  stream  Mach  number  on  the 
heating  distr.bution  for  a  fully  turbulent  condition.  We 
then  wanted  to  determine  the  effects  of  surface  roughness 
or.  the  distribution  of  heating  and  flow  field  properties 
for  both  high  and  low  Reynolds  number  free  stream 
conditions. 

Models  and  Instrumentation 

The  model  used  In  this  study  was  as  tri-radius 
representation  (see  Figure  4b)  of  an  ablated  indented 
nosetip.  This  configuration  is  one  of  a  family  of  nose 
shapes,  a  number  of  which  have  been  examined 
experimentally  in  earlier  studies  conducted  at  Calspan  to 
examine  the  stability  boundaries  for  the  flow  over  indented 
nosetips.  From  the  earlier  studies^  conducted  at  Calspan, 
we  would  predict  that  provided  the  boundary  layer 
remained  turbulent  over  the  nosetip  ahead  of  the  concave 
section  of  the  nosetip,  the  boundary  layer  should  remain 
attached  or  slightly  separated  from  the  base  of  the  body 
shock.  Our  earlier  studies  with  rough  and  smooth  45 
degree  biconic  nosetips4  suggested  that  a  fully  turbulent 
boundary  layer  could  be  anticipated  for  free  stream 
Reynolds  numbers  of  10  x  lOfyft  and  5  x  10^/ft  at  Mach 
11  and  13,  respectively,  while  for  Reynolds  numbers  of  2 
x  10^/ft,  the  upstream  of  the  body  shock  would  be  laminar 
or  transition.  In  this  latter  case,  we  would  expect  laminar 
separation  to  occur  and  spread  forward  to  the  tip  causing 
the  flow  field  to  differ  markedly  from  the  equivalent  fully 
turbulent  case. 

The  mode!  of  the  indented  nosetip  used  in  these 
experimental  studies  is  shown  in  both  its  smooth  and  rough 
configuration  in  Figure  6.  The  rough  model  had  15  mil 
sand  grain  roughness  bonded  to  its  surface  with  pressure- 
sensitive  adhesive  and  was  instrumented  with  rough 
calorimeter  heat  transfer  gages  (see  Ref.  6).  The  gage 


positions  were  selected  to  provide  detailed  coverage  in 
the  area  of  greatest  interest— the  recompres- 
sion/reattachment  region  beneath  the  body  shock.  The 
gradients  of  heating  in  this  area  are  such  that  isolated 
gages,  or  gages  placed  on  a  low  diffusivity  substrate  are 
required  to  prevent  inaccuracies  which  will  result  from 
streamwise  heat  conduction.  Employing  "thin  skin"  heat 
transfer  techniques  with  a  non-segmented  model  could 
lead  to  serious  inaccuracies  in  the  heat  transfer 
measurements  in  the  shock  interaction  region  on  the  model 
(Ref.  7).  We  used  a  single  pass  schlieren  system  with  a 
focal  length  of  10  feet  for  flow  visualization  in  this  study., 
the  knife  edge  was  set  parallel  to  the  conical  segment 
which  precedes  the  concave  section  of  the  model,  and 
adjusted  to  give  a  15  percent  cut-off. 

Results  and  Discussion 

In  this  study  the  effects  of  the  Mach  and  Reynolds 
number  of  the  freestream  were  investigated  first  on  the 
smooth  indented  nose  shapes.  We  then  coated  same  nosetip 
with  15  mu  roughness  and  again  examined  the  effect  of 
Mach  number  and  Reynolds  number.  Heat  transfer  and 
schlieren  measurements  were  made  on  both  configurations 
to  examine  the  flow  mechanisms  upstream  and  in  the 
regions  of  body  shock-boundary  layer  interaction  over  the 
nosetip.  Particular  emphasis  was  placed  in  these  studies 
on  examining  how  the  scale  of  the  body  shock-boundary 
layer  interaction  and  the  distribution  and  peak  heating  is 
influenced  by  the  freestream  Mach  number  and  the  ratio 
of  roughness  height  to  momentum  thickness  (K 19)  upstream 
of  the  interaction  region. 

Studies  with  the  Smooth  Nose  shapes 

The  first  set  of  measurements  were  made  at  Mach  1 1 
for  a  Reynolds  number  (Regy.)  of  12  x  106  based  on  booy 
diameter.  Here  the  emphasis  was  first  on  determining 
whether  we  could  obtain  a  fully  turbulent  boundary  layer 
over  the  model  without  employing  boundary  layer  trips. 
We  then  examined  the  effects  of  boundary  layer  trips  on 
the  properties  of  the  body  shock  interaction  region.  The 
effects  of  Mach  number  and  Reynolds  number  were  then 
studied. 

The  distribution  of  heat  transfer  to  the  completely 
smooth  model  and  a  schderen  photograph  of  the  flow  at 
M=ll  Red=10^,  are  shown  in  Figures  7  and  8.  The  heat 
transfer  measurements  demonstrate  that  the  boundary 
layer  remains  turbulent  as  it  expands  onto  the  short  conical 
section  of  the  model  upstream  of  the  body  shock-boundary 
layer  interaction.  The  severity  of  the  shock  interaction 
is  such  that  a  small  bulent  separated  region  is  induced 
ahead  of  the  body  shock,  its  presence  being  marked  by 
separation  shocks  in  the  schlieren  photographs  and 
indicated  by  the  characteristic  rise  and  plateau  region  ,n 
heat  transfer  distribution  upstream  of  the  body  shoex.  A 
steep  rise  in  the  heating  rate  occurs  in  the  re-compression 
region  downstream  of  the  body  shock  resulting  ,n  peak 
heating  rates  of  over  three  times  those  measured  at  the 
stagnation  point.  The  inviscid  flow  above  the  boundary 
layer  downstream  of  the  body  shock  is  dominated  by  the 
interaction  between  the  nosetip,  body  and  separation  and 
re-utachment  shocks.  The  vorticity  generated  in  the 
inviscid  flow  behind  the  body  shock  is  reflected  in  the 
large  density  variations  which  can  be  observed  in  the 
schlieren  photographs.  The  complexity  of  interacting  flow 
m  this  important  region  )ust  downstream  of  the  body  shock 
makes  a  numerical  analysis  of  these  flows  very  difficult. 

We  then  examined  whether  employing  an  annular  15 
mil  trip  ring  around  the  sonic  point  on  the  nosetip  would 
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disturb  the  flow  structure  and  distribution  of  heating  to 
the  model  under  the  fully  turbulent  conditions  found  in 
Run  L.  The  schliercn  photographs  and  distribution  of 
heating  for  this  tripping  experiment  are  shown  in  Figure  S. 
There  is  little  difference  between  the  flow  patterns 
obtained  in  Runs  1  and  2;  and  the  comparisons  between 
the  measured  heating  distributions  shown  in  Figure  7 
indicates  that  with  the  exception  of  the  region 
immediately  ahead  and  behind  the  trips,  the  distribution 
of  heating  is  relatively  uninfluenced  by  the  presence  of 
the  trips.  The  effects  of  freestream  Mach  number  were 
explored  by  repeating  the  studies  at  Mach  13.  The 
schlieren  photographs  and  distribution  of  heating  obtained 
are  shown  in  Figures  9  and  10.  Both  the  schlieren 
photographs  and  the  heat  transfer  distribution  indicate 
that  the  flow  remains  attached  upstream  of  the  body 
shock.  Comparing  the  distribution  of  heat  transfer 
coefficients  for  the  Mach  11  and  13  flows  (shown  in 
Figure  9),  it  is  clear  that  the  peak  heat  transfer 
coefficients  for  the  Mach  13  flow  are  almost  twice  those 
at  Mach  11.  We  believe  this  results  from  a  greater 
severity  of  the  compression  process  in  the  interaction 
region  of  the  Mach  13  flow.  It  is  clear  from  the  levels 
of  heating  through  the  body  shock-boundary  layer 
interaction  region  that  the  boundary  layer  was  fully 
turbulent  in  both  the  Mach  11  and  13  flows  over  the 


measurements  by  using  isolated  gages  or  a  poorly 
conducting  model. 


STUDIES  OF  THE  FLOW  OVER  RECOVERED  NOSE 
SHAPES 


Objective 


The  Nose  Tip  Recovery  Vehicle  (NRV)  and  the  Role 
Torque  Experiment  Vehicle  (RTE)  nosetips  are  two  of  the 
few  nosetips  which  have  been  recovered  from  flight  during 
re-entry®.  These  nosetips  are  of  particuia,  interest 
because  they  were  "caught"  during  the  portion  of  the 
trajectory  where  transition  was  spreading  over  the  noscti^, 
and  a  non-similar  shape  change  was  occurring.  The 
nosetips  shown  in  Figure  15,  display  two  characteristics 
which  are  of  key  importance  to  the  nosetip  designer.  The 
first  is  that  the  nosetip  is  indented  to  the  point  where 
three-dimensional  separated  regions  of  significant 
proportions  must  have  been  formed.  The  second  feature 
is  *h'  distinctive  three-dimensional  grooved  shape  of  the 
NRV  and  RTE  nosetips,  shapes  which  could  not  have  been 
predicted  from  shape  change  codes. 


nosetip. 

Studies  with  the  Rough  Nosetip 

The  distribution  of  heating  to  and  the  flow  pattern 
of  the  rough  nosetip  at  Mach  11  for  RE^pll  x  10”  are 
shown  in  Figures  11  and  12.  While  the  position  of  the 
body  shock  remains  essentially  unchanged  by  surface 
roughness,  separation  has  spread  further  forward  than  on 
the  smooth  model  at  the  same  freestream  conditions.  We 
believe  this  forward  moment  of  separation  results 
principally  from  a  significant  thickening  of  the  boundary 
layer  upstream  of  the  interaction  combined  with  the 
decreased  velocity  adjacent  to  the  wall  both  of  which 
result  from  surface  roughness.  The  heat  transfer 
measurements  on  the  rough  and  smooth  configurations, 
which  are  compared  in  Figure  11,  show  that  while  the 
position  of  the  region  of  peak  heating  remains  unchanged 
by  surface  roughness,  the  recompression  heating  is  nearly 
doubled  by  roughness  effects.  The  measurements  made 
on  the  nosetip  at  Mach  13  are  shown  in  Figures  13  and 
14.  As  we  observed  in  the  smooth  wall  studies,  increasing 
Mach  number  increases  the  nondimensiona!  heating  rate, 
and  the  comparison  shown  in  Figure  18  indicates  that  this 
is  also  true  for  rough  walls. 

It  is  clear  from  these  studies  that  the  scale  of  the 
body  shock-boundary  layer,  as  well  as  the  peak  interaction 
heating,  are  strongly  controlled  by  the  scale  of  the  surface 
roughness  anc  the  freestream  Mach  number.  We  measured 
increases  in  the  heating  coefficients  of  close  to  100 
percent  from  either  increasing  Mach  number  or  adding 
surface  roughness,  resulting  in  a  combined  enhancement 
effect  of  close  to  300  percent.  Because  of  the  formation 
of  separated  embedded  interaction  regions  on  these 
indented  noseshapes,  the  blunt  body  Mach  number 
independence  principal  cannot  be  expected  to  apply 
directly  to  these  flows. 

This  study  underlined  the  importance  of  careful 
selection  of  the  scale  of  the  surface  roughness,  in  addition 
to  testing  at  the  highest  Mach  number  for  which  fully 
turbulent  conditions  can  be  generated..  The  large  gradients 
generated  in  the  interaction  region  makes  it  imperative 
that  streamwise  heat  transfer  conduction  along  the  model 
surface  be  minimized  to  obtain  accurate  heat  transfer 


The  purpose  of  the  present  study  was  to  obtain 
detailed  heat  transfer  and  pressure  measurements  on  a 
model  of  the  NRV  nosetip  at  a  Mach  number  and  Reynolds 
number  which  duplicated  flight  conditions.  The  model 
was  highly  instrumented  to  examine  the  heat  transfer  and 
pressure  in  the  re-attachment  compre.tsion  regions  and 
close  to  the  axial  ridges.  We  wished  to  examine  whether 
the  flow  was  steady,  to  look  for  evidence  of  the 
streamwise  vortices,  and  to  determine  the  sensitivity  to 
changes  in  Reynolds  number. 

Model  and  Instrumentation 

A  4X  scale  model  of  the  NRV  nosetip  was 
constructed  for  this  program  so  that  we  could  obtain 
detailed  spatially  resolved  measurements  and  the  ability 
to  duplicate  the  Reynolds  numbers  obtained  in  the  flight 
tests.  The  model  was  constructed  using  a  precision  three- 
dimensional  pantagraph  milling  machine  with  a  epoxy 
repiica  of  the  recovered  nosetip  as  the  original.  The 
model  was  instrumented  with  100  heat  transfer  gages  and 
43  pressure  gages.  The  model  is  shown  mounted  in  the 
96-Inch  Shock  Tunnel  in  Figure  16. 

Results  and  Discussion 

A  schlieren  photograph  of  the  flow  over  the  NRV 
model  is  shown  in  Figure  17.  This  photograph  illustrates 
the  two  classes  of  flow  with  embedded  shocks  which  are 
typical  of  the  flow  over  indented  nose  shapes.  The  flow 
over  the  top  of  the  model  separates  as  it  expands 
downstream  of  the  spherical  cap  and  a  small  three- 
dimensional  separated  region  is  formed  in  the  cavity 
formed  by  the  grooves.  The  compression  waves  generated 
in  the  reattachment  region  coalesce  to  form  a  single 
shock— the  re-compression  shock.  The  bow  shock  and  re¬ 
compression  shock  interact  with  the  formation  of  a  single 
shock  and  a  shear  layer.  The  flow  over  the  bottom  of 
the  model  remains  attached:  however,  the  flow  over  the 
roughly  conical  section  of  the  nosetip  was  again 
characterized  by  a  re-compre.ision  process.  Here, 
however,  re-compression  takes  pl.ice  through  a  series  of 
oblique  shocks  as  a  wall  jet  is  formed.  It  is  the  structure 
and  development  of  the  boundary  layer  at  the  base  re¬ 
compression  processes  which  controls  the  large  heat 
transfer  rates  which  are  developed  in  these  regions. 
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The  distributions  of  heat  transfer  and  pressure  over 
the  NRV  nosetip  made  in  this  study  are  shown  in  Figures  IS 
through  21.  A  typical  heat  transfer  distribution  for  free 
stream  conditions  closest  to  those  encountered  in  flight 
are  shown  in  Figure  18a.  The  heat  transfer  measurements 
on  the  sperical  cap  indicate  that  the  boundary  layer 
remains  laminar  over  this  part  of  the  model,  with 
transition  occurring  in  the  free  shear  layer  downstream 
of  the  shoulder  expansion.  While  the  heat  transfer  rates 
drop  below  tne  stagnation  point  value  in  the  separated 
regions  just  downstream  of  the  shoulder,  heat  transfer 
rates  in  the  reattachment  and  re-compression  regions  over 
the  model  can  rise  to  close  to  three  times  the  stagnation 
poirt  value.  This  increase  is  influenced  by  transition,  the 
inci  eased  density  and  thinning  01  the  boundary  layer  in 
the  re-compression  region  and  entropy  swallowing. 
However,  the  measurements  clearly  indicate  that  the 
surface  recession  would  be  largest  on  the  conical  frustum 
of  this  nosetip.  The  pressure  distribution  over  the  nosetip 
for  the  conditions  discussed  /.,  shown  in  Figure  18b.  As 
we  would  anticipate  from  our  observations  on  the  schlieren 
photographs  and  the  heat  transfer  measurements,  the 
pressure  rises  along  the  roughly  conical  portion  of  the 
i.osetip  as  a  result  of  the  re-compression  process  and 
reaches  close  to  89  percent  of  the  stagnation  pressure 
toward  the  base  of  the  nosetip.  However,  there  is  a 
significant  variation  in  the  pressure  around  the  conical 
segment  of  the  nosetip  as  a  result  of  the  different  re- 
compress'on  processes  asso'- iated  ruth  the  circumferential 
variations  in  surface  geometry. 

During  the  design  of  this  experiment,  it  was 
speculated  that  the  flu  over  the  NRV  nosetip  would  be 
very  sensitive  to  the  Reynolds  umber  and  Mach  number 
of  the  freestrtam,  because  boundary  layer  transition  and 
separation  were  intrinsic  feat  .ires  of  the  flow  field. 
Therefore,  a  large  scale  replica  of  the  NRV  nosetip  was 
constructed  so  that  the  flight  Reynolds  number  could  be 
duplicated  in  the  Shock  Tunnel.  Thus,  it  was  of  interest 
to  examine  the  effects  of  Reynolds  number  on  the  flow 
structure  and  the  distribution  of  heat  transfer  and  pressure 
over  the  nosetip.  Figures  19a  and  19b  show  the 
distribution  of  heat  transfer  and  pressure  over  the  NRV 
nosetip  for  conditions  where  we  have  lowered  the  Reynolds 
number  to  1  x  10^.  This  results  in  an  mcrease  in  the 
size  of  the  separated  region  and  the  aft  movement  of  the 
re-compression  regions.  It  also  results  in  a  decrease  in 
pressure  rise  (relative  to  the  stagnation  point  pressure) 
on  the  conical  segment  of  the  nosetip  and  a  dramatic 
reduction  (a  factor  of  3  relative  to  the  stagnation  point 
heat  transfer)  in  the  heat  transfer  to  the  conical  frustum. 
These  results  underline  the  importance  of  conducting 
experimental  studies  of  the  flow  over  complex  indented 
n  se  shapes  at  Reynolds  numbers  and  Mach  numbers  as 
close  to  the  flight  :onditions  as  possible.  Further 
measurements  made  at  Mach  11  and  13  with  Reynolds 
numbers  of  5  x  l06/ft,  which  are  shown  in  Figures  20  and 
21,  demonstrate  that  the  heat  transfer  and  pressure 
distributori  are  sensitive  to  the  position  of  the  boundary 
layer  transition  and  the  flow  pattern,  which  are  again 
controlled  by  the  Mach  number  and  Reynolds  numbei  of 
the  freestream. 

The  measurements  made  in  this  study  clearly 
demonstrate  that  the  flow  field  and  the  distribution  of 
properties  on  "real"  ablated  rose  shapes  differ  significantly 
f.om  the  .deal.zed  bteorue  configurations  assumed  by  man/ 
analyst?  The  flow  over  NRV  configuration  if.  intrinsically 
three-dimensional  with  r  .gior.s  of  embedded  separated 
flows  as  well  as  embedded  shocks.  Experimental  studie. 
oi  the  flows  over  such  configurations  nw-t  be  r^nductsd 
at  Reynolds  numbers  and  Mach  numbers  which  e  close 


to  duplicating  those  encountered  in  flight  for  meaningful 
results  to  be  obtained. 

LAMINAR  FLOWS  OVER  INDENTED  NOSE  SHAPES 
Introduction 

When  indented  nose  shapes  occur  on  ablating  nosetips 
they  do  sc  because  the  flow  over  the  major  segment  of 
the  nosetip  is  transitional  and  turbulent.  However,  there 
have  b.  en  a  number  of  theoretical  studies  in  which 
predictions  have  been  made  of  laminar  flow  over  indented 
nose  shapes.  It  was  for  this  reason  that  the  measurements 
made  in  this  segment  of  the  study  were  obtained  to  provide 
a  detailed  set  of  measurements  in  laminar  flows  with 
which  to  examine  the  validity  of  such  solutions.  We 
studied  the  laminar  flow  over  two  ablated  nose  shapes  for 
which  Widhopf  and  Victoria^  obtained  detailed  numerical 
solutions  to  the  complete  flow  field  using  the  time- 
dependent  Navier-Stokes  equations.  These  nose  shapes 
termed  the  1-10-10  and  1-10-1  in  the  tri-radius 
terminology,  (or  the  Widhopf  1  and  2)  are  shown  in 
Figure  22. 

Model  and  Instrumentation 

Detailed  heat  transfer  and  pressure  measurements 
were  made  over  both  of  these  configurations  at  Mach  10 
for  a  Reynolds  number  (Rejq)  of  300.  The  instrumentation 
and  recording  equipment  were  designed  to  enable 
frequencies  up  to  1  Mz  to  be  recorded  with  little 
attenuation.  Widhopf's  calculations  for  the  1-10-10 
configuration  showed  a  high  frequency  instability  occurred 
in  the  flow  field  close  to  the  point  of  surface  inflexion 
on  this  body.  This  localized  oscillation  was  predicted  to 
have  a  major  frequency  of  approximately  600  KHz.  Thus, 
the  heat  transfer  gages  were  constructed  so  that  their 
physical  dimensions  were  smaller  than  the  wavelength  of 
the  oscillation. 

Results  and  Discussions 

The  heat  transfer  distribution  to  the  1-10-10 
configuration  is  shown  in  Figure  23  together  with  the 
predictions  from  the  Navier-Stokes  elutions  of  Widhopf. 
Although  the  measurements  are  m  good  agreement  with 
the  mean  values  from  these  calculations,  we  were  unable 
to  detect  the  flow  instabilities  observed  by  Widhopf  m 
his  numerical  *  k  ons.  Measurements  were  made  at  a 
diffe>ent  Mach  nur,  ;r  and  Reynolds  number  to  determine 
whether  an  oscihat.on  could  be  induced  at  slightly 
different  flow  conditions,  but  again  we  found  a  stable 
flow.  However,  one  of  the  most  interesting  aspects  of 
the  aerodynamics  of  this  configuration  is  the  large  heating 
rates  generated  at  the  inflexion  of  the  concave  surface 
at  the  base  of  the  embedded  shock  (see  Tigure  24).  Studies 
of  shock  wave-turbulent  boundary  layer  interaction  on 
concave  bodies  suggest  that  an  even  larger  increase  m 
heat  transfer  would  occur  on  the  1-10-10  configuration  if 
tne  boundary  layer  were  transitional  or  turbulent.  An 
ablati  body  upon  which  these  heating  levels  were 
impressed  would  tend  to  become  more  highly  indented  and 
possibly  develop  onto  the  shape  of  the  1-10-1 
configuration. 

Heat  transfer,  pressure  and  schheren  records  from 
measurements  made  on  the  1-10-1  configuration  at  Mach 
10  demonstrate  that  the  flow  is  fully  pulsating  m  the  E 
oscillation  mode.^  The  schheren  photographs,  shown  m 
Figure  25  indicate  the  oscillation  is  symmetric  with  a 
frequency  which  gives  a  Strouhai  number  (FD IV)  of  0.19. 
The  heat  transfer  and  pressure  to  the  model  (shown  ,n 
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i  gure  26)  undergoes  a  cyclic  variation  with  the  heat 
transfer  varying  between  a  value  which  is  approximately 
equal  that  to  the  stagnation  point  on  a  flat  faced  cylinder 
with  the  same  body  diameter,  to  a  level  six  times  as 
large.  This  later  condition  is  attained  during  the  part  of 
the  cycle  where  the  separated  shear  layer  impinges  onto 
the  almost  vertical  face  of  the  afterbody.  A  small  angle 
of  attack  can  induce  a  second  asymmetric  oscillational 
mode  during  which  a  circumferential  movement  of  the 
separated  region  occurs.  While  in  laminar  flows,  the  basic 
frequency  of  the  E  oscillational  mode  and  the  higher 
harmonics  appear  to  be  little  influenced  by  the  Mach 
number  or  Reynolds  number  of  the  freestream,  for 
turbulent  flows  over  configurations  which  are  close  to  the 
oscillation  boundaries  these  effects  do  become  of 
importance.  In  Part  11  of  this  study,  we  describe  an 
investigation  of  the  occurence  and  structure  of  unsteady 
flows  over  symmetrical  and  asymmetric  indented 
noseshapes. 

CONCLUSIONS 

Experimental  studies  have  been  conducted  to 
examine  the  structure  ar.',  properties  of  steady  flows  over 
indented  nosetips.  In  this  three-part  study  we  first 
investigated  the  influence  of  boundary  layer  trips,  surface 
roughness  and  the  Mach  number  and  unit  Reynolds  number 
on  the  flow  field  and  distribution  of  heat  transfer  and 
pressure  over  an  axisymmetric  indented  noseshape.  At 
Mach  1 1  it  was  found  that  the  surface  properties  and  flow 
field  was  uninfluenced  by  tripping,  but  introducing  uniform 
surface  roughness  increased  the  size  of  the  separated 
region  at  the  base  of  the  body  shock  and  doubled  the 
maximum  heating  rate  in  this  region.  Both  on  the  rough 
and  smooth  configuration  the  maximum  heat  transfer 
coefficient  was  doubled  when  the  Mach  number  was 
increased  from  11  to  12.  It  is  clear  that  to  correctly 
simulate  the  characteristics  of  the  flow  over  indented 
noseshapes,  the  surface  roughness  of  the  nosetip,  the  Mach 
number  and  Reynolds  number  of  the  freestream,  the 
incidence  of  the  nosetip  must  be  closely  duplicated. 

In  the  second  investigation,  in  which  a  4X  replica 
of  the  recovered  NRV  nosetip  was  studied  at  conditions 
which  duplicated  those  obtained  in-flight,  it  was 
demonstrated  that  the  flow  field  and  distribution  of 
pressure  and  heating  differed  significantly  from  the 
axisymmetric  analog.  The  flow  over  the  NRV  nosetip  is 
intrinsically  three-dimensional  and  contains  regions  of 
embedded  separated  flows  and  complex  viscous/inviscid 
interaction. 

The  study  of  the  laminar  flow  over  indented 
noseshapes  was  performed  specifically  for  a  validation  of 
a  solution  generated  by  Widhopf  and  Victoria  using  their 
Navier-Stokes  code.  While  we  did  not  find  the  high 
frequency  oscillation  predicted  by  the  numerical  solution, 
the  time  averaged  predictions  were  in  good  agreement 
with  the  experimental  measurements  made  in  the  current 
studies. 

In  part  two  of  this  study  we  will  examine  the 
occurrence  and  properties  of  unsteady  flows  over  indented 
noseshapes. 
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(a)  SMOOTH  MODEL 


(b)  MODEL  COATED  WITH  15-MI L-GRIT  ROUGHNESS 
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A  REVIEW  OF  AEROTHERMAL  PROBLEMS  ASSOCIATED 
WITH  HYPERSONIC  FLIGHT 


Michael  S.  Holden* 
Calspan/UB  Research  Center 
Buffalo,  NY  14225 


ABSTRACT 


A  review  is  presented  of  some  important  aerothermal 
problems  associated  w'th  hypersonic  flight  with  emphasis 
on  those  involved  with  viscous/inviscid  interaction.  The 
flow  regimes  of  importance  to  hypersonic  vehicles  of 
current  interest  are  discussed  and  the  performance  of 
experimental  facilities  currently  available  to  study  real 
gas  phenomena  associated  with  hypervelocity  flight  are 
briefly  reviewed.  Viscous  interaction  and  real  gas  effects 
in  hypervelocity  flight  at  high  altitudes  are  discussed 
with  reference  to  the  measurements  which  are  required 
to  support  the  development  of  improved  prediction 
capabilities  in  this  regime.  Theoretical  and  experimental 
research  on  2D  and  3D  shock  wave/laminar  boundary 
layer  interaction  relevant  to  the  performance  of  control 
surfaces  and  inlets  in  hypersonic  flight  at  high  altitudes 
is  reviewed.  The  characteristics  of  transition  regions  in 
hypersonic  flows  are  discussed  with  particular  reference 
to  transitional  flows  over  nosetips  and  swept  leading 
edges.  Studies  of  shock  interaction  phenomena  on 
indented  nosetips  are  discussed,  followed  by  a  review  of 
the  surface  roughness  and  blowing  effects.  Research  to 
describe  the  characteristics  of  2D  and  3D  regions  c-f 
shock  wave/turbulent  boundary  layer  interaction  in 
hypersonic  flow  is  reviewed  with  particular  rofcrcime  to 
intrinsic  unsteadiness  and  potential  compressibility 
effects  in  these  flows.  This  review  is  concluded  with  a 
brief  discussion  of  observations  suggesting  a  decrease  ir, 
the  spreading  rate  in  shear  layer  mixing  in  hypersonic 
flow. 


I.  INTRODUCTION 


With  the  current  interest  in  the  design  of  trans- 
atmospheric  vehicles,  maneuvering  re-entry  vehicles, 
orbital  transfer  vehicles  and  smaller  space  shuttles,  ‘here 
has  been  ah  Increased  interest  within  the  aerospace 
community  in  aerothermodynamic  problems  associated 
with  hypersonic  flight.  Following  Jim  extensive  programs 
conducted  in  the  sixties  and  early  seventies  most  subse¬ 
quent  research  in  hypersonic  flows  has  been  limited  to 
the  support  of  specific  and  vsiy  conservatively  designed 
systems  such  as  the  Space  Shuttle,  the  lovirn  entry 
vehicle,  and  ballistic  and  maneuverable  re-ei  ,-y  vehicles. 
With  the  advent  of  interest  in  the  mote  complex, 
sophisticated,  and  by  necessity,  more  efficient  vehicles 
such  as  terminally  guided  maneuvering  re-entry  vehicles, 
has  come  the  realization  that  our  current  predictive 
capability  is  not  able  to  support  the  design  of  such 
systems,  Some  of  the  fundament  d  problem  area'*  which 
must  be  addressed  ares  the  fundamental  understanding 
and  prediction  of  flows  in  the  low  density  transition 
regime  in  the  presence  of  real  gas  effects  and  surface 
catalysis;  the  control  >f  lamina,  viscous/inviscid  inter¬ 
action  regions  of  non-vv,ui  librium  airflows  on  inlets  and 
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on  control  surfaces  of  high  flying  hypersonic  vehicles; 
the  prediction  of  transition  on  complex  (and  simple) 
vehicle  configurations;  th*>  prediction  of  regions  of  shock 
wave/ turbulent  boundary  interaction,  chemically  reacting 
turbulent  shear  flows  and  the  structure  of  plasmas  in 
the  turbulent  boundary  layer  over  an  ablating  or 
transpiration-cooled  vehicle. 


Hypersonic  flow  regions  of  viscous/inviscid  inter¬ 
action  in  laminar  and  turbulent  flows  represent  one  of 
the  most  taxing  group  of  aerothermal  problems  for  the 
designer  and  certainly  are  the  most  difficult  to  predict 
in  detail  with  any  accuracy.  The  severity  of  the  gradients 
in  heat  transfer  and  pressure,  coupled  with  the  intrinsic 
unsteadiness  of  separated,  turbulent  regions  of  shock 
wave/boundary  layer  interaction  make  the  computation 
of  such  flows  one  of  most  challenging  in  high  speed 
flows.  During  the  past  two  decades  there  has  been 
significant  advances  resulting  from  fundamental 
theoretical  and  experimental  research  in  shock  wave¬ 
boundary  layer  interaction  in  supersonic  and  hypersonic 
flows.  Pi  section  techniques  have  1  .n  successfully 
developed  to  describe  laminar  separated  .lows,  first  using 
Integra!  techniques  to  solve  the  boundary  layer  equations 
and  then  through  use  of  the  full  Navier-Stokes  equations, 
solutions  to  these  eauations  have  yet  to  be  obtained  for 
laminnv  separated  flew  with  non-equilibrium  chemistry. 
However,  two-dimensional  turbulent  separated  flows 
induced  in  regions  of  shock  wave  boundary  layer 
interaction  have  not  yet  yielded  to  a  detailed  numerical 
analysis  using  the  full  Navier-Stokes  equations.  In  many 
cases  the  severe  heating  loads  and  gradients  generated 
ir  regions  of  shock/boundary  layer  interaction  mandate 
the  use  of  thermal  protection  systems  employing  ablative 
or  transpiration  cooling.  In  each  case  surface  roughness 
is  intrinsically  a  factor  resulting  from  either  surface 
structure  of  the  ablated  surface  or  the  geometric  design 
of  the  transpiration-cooled  surface.  The  effects  el 
surface  roughness  arid  biowing  will  have  a  large  e'fect 
on  the  size  and  characteristics  of  regions  of  shock  wave- 
turbulent  boundary  layer  interaction  in  hypersonic  flow. 


In  the  following  paper  wo  review  some  important 
aerothermal  problems  associated  with  flight  at  hypersonic 
speed,  with  particular  emphasis  on  t  'gions  of 
viscous/inviscid  interaction  in  hypersonic  flows.  We  first 
discuss  the  flow  regimes  which  are  of  importance  lor 
the  Hypersonic  vehicles  of  current  interest  and  facilities 
which  are  currently  available  to  study  hypervelocity 
airflows.  Viscous  interaction  and  rea)  gas  effects  in  the 
high  altitude  transition  regime  are  briefly  reviewed. 
Studies  of  laminar  2D  and  3D  viscous  interaction  regions 
are  then  discussed  with  emphasis  on  comparisons  between 
computational  methods  and  experimental  measurements. 
Research  on  boundary  layer  transition  and  transitional 
Hows  over  nosetips  and  leading  edges  at  hypersonic  speeds 
are  briefly  reviewed.  The  regions  of  shork/boundary 
layer  interaction  which  control  the  nose  shaping  ol 
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-ablating  nosetips  are  described  with  particular  reference 
to  the  importance  of  Mach  number  and  Reynolds  number 
simulation  in  experimental  studies  of  these  flows.  The 
characteristics  of  2D  and  '3D  shock  wave/turbulent 
interactions -in  hypersonic  flows  are  then  described  again 
with  particular  reference  to  solutions  obtained  to  the 
Navier-Stokes  equations.  We  conclude  with  discussion  of 
compressibility  and  unsteady  effects  in  hypersonic  flows. 

H.  FLIGHT  REGIMES  AND  PHENOMENA  OF  INTEREST 
AND  GROUND  TEST  SIMULATION  CAPABILITIES 

While  the  accurate  prediction  of  three-dimensional 
turbulent  regions  of  shock  wave/boundary  layer 
•interaction-with  embedded,  separated  regions  over  rough 
blowing  surfaces  represents  one  of  the  most  difficult 
problems  from  the  viewpoint  of  analysis,  it  is  the 
occurrence  of  extensive  laminar  and  transitional  regions 
of  viscous/irviscid  interaction  which  may  present  the 
designer  of  high  flying  hypersonic  vehicles  with  the 
greatest  design  problems.  A  plot  of  the  Mach 
number/Reynolds  number  and  altitude/  velocity 
capabilities  of  a  number  of  high  performance  hypersonic 
facilities  are  shown  together  with  the  flight  envelopes 
of  a  number  of  proposed  vehicles  in  Figures  1  and  2. 
Because  of  the  increased  stability  of  the  laminar  boundary 
layer  at  high  Mach  numbers,  the  trans-atmospheric  and 
the  orbital  transfer  vehicles  will  be  immersed  in  laminar 
boundary  layers  for  a  significant  part  of  their  flight. 
Laminar  viscous/inviscld  interaction  regions  for  large 
values  of  the  viscous  interaction  parameter  M  Vc*/>5o  will 
cause  a  significant  reduction  in  the  control  forces  that 
can  be  developed  by  elevons  or  reaction  jet  controls. 
As  illustrated  in  Figure  3,  extensive  regions  of  laininar 
attached  or  separated  flow  can,  In  essence,  "fair-in"  the 
compression  surface,  reducing  their  effectiveness. 
Likewise,  the  boundary  displacement  effects  may 
significantly  modify  the  effective  shape  of  the 
compression  tamp  on  a  hypersonic  intake.  Because  of 
the  intrinsic  structure  of  hypersonic  laminar  boundary 
layers  over  cooled  surfaces  it  may  be  extremely  difficult 
to  develop  an  affective  boundary  layer  control. 

For  at  least  two  concepts  currently  under  study,  the 
orbital  transfer  vehicle  and  the  trans-atmospheric 
vehicle,  designed  to  operate  at  very  high  altitudes,  it 
has  yet  to  be  established  that  the  Navier-Stokes  codes 
wlil  successfully  handle  flows  on  the  boundaries  of 
continuum  flow,  because  of  the  high  re-entry  velocities 
associated  With  the.yr  missions  it  is  clear  that  the  airflows 
over  the  vehicles  v/ili  be  in  chemical  non-equilibrium, 
and  again  incorporat-,ng  non-equilibrium  chemistry  into 
the  Navier-Stokes  has  yet  to  be  accomplished  successfully 
for  the  hypersonic  flow  regime.  Surface  catalysis  must 
be  considered,  and  here  errors  of  up  to  300  percent  can 
result  from  uncertainties  in  the  current  prediction 
techniques  for  heat  transfer. 

The  prediction  of  boundary  layer  transition  remains 
one  of  the  most  important  and  difficult  problems  in 
hypersonic  flow.  The  theoretical  studies  have  failed  to 
provide  either  definitive  prediction  methods  or  the 
framework  within  which  to  correlate  the  results  from 
the  experimental  studies.  Of  the  many  ways  which  have 
been  suggested  to  correlate  transition  measurements  from 
flight  and  ground  test  facilities,  describing  the  beginning 
of  transition  in  terms  of  a  local  Reynolds  number  based 
on  either  the  local  momentum  thickness  or  streamwise 
distance  and  a  local  Mach  number  remains  the  most 
popular.  A  typical  example  of  a  transition  correlation 
of  flight  and  ground  test  measurements  is  shown  in 
Figure  <f.  The  ground  test  measurements,  specifically 
those  made  on  swept  cylinders,  are  in  surprisingly  good 


agreement  with  those  made  in  flight,  plotted  in  this 
format.  This  good  agreement  may  result  in  part  from 
the  relative  insensitivity  of  stagnation  line  transition  to 
facility-related  disturbances.  in  many  cases  surface 
roughness  can  play  an  important  role  in  controlling 
transition  on  flight  vehicles,  and  correlations  employing 
the  Reynolds  number  based  on  momentum  thickness  and 
the  roughness  height  in  a  form  similar  to  that  shown  in 
Figure  5  can  be  used  to  extrapolate  ground  test  data  to 
flight.  On  more  complex  aerodynamic  configurations 
however,  where  adverse  pressure  gradients  are  generated 
in  viscous  interaction  regions  over  controls  or  inlets, 
transition  can  be  expected  to  occur  within  these  regions 
at  Reynolds  numbers  of  an  order  of  magnitude  less  than 
those  presented  above. 

The  largest  and  most  difficult  to  predict  aerothermal 
loads  which  are  experienced  by  an  advanced  hypersonic 
vehicle  are  encountered  under  low  altitude,  turbulent 
flow  conditions.  The  severity  of  the  heat  transfer  rates 
and  gradients  encountered  in  regions  of  shock/shock  and 
shock/boundary  layer  interaction  make  it  essential  for 
the  designer  to  have  prediction  techniques  available  which 
will  describe  these  regions.  As  discussed  later  in  the 
context  of  flows  over  indented  nosetips,  to  obtain  an 
accurate  simulation  of  shock/boundary  layer  interaction 
phenomena  in  a  ground  test  facility,  it  is  important  that 
both  the  Mach  number  and  Reynolds  number  be  simulated. 
However,  if  the  boundary  layers  and  shear  layers  in  the 
experiments  are  fully  turbulent  (without  resorting  to 
tripping)  the  characteristics  of  the  interaction  regions 
do  not  appear  to  exhibit  a  strong  Reynolds  number 
dependence.  There  is,  however,  a  strong  effect  of  Mach 
number  on  the  structure  of  viscous  interaction  regions 
and  Mach  number  independence  is  not  valid  for  blunt 
bodies  with  embedded  shock  systems.  Also  for  the 
correct  simulation  of  hypervelocity  flows,  the  effects  of 
real  gas  chemistry  must  be  included. 

For  velocities  above  10,000  ft/scc  where  oxygen  is 
fully  dissociated  and  nitrogen  dissociation  and  ionization 
start  to  occur,  surface  related  reactions  will  begin  to 
influence  heating,  and  real  gas  effects  will  start  to 
influence  the  structure  and  characteristics  of  regions  of 
shock  wave/boundary  layer  interaction.  The  first  order 
effects  of  air  chemistry  are  to  change  the  specific  heat 
of  the  gas  and  thus  the  pressure/Mnch  relationships  in 
the  viscous  and  inviscid  flow  equations.  It  can  be  seen 
from  Figures  1  and  2  that  with  the  exception  of  Shock 
Tunnels  most  of  the  high  Reynolds  Number  hypersonic 
facilities  employ  nitrogen  ns  a  test  gas  and  cannot 
generate  velocities  above  6000  ft/sec. 

There  are  very  few  facilities  which  are  currently 
operational  in  which  detailed  studies  can  be  performed 
to  study  real  gas  effects  in  hypersonic  flows.  With  the 
exception  of  ballistic  ranges,  the  shock  tunnels  at  Calspan 
and  piston  driven  shock  tunnel  (or  Stalker  Tube)  in 
Australia  are  the  only  operational  high  Mach  number 
facilities  which  can  be  used  to  generate  clean  airflows 
at  the  temperatures  of  above  7000°K  (see  Figure  6) 
necessary  to  excite  both  nitrogen  and  oxygen.  The 
tremendous  power  that  is  being  expended  during  the 
running  time  of  these  facilities  (typically  a  million 
horsepower)  underscores  why  the  construction  of  a 
facility  to  operate  at  similar  conditions  with  test  times 
of  over  I  second  is  a  formidable,  if  not  unrealistic,  task. 
However,  there  is  a  clear  need  for  facilities  capable  of 
generating  high  enthalpy  airflows  with  low  ambient 
freestream  dissociation  levels  and  sufficient  test  time 
for  combustion  and  transpiration  cooling  processes  to 
become  established.  Our  experience  suggests  that  most 
simple  boundary  layer  flows  are  established  in  a  fraction 
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of  a  millisecond  while  more  complex  flow  fields,  involving 
blowing  and. flow  separation,  can  take  up  to  3  milliseconds 
to  establish.  A  high  enthalpy  facility  with  run  times  of 
up  to  10  milliseconds  would  handle  most  flows  of  interest 
with  ease. 

Employing  an  isentropic  compression  process  to 
compress  the  test  gas  (rathe;  than  a  reflected  shock) 
results  in  a  reservoir  which  has  a  lower  level  of 
dissociation,  ar..J  for  large  drive  pressure  ratios  can  lead 
to 'higher  stagnation  pressures  and  temperatures.  The 
isentropic  compression  tube  was  one  of  the  concepts 
under  study  in  die  late  sixties  along  with  the  MHD 
augmented  shock  tunnel*,  and  the  free  piston  driven 
shock  tunnel  and  expansion  tube**.  It  was  initially  hoped 
that  stagnation  temperatures  of  up  to  10,0G0oK  could  be 
generated  in  the  Gun  and  Longshot  tunnels,  where  the 
gas  is  compressed  by  a  free  piston  motion  rather  than 
a  program  piston  trajectory  for  the  isentropic  tube  (sec 
Figure  7).  The  aerotherinal  loads  on  the  piston  have 
been  such  that  neither  facility  has  achieved  its  potential. 
Currently  stagnation  temperatures  of  less  than  2000°l< 
are  obtained  in  Gun  Tunnels  and  the  Longshot  facilities 
(see  Figures  8  and  9);  their  value  lies  in  their  ability  to 
generate  relatively  high  Mach  number/high  Reynolds 
number  with  little  cost.  While  interest  in  the 
development  of  high  enthalpy  hypersonic  facilities 
dwindled  in  the  U.S.  during  the  late  sixties  and  early 
seventies,  Stalker  continued  work  on  a  facility  to 
generate  real  gas  flow  at  temperatures  up  to  9000°K. 
As  shown  in  Figure  10,  the  Stalker  tunnel  is  in  essence 
a  shock  tunnel  with  a  compression  tube  used  to  supply 
the  driver  gas.  This  feature  eliminates  the  problems 
with  operating  a  high  temperature,  high  pressure  driver 
section  at  the  expense  of  a  more  complex  device  and 
one  with  shorter  run  times  than  the  shock  tunnel,  which 
at  the  high  enthalpy  conditions  are  less  than  1 
millisecond.  ludicious  use  of  current  facilities  can 
provide  insight  to  many  of  these  nonequilibirum  problems, 
but  it  is  clear  that  further  inventive  development  of 
hypervelocity  facilities  are  required  if  vehicles  like  the 
hypersonic  transport  are  to  be  actively  developed. 
However,  in  the  short  term,  problems  with  real  gas 
effects,  similar  to  those  discussed  in  the  following 
section,  can  at  best  be  addressed  with  the  limited  results 
achievable  in  flight  tests. 

IN-  VISCOUS  INTERACTION  AND  REAL  GAS  EFFECTS 
IN  HYPERVELOCITY  FLIGHT 

The  recent  interest  in  hypersonic  flight  at  the  edges 
of  the  atmosphere  has  rekindled  interest  in  hypervelocity 
low  density  flows  where  viscous/inviscid  interaction  and 
chemical  non-equilibrium  effects  are  important.  The 
success  of  a  hypersonic  air  breathing  propulsion  system 
designed  to  operate  at  very  high  altitudes,  and  the 
aerothermal  performance  of  the  AOTV  will  depend  on 
how  well  viscous  interaction,  real  gas  and  surface  blowing 
effects  can  be  predicted  and  controlled.  Within  the  past 
decade  great  strides  have  been  made  in  the  prediction 
of  free  molecular  flows  and  work  in  this  area  has  recently 
been  reviewed  by  Bird^.  For  continuum  flows  the 
tremendous  increase  in  computational  power  has  allowed 
boundary  layer  displacement  and  leading  edge  bluntness 
effects  to  be  computed  with  the  full  Navier-Stokes 
equations.  Fot  temperatures  up  to  6,000  °K  the  reaction 
rates  of  air  are  known  with  sufficient  accuracy**  to  permit 
solutions  to  the  full  Navier-Stokes  equations  for  laminar 
flows  with  non-equilibrium  air  chemistry  to  be  obtained; 
however,  there  remains  some  significant  numerical 
problems  associated  with  the  finite  rate  chemistry.  Shock 
layer  and  boundary  layers  solutions  with  finite  rate 
chemistry  are  now  being  routinely  obtained^.  However, 
until  there  is  significant  improvement  in  the  under¬ 
standing  of  surface  catalysis,  such  solutions  will  be  valid 


only  for  either  fully  catalytic  or  non-catalytic  surfaces, 
and  the  differences  in  heat  transfer  between  these  cases 
can  amount  to  factors  of  up  to  2.  One  case  in  point  is 
current  anomalies  associated  with  measured  and  com¬ 
puted  heat  transfer  rates  from  STS  flights  1,  2  and  3. 
While  for  altitudes  above  150,000  feet  the  codes  predict 
nonequitibrium  effects  associated  with  a  lag  in  the 
nitrogen  recombination  rate  (oxygen  remains  fully 
dissociated),  the  best  comparisons  have  been  obtained 
with  calculations  based  on  equilibrium  air  properties*’ 
(see  Figure  11).  Since  in  hypersonic  flows  the  windward 
pressures  reflect  the  momentum  of  the  incident  molecules 
(Newtonian  theory)  and  not  their  chemical  state,  they 
are  only  weakly  influenced  by  nonequilibrium  effects. 
When  such  effects  become  large  they  are  generally 
associated  with  regions  of  strong  expansion  or  smooth 
compression  process  where  their  effects  on  2'  are 
important.  Under  such  circumstances  the  pressure  is 
also  strongly  influenced  by  viscous/invisicid  interaction 
effects  as  illustrated  by  wedge/flat  plate  shock  tunnel 
studies  of  Vidal  and  Stoddard7  where  the  pressures  on 
the  flat  plate  downstream  of  a  strong  expansion  were 
influenced  by  viscous  effects  to  a  greater  extent  than 
by  chemical  non-equilibrium.  ‘  For  the  flows  over  the 
Space  Shuttle,  in  the  absence  of  viscous  interaction,  real 
gas  effects  would  result  in  higher  pressures  on  the  lecside 
and  slightly  lower  pressures  on  the  windward  side. 
However,  the  (low  over  the  control  surfaces  are 
significantly  influenced  by  viscous  interaction  effects  as 
are  those  on  the  leeside,  and  to  reconcile  the  differences 
between  force  measurements  made  in  wind  tunnels  and 
flight  solely  on  the  basis  of  an  equilibrium  inviscid  code 
to  correct  for  real  gas  chemistry  (as  suggested  by  Maus 
and  Griffith**)  might  be  an  interesting  but  slightly 
misleading  first  step.  It  is  of  interest  to  note  that  the 
greater-than-expccted  positive  pitching  moment 
experienced  in  flight  could  also  be  explained  directly  in 
terms  of  a  loss  in  flap  effectiveness  as  a  result  of 
increased  size  of  the  interaction  region  at  the  lower 
Reynolds  numbers  encountered  in  flight.  The  force 
measurements  made  on  the  shuttle  model  at  low  Reynolds 
numbers  in  the  Calspan  Shock  Tunnel,  but  apparently  not 
included  in  the  design  data  book  (shown  in  Figure  ’  2) 
add  support  to  this  contention. 

While  in  principle  it  is  possible  to  describe  laminar 
non-equilibrium  interacting  flows  exactly  within  the 
framework  of  the  Navier/Stokcs  equations  once  catalytic 
wall  effects  are  understood,  at  the  higher  altitudes  whore 
wall  slip  effects  become  important  the  governing 
equations  are  not  as  well  defined.  However  the  general 
success  obtained  with  extending  simple  viscous  layer 
methods  into  the  transition  regime  using  simple  slip 
models9  (see  Figure  13)  suggests  that  extending  the 
Navier/Stokes  codes  to  describe  transitional  flows  using 
a  similar  approach  could  be  productive.  It  is  in  this 
flow  regime  where  nonequi librium,  viscous  interaction 
and  wall  slip  effects  are  important  that  there  is  a  signifi¬ 
cant  need  for  detailed  theoretical  and  experimental 
research.  While  there  are  surface  measurements  of 
pressure  and  heat  transfer  available'**  there  is  a  clear 
need  for  measurements  in  high  velocity  airflows  of  the 
detailed  flow  structure,  particularly  near  the  wall,  to  aid 
in  modeling  the  slip  flow.  These  measurements  must  be 
made  in  high  enthalpy  airflows,  where  real  gas  effects 
are  important,  to  capture  some  of  the  essential  physics 
needed  to  check  the  modeling.  Complementary  studies 
employing  NASA's  proposed  Entry  Research  Vehicle  and 
high  enthalpy  Shock  Tunnels  and  Plasma  .lets  could 
provide  the  required  information  on  catalytic  wall, 
viscous  interaction  and  slip  effects  in  non-equilibrium 
airflow  that  is  required  to  obtain  an  accurate  predictive 
capability  for  hypervelocity  flight  at  high  altitudes. 


IV.  LAMINAR  V1SCOUS/INV1SC1D. INTERACTION 
IN  HYPERSONIC  FLOW 

In  high  Mach  number,  low  Reynolds  number  flows, 
the  interaction  between  the  growth  of  the  laminar 
boundary  layer  and  the  outer  inviscid  flow  over  the 
nosetip  or  leading  edge,  and  subsequent  interaction  over 
control  surfaces  or  inlets  can  play  an  important  role  in 
determining  the  aero-thermodynamic  performance  of  a 
high  flying  hypersonic  vehicle.  Also,  because  most 
practical  vehicle  designs  employ  blunt  leading  surfaces, 
the  entropy  layer  can  have  a  strong  influence  on  the 
structure  of  the  boundary  layer  upstream  of  the 
interaction,  and  on  the  pressures  and  heating  rates  on 
the  compression  surfaces.  As  illustrated  in  the  schlieren 
photographs  shown  in  Figures  14  and  15  for  flow  over 
flat  plate/wedge  and  curved  compression  surface,  at  high 
Mach  numbers  it  is  difficult  to  distinguish  between  the 
separation  shock  and  the  edge  of  the  viscous  layer,  while 
in  supersonic  flow  these  features  are  distinct.  An 
important  feature  of  the  structure  of  a  hypersonic 
laminar  boundary  layer  over  a  cooled  wall  is  that  most 
of  the  mass  and  momentum  is  contained  at  the  outer 
edge  of  the  boundary  layer,  thus  making  it  difficult  to 
employ  boundary  layer  control. 

The  importance  of  viscous  interaction  and  leading 
edge  bluntness'  on  the  flow  over  a  two-dimensional 
compression  corner  are  illustrated  by  the  measurements 
of  heat  transfer,  skin  friction  and  pressure  shown  in 
Figure  15.  Here  on  the  configuration  with  the  sharp 
leading  edge,  the  viscous  interaction  extends  20  boundary 
layer  thicknesses  downstream  of  the  corner,  at  which 
point  the  pressure  has  risen  two  orders  of  magnitude  to 
the  inviscid  wedge  value.  In  contrast,  the  flap  pressure 
(and  heat  transfer)  on  the  configuration  with  the  blunt 
leading  edge  are  dominated  by  entropy  swallowing  effects 
and  rise  to  no  more  than  five  times  the  values  just 
upstream  of  the  corner  interaction.  In  this  study  it  was 
found  that  while  bluntness  was  found  to  decrease  the 
size  of  the  reverse  flow  embedded  within  the  interaction 
region,  flow  separation  occurred  at  approximately  the 
same  wedge  angle.  In  a  subsequent  study  it  was  also 
shown  that,  in  contrast  with  turbulent  flows,  the  angle 
through  which  the  flow  could  be  turned  without  inducing 
separation  could  not  be  changed  by  radi  using  the  corner 
until  the  corner  radius  approached  at  least  15  initial 
boundary  thicknesses.  While,  as  discussed  below, 
Holden1  1  was  able  to  obtain  predictions  in  reasonable 
agreement  with  measurements  made  in  these  studies 
employing  boundary  layer  equations  modified  to  incor¬ 
porate  the  normal  pressure  gradient,  it  is  clear  that  such 
predictions  are  now  best  obtained  with  numerical 
solutions  to  the  full  Navier-Stokes  equations.  However, 
without  the  massive  computational  facilities  currently 
available,  the  devisors  of  earlier  prediction  methods 
employed  the  boundary  layer  equations,  coupled  through 
equation  describing  a  displacement  surface,  to  the  "outer 
inviscid"  flow. 

Laminar  Two-Dimensional  Viscous  Interaction 

The  increased  stability  of  the  laminar  boundary  at 
high  Mach  numbers,  and  the  ease  with  which  laminar 
boundary  layers  separate,  coupled  with  interest  in  high 
altitude  hypersonic  flight,  has  made  laminar  flow 
separation  in  two-  and  three-dimensional  interacting 
flows  of  considerable  practical  importance.  Experimental 
studies  and  flight  tests  show  that  the  very  large  heat 
transfer  rates  and  gradients  generated  in  the  reattach¬ 
ment  regions  of  laminar  separated  flows  in  high  speed 
flight  are  of  considerable  importance  to  the  heat  shield 
designer.  The  simple  viscous  interaction  flow  models 


that  have  been  developed  to  describe  laminar  two- 
dimensional  separated  regions  have  their  origins  in  the 
earlier  studies  of  Howarth12,  Lighthill13  and  Oswatitsch 
and  Wieghardt 1,(,  who  proposed  models  for  the  mechanism 
of  upstream  influence  and  boundary  layer  separation.  It 
was  then  shown  that  upstream  influence  and  flow 
separation  could  be  described  with  good  accuracy  by  a 
model  in  which  the  viscous  layer  grew  by  mutual 
interaction  with  the  outer  inviscid  flow.  Cilick  • ,  using 
a  semi-empirical  modification  of  the  Crocco-Lces1” 
method  based  on  this  "free  interaction"  model,  was  able 
to  describe  successfully  the  properties  of  a  complete 
shock-induced  separated  region.  Honda17  and  later  Lees 
and  Reeves111  added  a  third  equation,  the  moment  of 
momentum  equation  and  obtained  generally  good  agree¬ 
ment  with  measurements  in  supersonic  separated  flows 
over  adiabatic  walls.  To  describe  separation  in  laminar 
viscous  interaction  regions  under  highly  cooled  wali 
conditions  and  to  predict  reattachment  heating,  Holden19 
added  the  integral  form  of  the  energy  equation  to  the 
equations  for  mass,  momentum,  and  moment  of  momen¬ 
tum  employed  by  Honda,  using  Cohen  and  Reshotko's20 
similar  solutions  as  proposed  by  Lees  and  Reeves.  Good 
agreement  was  found  between  theory  and  experiment  at 
Mach  numbers  in  the  low  hypersonic  range;  however,  for 
highly  cooled  flow  above  Mach  1 1  the  laminar  boundary 
layer  exhibited  a  "supercritical  response".  At  high  Mach 
numbers,  Holden11  showed  that  the  normal  pressure 
gradients  must  be  included  in  the  description  of 
hypersonic  interaction  regions  to  describe  such  flows 
correctly.  A  comparison  between  Holden's  theory  and 
the  measurements  of  skin  friction,  pressure  and  heat 
transfer  made  on  the  flat-plate,  18°-wedge  is  shownin 
Figure  16.  Incorporating  Sp/dy  into  the  formulation 
enabled  a  solution  to  be  obtained  where  previously 
artifices,  such  as  the  supcr-subcritical  jump,  would  have 
to  be  employed  to  overcome  the  problems  inherent  in 
conventional  boundary  layer  theory  for  cooled  wall 
hypersonic  flows. 

While  Carter21  obtained  one  of  the  first  solutions 
to  the  Navier-Stokes  equations  for  separated  flows,  the 
technique  developed  by  Hung  and  MacCormack22 
represents  the  first  relatively  efficient  method  devised 
to  predict  the  characteristics  of  laminar  separated  flows. 
Using  tjte  Navier-Stokes  code  developed  by  Hung  and 
MacCormack,  laminar  solutions  were  obtained  for  com¬ 
parison  with  the  experimental  measurements  at  Mach  16. 
Comparisons  between  the  Navier-Stokes  solutions  and 
measurements  in  attached,  incipient  separated  flows  arc 
shown  in  Figures  17a,  b,  and  c.  For  attached  flow  over 
the  flat-plate  wedge  configuration  (Figure  17ai,  the 
theoretical  predictions  are  in  excellent  agreement  with 
the  experimental  data.  In  particular,  the  form  of  the 
heat  transfer  and  skin  friction  in  the  region  of  minimum 
heat  transfer  is  well  reflected  by  the  theoretical  predic¬ 
tions.  Both  theory  and  experiment  display  a  minimum 
skin  friction  downstream  of  the  corner  or  the  wedge. 
As  the  strength  of  the  interaction  is  increased  to  obtain 
a  small  separated  region,  the  heat  transfer  dips  sharply 
ahead  of  the  corner,  displaying  a  cusp-like  form  in  the 
presence  of  a  separated  region  (See  Figure  17b).  Again, 
the  skin  friction  distribution  is  well  predicted  by  theory. 
For  the  well-separated  flow  (Figure  17c),  the  Navier- 
Stokes  solution  significantly  underpredicts  the  size  of  the 
separated  region.  While  the  form  of  the  distributions 
and  the  absolute  levels  in  the  plateau  and  reattaclunent 
region  are  well  |x-edicted,  the  separated  region  is 
approximately  0.6  of  the  length  found  in  the  experiment. 
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Laminar  Three-Dimensional  Interactions 

Most  studies  of  three-dimensional  interaction  have 
been  conducted  to  examine  flow  in  the  corner  of  two 
intersecting  wedges  with  sharp  leading  edges  aligned  at 
an  angle  of  90  degrees  with  each  other.  Stainback22’2^ 
made  detailed  heat  transfer  and  pressure  measurements 
in  laminar  flow  over  such  a  configuration  at  Mach  5, 
and  later  at  Mach  8,  for  a  range  of  Reynolds  number  to 
identify  transition.  These  experiments  led  Stainback  to 
distinguish  between  two  flow  regions:  the  near  corner 
with  which  most  previous  theoretical  work  had  been 
concerned  and  the  far  corner  which  had  been  studied 
experimentally  (see  Figure  19).  The  near  comer  was  a 
region  of  mutual  interaction  of  two  boundary  layers  and 
resulted  in  lower  heating.  The  far  corner  was  the  region 
of  mutual  interaction  of  the  two  inviscid  flow  fields 
caused  by  the  intersection  of  the  two  flat  plates  forming 
the  corner.  In  subsequent  studies,  Stainback  and 
Weinstein2?  further  observed  that  interaction  between 
the  mutual  boundary  layers  in  the  corner  results  in  a 
decrease  in  heat  transfer  in  the  very  near  corner.  Away 
from  the  corner,  the  vortex  system  and  reattachinent  of 
the  boundary  layer  downstream  of  the  shock  induced 
separation  (see  Figure  19)  results  in  an  increase  of  heat 
transfer  outboard  of  the  near  corner  region.  The  basic 
structure  of  laminar  interacting  flows  was  defined  in 
studies  by  Charwat  and  Redekeopp30  at  Mach  2  to  4. 
Studies  in  hypersonic  flow  by  Watson  and  Weinstein3! 
demonstrated  similar  features.  Figure  20  shows  the  flow 
structure  established  from  pitot  pressure  and  flow 
visualization  measurements.  The  shock  waves  generated 
by  each  wedge  are  joined  by  a  third  shock  wave,  bordering 
Zone  I  and  the  freestream  flow.  Slip  surfaces  pass  from 
the  ends  of  this  third  shock  towards  the  corner.  Shock 
legs  extend  to  the  surface  —  a  curved  shock  between 
Zones  11  and  HI  and  a  spread  of  the  corner  disturbance 
outside  of  the  inner  shocks.  Kutlcr  et  al.,32  have  shown 
that  shock  capturing  techniques  describe  the  inviscid  (low 
with  good  accuracy;  however,  viscous  characteristics, 
particularly  with  laminar  flows,  where  viscous  interaction 
is  important,  are  predicted  with  significantly  less 
accuracy. 

V.  TRANSITION  TO  TURBULENCE  IN  HYPERSONIC 
FLOWS 

While  boundary  layer  transition  is  one  of  the  most 
'important  parameters  in  the  design  of  hypersonic 
vehicles,  there  remains  a  considerable  gap  between 
"engineering"  efforts  to  correlate  the  occurrence  of 
transition  and  fundamental  theoretical  studies.  The 
engineering  studies  have  concentrated  on  correlating 
experimental  measurements  of  transition  obtained  in 
flight  tests,  ballistic  ranges,  and  wind  tunnels  against 
almost  every  conceivable  parameter,  while  the 
fundamental  studies  have  been  aimed  principally  at 
exploring  the  modes  of  instability  of  the  laminar 
boundary.  The  basic  problem  is  that  boundary  layer 
transition  is  controlled  by  the  detailed  aerodynamic 
environment  as  well  as  the  Reynolds  number  as 
demonstrated  by  Osborn  Reynolds33  in  his  classic  studies. 

Both  theoretical  and  experimental  studies  of 
transition  in  incompressible  flows  have  suggested  that, 
in  the  absence  of  large  freestream  disturbances  or 
disturbances  from  within  the  boundary  layer,  the  two- 
dimensional  Tollmien-Schlichting3^  model  is  a  good 
representation  of  the  initial  breakdown  of  laminar  flow. 
However,  when  velocity  fluctuations  of  over  10  percent 
are  Introduced  into  the.  freestream,  "by-pass"  inodes 
involving  the  generation  of  three-dimensional,  stream- 


wise,  vortex-like  instability  may  be  the  principal 
mechanisms  involved  in  the  transition  process.  Mack33 
among  others,  has  speculated  that  a  helical  instability 
mode  should  be  the  dominant  feature  of  the  flow 
upstream  of  the  non-linear  breakdown  into  turbulence  of 
a  laminar  boundary  layer  in  supersonic  flow.  However, 
the  lack  of  direct  evidence  to  support  this  prediction 
has  slowed  theoretical  developments  of  this  nature.  An 
tinderstanding  of  how  disturbances  in  the  freestream 
influence  transition  requires  knowledge  of  the  coupling 
mechanisms  between  the  laminar  boundary  and  the  outer, 
inviscid' flow.  While,  in  most  instability  problems,  it  is 
necessary  only  to  calculate  the  conditions  for  the 
existence  and  growth  of  a  disturbance,  transition 
prediction  using  linear  stabilty  theory  rests  almost 
completely  on  the  ability  to  trace  the  origin  of  the 
instability3*’. 

Experimental  studies  of  transition  are  of  little 
general  value  unless  it  can  be  shown  that  the  instability 
modes  in  the  boundary  layer  remain  unexcited  by 
disturbances  in  the  freestream  or  respond  to  stronger 
disturbances  generated  within  the  boundary  layer,  for 
example,  by  surface  roughness.  As  discussed  earlier,  the 
author  remains  openly  skeptical  of  the  general  usefulness 
of  any  transition  data  gathered  in  any  wind  tunnel. 
However,  the  agreement  between  the  measurements  of 
transition  Reynolds  number  made  in  large,  high  Mach 
number  facilities,  like  the  Calspan  Shock  Tunnel,  with 
free  flight  measurements  suggest  a  low  level  of  tunnel 
noise.  Studies  at  AEDC32  and  NASA-Langley3*  have 
demonstrated  that,  in  conventional  supersonic  wind 
tunnels,  the  Reynolds  number  at  which  transition  occurs 
is  strongly  influenced  by  the  fluctuating  pressure  level 
in  the  freestream.  Pate  and  Schueler3'  have  demon¬ 
strated  that  the  level  of  pressure  fluctuation  can  be 
related  to  the  geometric  features  of  the  tunnel  used  and 
to  the  characteristics  of  the  boundary  layer  on  the  tunnel 
walls,  and  a  decrease  in  tunnel  size  for  the  same 
freestream  conditions  should  result  in  a  decrease  in  the 
transition  Reynolds  number.  Comparing  measurements 
made  in  the  Calspan  Shock  Tunnel,  shown  in  Figure  21, 
we  observed  Re&°'z  variation  and  not  a  jump  in  Reg, 
with  tunnel  size. 

Correlations  of  the  transition  measurements  made 
on  sharp  cones  and  flat  plates  in  studies  at  Calspan, 
flight  measurements  reported  by  TRW,  and  measurements 
in  the  ballistic  ranges  at  AEDC  and  NSWC  are  shown  in 
Figure  22.  A  further  comparison  between  tunnel 
measurements  and  those  made  in  the  more  recent  studies 
of  Reda3?,  plotted  in  terms  of  the  unit  Reynolds  number, 
are  shown  in  Figure  23.  We  find  relatively  good 
agreement  between  the  two  sets  of  measurements,  and, 
in  common  with  the  studies  of  Potter?®,  Shcctz?*  and 
Reda,  we  observe  a  unit  Reynolds  number  effect.  The 
source  and  significance  of  the  unit  Reynolds  number 
effect  have  been  the  subject  of  extensive  debate.  The 
analyses  of  Morkovin?2  and  Reshotko?3  have  suggested 
that  the  unit  Reynolds  number  effects  may  be  traced  to 
a  sensitivity  to  the  non-dimensional  frequency  (Ue2/^e) 
or  to  the  wavelength  of  the  disturbance  (Ue/z^);  however, 
in  reality,  the  disturbance-inducing  transition  may  stein 
from  superposition  of  a  number  of  different  mechanisms. 

A  composite  of  the  schlicren-photograph  and  heat 
transfer  records  obtained  at  Mach  13,  shown  in  Figure 
24,  presents  a  "good  picture"  of  the  general  structure  of 
the  transition  region.  The  "spikes"  in  the  laminar  heat 
transfer  trace,  which  mark  the  beginning  of  the  transition 
process,  can  be  described  as  turbulent  bursts,  which  are 
fairly  well  structured,  and  which  travel  at  very  close  to 
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the  velocity  of  the  inviscid  flow.  As  these  bursts  are 
converted  downstream,  they  begin  to  break  up  into  large- 
scale  instabilities.  The  convection  of  these  large-scale 
eddies  past  a  thin-film  gage  gives  rise  to  the  intermittent 
characteristics  shown  in  Figure  24.  As  the  spectrum  of 
turbulent  scale  sizes  fills  out,  the  heating  to  the  surface 
approaches  the  turbulent  heating  range.  Photographic 
examples  of  "wave-like”  instabilities  which  preceded  the 
large-scale  breakdown  into  turbulence  are  shown  in  Figure 
2. .  In  many  cases,  wave  trains  would  be  a  more  accurate 
description  for,  as  shown  in  Figure  25,  there  were  many 
instances  where  these  instabilities  did  not  undergo  a 
significant  change  in  structure  as  they  were  converted 
along  the  cone.  High-speed  photographs  obtained  with 
the  Cordin  camera  indicate  that,  in  fart,  these 
instabilities  "twisted"  as  they  were  converted 
downstream,  suggesting  a  helical  structure.  Following 
the  region  of  "wave-like"  instability,  the  laminar  boundary 
layer  began  to  develop  large-scale  instabilities  at  an  edge 
(see  Figure  26);  then,  transition  began.  Photographic 
examples  of  the  development  of  the  transition  region 
following  the  "wave-like"  instabilities  discussed  in  the 
previous  paragraph  are  shown  in  Figures  26a  and  b.  In 
Figures  26c  and  26d,  we  see  the  gross  instabilities 
developing  as  they  were  swept  downstream,  culminating 
in  some  very-large-scale  turbulent  bursts.  These  dis¬ 
turbances  traveled  down  the  cone  at  convective  velocities 
between  0.5  and  0.8  of  the  freestream  velocity.  The 
gross  size  of  the  eddies  generated  in  this  region  is  of 
the  same  order  as  the  boundary  layer  thickness.  In 
Figures  27a  and  27b,  are  shown  examples  of  the  gross 
mixing  in  the  transition  region  which  precedes  the 
development  of  the  turbulent  boundary  layer.  These 
photographs  show  the  "wave-like"  instability,  followed  by 
by  the  gross  instability,  followed  by  a  region  in  which 
there  is  large-scale  turbulent  mixing.  In  the  latter  region, 
it  appears  that  the  scale  of  the  turbulence  exceeds  that 
in  the  following,  fully  turbulent  boundary  layer;  it  is  in 
this  region  that  the  heat  transfer  overshoot  occurs.  In 
a  large  number  of  studies  in  hypersonic  flows,  it  takes 
many  boundary  layer  thicknesses  before  the  scale  sized 
of  the  turbulence  generated  in  the  transition  process  are 
no  longer  of  importance.  Figure  28  shows  a  compilation 
of  measurements  made  of  the  velocity  exponent  (n)  in 
the  relationship  ^  =  (-£•)  'ft  obtained  downstream  of 
the  transition  region  in  hypersonic  flow.  Immediately 
downstream  of  transition  n  reaches  a  peak,  reflecting 
the  existence  of  a  predominance  of  large  scale  eddies 
in  the  boundary  layer.  These  eddies  break  into  smaller 
eddies  as  they  are  transported  downstream.  However 
transition  lengths  from  50  to  100  boundary  layer  thickness 
are  required  before  the  turbulent  boundary  layer 
approaches  equilibrium. 

VI.  TRANSITIONAL  AND  INTERACTING  FLOWS 
OVER  NOSETIPS  AND  LEADING  EDGES 

The  occurrence  of  boundary  layer  transition  and  the 
subsequent  turbulent  flow  in  stagnation  regions  or  on  the 
attachment  line  of  a  highly  swept  leading  edge  of  a  wing 
or  intake  will  have  important  impact  on  the  design  and 
aerothermal  performance  of  hypersonic  vehicles.  In  fact 
the  trajectory  flown  by  such  vehicles  may  be  selected 
in  part  on  the  basis  of  transition  criteria.  In  general 
the  application  of  *he  measurements  of  the  occurence 
and  mechanism  of  boundary  layer  transition  made  in 
hypersonic  facilities  to  predict  transition  on  flight 
vehicles  has  in  fact  met  with  little  success  except  where 
transition  is  tripped  by  roughness  or  a  disturbance  induced 
on  the  surface  of  the  body.  Recently,  however,  it  has 
been  shown  by  Poll^,^,  and  others*^  that  transition  on  the 
surfaces  of  flight  vehicles  like  the  Space  Shuttle,  with 


highly  swept  leading  edges,  may  be  predicted  with  good 
accuracy  from  transition  measurements  made  in  ground 
test  facilities.  One  explanation  for  the  success  of'  the 
wind  tunnel  based  predictions  in  this  instance  may  lie  in 
understanding  that  the  underlying  mechanisms  of 
transition  is  a  "by-pass"  disturbance  that  is  generated  on 
he  body  rather  than  external  to  it.  For  models  whose 
svface  roughness  generates  disturbances  large  enough  to 
exceed  those  generated  externally,  transition  may  be 
controlled  principally  by  the  unit  Reynolds  number  of 
the  freestream  and  the  modcis  geometric  characteristics 
rather  than  those  of  the  tunnel  in  which  the 
measurements  were  made. 

While  there  have  been  a  number  of  correlations 
suggested  over  the  years  to  correlate^-''*  the  occurrence 
of  boundary  layer  transition  on  rough  surfaces,  the  most 
recent  work  by  Batt  and  Legner5^  appears  to  be  one  of 
the  most  comprehensive  and  useful.  As  shown  in  Figure 
5  the  occurrence  of  transition  was  found  to  be  best 
correlated  in  terms  of  a  Reynolds  number  based  on 
momentum  thickness  at  the  start  of  transition  and  non- 
dimensional  roughness  height  «'/{*'  ,  where  fc'=kO«35o«</e„) 
and  Wre  •  This  correlation  is  supported 

by  range  measurements  made  by  Reda1*”  who  however 
suggests  a  transition  criteria  in  terms  of  a  roughness 
Reynolds  number  fiK  i/K  k/m^  ,  based  on  the  local 
properties  at  the  edge  of  the  roughness.  Using  either 
criteria  it  is  easily  shown  that  for  a  significant  part  of 
the  trajectories  proposed  for  a  transatmospheric  vehicle, 
surface  irregularities  of  the  order  of  I/I 000  inch  would 
trip  transition. 

As  shown  by  Poll 5£f,  the  junction  between  a  highly 
swept  wing  or  intake  and  forebody  can  act  as  a  large 
trip  forcing  transition  to  turbulence  on  the  attachment 
line  of  a  leading  edge  at  relatively  low  freestream 
Reynolds  numbers.  For  these  flows  it  has  been  found 
that>?„-(iyUe/d>d/v*)'”  is  a  characteristic  dimension,  and 
the  crossflow  Reynolds  number  q»  o/v  =  no  at  which 
transition  occurs  is  independent  of  Mach  number  (sec 
Figure  29).  Clearly  further  studies  in  hypersonic  flows 
are  required  to  examine  tire  validity  of  tins  relationship 
and  determine  Mach  number  and  wall  cooling  effects. 
However,  as  shown  in  Figure  39,  the  Mach  number 
dependence  suggested  by  this  relationship  is  not 
inconsistent  with  the  transition  measurements  on  sharp 
cones  and  rough  nosetips  in  hypersonic  flow.  One  of  the 
most  severe  aerothermal  loading  conditions  could  arise 
if  a  strong  shock,  from  an  external  component  or  from 
an  engine  or  pod  mounted  on  the  swept  wing,  is  incident 
on  the  turbulent  boundary  layer  over  tire  attachment  line 
as  discussed  in  Section  VIII.  Currently  the  only  measure¬ 
ments  available  in  regions  of  shock/shock  interaction 
have  been  made  for  laminar  flow  over  smooth  surfaces. 
As  in  the  case  of  shock/shock  interaction  on  the  indented 
nose  shapes  discussed  in  the  following  section,  we  believe 
that  tiro  largest  heat  loads  are  generated  in  turbulent 
interaction  regions  over  rough  surfaces. 


Shock/Boundary  Laver  Interactions  on  Indented 
Nosetips 

Measurements  made  in  full-scale  flight  tests  and 
ground  test  facilities  with  ablating  nosetips  have 
demonstrated  that  under  many  re-entry  conditions 
concave  or  indented  noseshapes  can  be  developed5’  as 
the  nosetip  ablates.  The  increased  heating  resulting  from 
boundary  layer  transition  close  to  the  sonic  region  is 
responsible  for  the  initial  development  of  the  indentation. 
Subsequently,  the  large  heating  rates  which  are  developed 
in  the  recompression  region  resulting  from  shock- 
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blowing-  are  not  additive,  in  most  codes  tire  magnitude 
of  rough-wall  heating  and  surface  blowing  are  calculated 
by  essentially  superimposing  calculations  based  bn  models 
of  these  two  effects  considered  separately.  In  these 
codes  the  assumption  that  surface  roughness  plays  the 
primary  role  in  controlling  heating, augmentation  may  be 
significantly  in  error  when  applied  to  highly  ablating 
surfaces.  If,  as  we  deduce  from  the  experimental  studies, 
surface  blowing  effectively  smooths  the  rough  surfaces, 
then  selection  of  an  effective  roughness  height  as  the 
single  correlating  parameter  connecting  flight 
measurements  with  code  predictions  may  be  in  error. 
Also,  since  most  predictive  techniques  employ  an 
effective  sand-grain  roughness  as  the  single  length  scale 
characterizing  roughness  size,  there  continues  to  remain 
a  key  problem  in  defining  the  surface  topography  of  a 
rough  hypersonic  vehicle. 


however,  no  single  parameter  or  combinations  of 
parameters  (e.g.,  Rc^K/O  )  has  been  used  with  any 
great  success  to  describe  the  general  similitude  of 
turbulent  boundary  layers  in  supersonic  and  hypersonic 
flows  over  rough,  highly  cooled  walls. 

The  studies  of  Dvorak^**,  Bcttermann^,  Lewis*’*’, 
Simpson**^  and  more  recently  Lin”®  and  Finson”°  have 
provided  further  insight  into  the  basic  effects  of 
roughness  shape  and  spacing  on  the  characteristics  of 
the  rough  wall  boundary  layer  and  skin  friction  and 
heating  to  a  rough  surface.  Dvorak  combined  the  effects 
of  roughness  shape  and  spacing  into  a  single  parameter  A 
(the  roughness  density),  defined  as  shown  in  Figure  40. 
He  linked  the  downward  shift  in  the  velocity  profile 
&ulu?  ^  to  a  combination  of  roughness  Reynolds 
number  RcK(uf  k}Vuj)  and  A  through  the  relationship 


The  experimental  studies  of  Nikuradse^  and 
Schllchtlng^,  both  hydraulic  pipe  flow  studies,  were 
principally  responsible  for  the  selection  of  sand-grain 
roughness  as  the  standard  against  which  to  measure 
relative  effects  of  other  types  of  roughness.  Although 
this  standard  has  been  frequently  employed,  the  topo¬ 
graphical  characteristics  of  a  sand-grain  surface  have 
yet  to  be  defined.  In  fact,  because  of  the  experimental 
difficulties  involved  in  the  preparation  and  inspection  of 
rough  surfaces  inside  a  small-diameter  pipe,  it  is 
surprising  that  Nlkuradse's  results  are  as  consistent  as 
reported.  For  fully  developed  pipe  flows,  Nikuradse 
established  that  the  parameter  controlling  the  similitude 
of  the  flows  is  the  roughness  Reynolds  number  (veK Ify). 
This  parameter  (gen)  was  selected  in  many  subsequent 
studies  to  characterize  boundary  layer  flows  where  other 
non-dlmenstlonal  groupings  (such  as  ks/f*,  Ks/«  ,  etc.) 
might  have  been  considered  more  valid.  The  Schlichting 
studies,  conducted  with  roughness  of  well-defined 
geometric  shapes,  provided  the  first  set  of  measurements 
which  could  be  reproduced  in  both  experimental  and 
theoretical  studies.  The  results  from  these  studies, 
together  with  those  from  a  number  of  subsequent 
Investigations  in  subsonic  adiabatic  flows,  were  correlated 
by  researchers  to  yield  relationships  between  an 
"effective  sand-grain  height"  and  parameters  which 
describe  the  geometric  features  of  the  surface— a  step 
which  further  perpetuated  the  use  of  sand-grain  roughness 
as  a  standard.  The  Dirling  correlation,  see  (Figures  39 
and  40)  is  one  such  plot,  from  which  an  effective  sand- 
grain  roughness  height  can  be  determined  from  knowledge 
of  peak-to-valley  roughness  height  together  with  the 
shape  and  spacing  of  the  roughness  elements.  While 
there  is  little  direct  supporting  experimental  evidence, 
the  "lambda"  form  of  the  correlating  curve  is  assumed 
to  reflect  a  sudden  change  in  flow  structure  from  an 
"open"  to  a  "closed"  cavity  flow  around  (between)  the 
roughness  elements  as  the  spacing  between  the  roughness 
elements  is  varied.  Do  open  and  closed  cavity  flows 
really  exist  on  rough  surfaces  constructed  from  three- 
dimensional  roughness  elements?  Is  not  the  structure  of 
the  flow  around  the  elements  also  dependent  upon  the 
local  Reynolds  number,  a  parameter  not  taken  into 
account  in  any  such  correlations?  If  an  effective 
roughness  height  can  be  accurately  determined  from  a 
'bump  curve',  this  dimension  must  be  combined  with  key 
fluid  dynamic  properties  to  yield  a  non-dimensional 
parameter  or  groups  of  parameters  with  which  to 
characterize  the  flow.  The  roughness  Reynolds  number 
Se^CVf  k/v^)  originally  used  by  Nikuradse,  and  the 
non-dimensional  roughness  heights  K!S*,^9m  and  fc/0r 
(whereS*,^  ,  and  Sr  are  the  displacement,  momentum, 
and  thermal  energy  thickness  respectively),  have  all  been 
used  to  correlate  the  aerothermal  effects  associated  with 
boundary  layers  over  rough  re-entry  vehicles.  To  date, 
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In  incompressible  flows  the  smooth  regime,  where  the 
surface  shear  is  entirely  due  to  viscous  shear  is  defined 
by  kU(.ji'<r> .  At  larger  ,(s< ReK  <  7o)  ,  the 

surface  shear  !c  composed  of  form  drag  on  the  roughness 
elements  combined  with  viscous  shear.  For  K-~  >io  , 
the  surface  shear  results  principally  hum  drag,  ai.d 
viscosity  is  no  longer  factor  in  controlling  the  velocity 
profile.  For  tullv  rough  flown,  liquation  I  can  be 
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which  Is  a  function  of  the  roughness  density.  Here  it 
should  be  noted  that  (*(/))* A  {//where  0%  is  the  velocity 
close  to  the  top  of  the  roughness  elements  and  A(^)  is 
the  smooth  wall  constant.  M«w  the  definition  of 
Nlkuradse's  sand-grain  roughness  is  basically 

17*  -  X  (kT)  4  . 3 

Hence,  combining  equations  2  and  3  we  obtain 
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which  is  generally  correlated  by  two  relationships.  The 
existence  of  two  regions  has  been  rationalized  on  physical 
grounds,  supported  by  experimental  studies  in  two- 
dimensional  flows.  There  is  serious  debate  about  whether 
a  significant  change  in  flow  structure  occurs  as  three- 
dimensional  roughness  elements  are  drawn  together. 
Certainly  the  measurements  with  stone  roughness  do  not 
exhibit  such  a  trend.  This  Is  particularly  unfortunate, 
since  Nikuradse's  data  falls  on  the  line  constructed 
principally  from  data  obtained  on  roughness  constructed 
by  two-dimensional  machined  grooves. 


Finson's  engineering  model  based  on  the  early 
concepts  of  Liepmann  and  Goddard^**  and  his  (Finson's) 
detailed  numerical  calculations  provides  a  good  basis  for 
interpreting  tire  physical  phenomena  of  key  importance 
in  rough  wall  heating,  as  well  as  a  relatively  simple 
prediction  technique.  The  shear  on  a  rough  wall  can  be 
expressed  as  the  sum  of  the  viscous  and  form  drag  of 
the  rough  surface:  rr  *  , 

where  BCp)  is  the  blockage  factor,  and  d(y)  and  0  are 
the  diameter  of  the  roughness  clement  and  the  spacing 
between  elements,  respectively.  From  his  detailed 
numerical  solutions,  Finson  showed  that  p  and  u.  were 
relatively  constant  between  the  base  and  top  of  the 
roughness  element  at  values  close  to  the  top  of 

the  roughness,  and  this  equation  becomes 
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where  A^,/ A9  is  the  ratio  of  projected  area  of  the 
roughness  element  in  the  direction  of  the  flow  to  total 
area  of  the  flow  or,  which  they  stand,  and  e  (•£•)  is-  the 
average  value  of  BCy)  . 

For  compressible  flows  Finson  found  that 

Vk/V,  =  O.  U7  +.Z34  fa  f AL ) 

*  \f>e  Af) 

Therefore,  assuming  Cf6  is  the  smooth  wall  heating 
level,  It  Is  possible  to  relate  the  rough  wall  skin  friction 
to  the  smooth  wall  value  in  the  generalized  form 

zfcZ  -  fa ■(*•  §  M  x  ®.  ( f ). 

ana  the  relationship  for  heat  transfer  is  assumed  of  the 


same  form. 
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If  we  assume  that  the  product  o!  the  blockage  factor 
and  Cp  are  invarient  with  roughness  shape  and  space, 
then  for  constant  local  free  stream  condition,  we  get 
the  Dvorak-Simpson  Parameter: 

£i!  <n>uin  _  p  / ) 
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A  slightly  different  form  can  be  obtained  by  the  subsonic 
blunt  body  approximation 

c*  =  CD  xer  (Af/A»s) 

and  using  d]k  rather  than  fa/Ap  to  obtain  the 
correlation  In  terms  of  the  Dirling71  parameter, 
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where  the  various  areas  are  illustrated  in  Figure  39. 


momentum  and  energy  exchange  in  high-speed  flows  over 
adiabatic  and  noh-adiabatic  surfaces. 


Voisiriet's6^  studies  of  the  combined  effects  of 
roughness  and  blowing  were  conducted  at  Mach  6  under 
adiabatic  wail  conditions.  The  measurements  made  in 
these  studies  have  clearly  demonstrated  that  the 
combined  effects  of  blowing  and  roughness  on  skin 
friction  cannot  be  described  in  a  simple  manner.  The 
effects  of  surface  roughness  alone  on  skin  friction  were 
shown  to  be  correlated  in  terms  of  the  roughness  Reynolds 
number  gen.',  a  result  consistent  with  earlier  measurements 
on  adiabatic  walls  by  Goddard'  ®  and  Reds7  .  In  contrast, 
measurements  oh  model  placed  in  the  flow  involving 
significant  levels  of  heating  (  /h0  *  0.5)  have  in 

general  correlated  better  with  parameters  like  K/S,  K/<? or 
KjSr  which,  as  shown  by  the  theoretical  studies  of  Dvorak 
and  Finson,  should  have  greater  relevance  to  roughness 
effects'  on  re-entry  vehicles.  Voisinet's  studies 
demonstrated  that  the  effects  of  surface  roughness  and 
blowing  on  skin  friction  cannot  he  deduced  from  simple 
expressions  derived  from  the  measurements  made  of  each 
of  the  separate  effects. 


In  recent  studies  of  effects  of  roughness  shape  and 
spacing  on  the  heat  transfer  and  skin  frlclion  to  the 
roughness,  nose  tips,  frusta  and  flaps  of  a  typical  MRV 
configuration,  Holden61  used  both  the  Dvorak/Shnoson 
parameter  As/Af,  and  the  Dirling  parameter/^  (^*?and 
achieved  reasonable  success  in  correlations,  In'lhese 
studies,  the  effects  of  roughness  shape  and  spacing  on 
the  heat  transfer  and  skin  friction  both  for  surfaces  with 
sand-grain  roughness  and  those  corjstructed  with 
geometrically  well-defined  hemispherical  and  conical 
roughness  elements  were  examined.  Also,  heat  transfer, 
skin-friction,  and  pressure-distribution  measurements 
were  obtained  on  spherical  and  ablated  noseshapes, 
conical  frusta,  ar,d  the  eontrol  surfaces  of  MRV  vehicle 
in  hypersonic  flow.  Typical  measurements  and  corre¬ 
lations  are  shown  in  Figures  41,  42,  43  The  results  ol 
this  Work  have  demonstrated  that  the  low-speed 
measurements  of  Nikuradsc62,  Schiichting63,  and  others, 
and  the  correlations  of  Dirling7 '/Simpson67,  cannot  be 
wed  directly  to  predict  rough-wall  heating  and  skin 
friction  in  supersonic  and  hypersonic  flows  over  non- 
adiabatic  surfaces.  Even  the  more  solidly  founded 
prediction  scheme  developed  by  Finson6^.  which  is  based 
upon  his  detailed  numerical  solutions,  consistently 
overpredicts  the  roughness-enhanced  heating  levels  in 
high-speed  flows.  While  direct  measurements  of  the  skin 
friction  of,  and  heat  transfer  to,  geometrically  well- 
defined  rough  surfaces  provide  the  opportunity  to  more 
closely  evaluate  the  accuracy  of  the  current  shape-change 
codes  and  the  more  fundamental  treatments  like  that  of 
Finson  in  high  Mach  number,  high  Reynolds  number  flows, 
there  remains  a  basic  need  to  define  more  closely  the 
fundamentals  of  the  fluid  dynamics  which  control 


The  measurements  made  in  studies  with  a 
transpiration  cooled  model  with  and  without  blowing  is 
shown  in  Figure  44  together  with  measurements  of  tire 
heating  rates  to  rough  and  smooth  non- porous 
hemispherical  models  at  the  same  condition.  The 
measurements  on  the  model  with  zero  blowing,  shown  in 
Figure  44,  clearly  show  that  the  intrinsic  roughness  of 
the  surface  causes  heating  enhancement  factors  of  over 
1.7.  fn  fact,  it  can  be  seen  by  comparing  Figures  44 
and  45  that  the  heat  transfer  measurements  on  the  rough 
hemisphere  are  In  good  agreement  with  those  obtained 
on  the  non-blowing  transpiration-cooled  nosetip  model 
such  that  they  would  be  biased  toward  tbe  larger  beating. 
When  a  small  amount  of  blowing  (m//5u.ew=  0.032)  was 
introduced,  the  heating  rates  ovc-r  a  major  part  of  the 
transpiration-cooled  model  dropped  to  levels  dose  to 
those  recorded  on  the  smooth  model,  as  shown  in  Figure 
44,  It  could  be  postulated,  on  the  basis  of  those 
measurements,  that  the  initial  effect  of  mass  addition 
from  a  rough  ablating  nosetip  is  to  modify  the  flow 
around  the  roughness  elements,  by  eliminating  the  cavity 
flows,  in  such  a  way  that  the  momentum  defect 
introduced  by  tlse  roughness  is  small  as  illustrated  in 
Figure  46.  How  does  roughness  shape  and  spacing 
influence  the  heat  transfer  to  the  surface  in  the  presence 
of  mass  addition?  Are  the  current  relationships  which 
have  been  developed  Irom  an  essentially  subsonic, 
adiabatic  flow  data  base,  to  describe  roughness  heating 
on  blunt  nosetips  applicable  to  the  hypersonic  flow 
transpiration-cooled  vehicles?  Holden's  studies  of  rough¬ 
ness  shape  and  spacing  effects  on  on-ablating 
configurations  suggest  that  not  only  are  the  subsonic 
studies  inapplicable  to  tire  heating  of  heat  shields  in 
hypersonic  flow,  but  further  that  the  basic  modeling  of 
the  roughness  drag  and  mechanisms  cf  heating  in  the 
theoretical  models  is  highly  questionable.  In  these  studies 
it  was  shewn  that  even  in  mildly  supersonic  How, 
shocklets  are  clearly  evident  (see  Figure  47)  around  each 
roughness  element.  Thus  not  only  is  the  model  of  drag 
incorrect,  the  entropy  layer  associated  with  each  shocklGi 
contributes  to  a  higher  temperature,/  lower  momentum 
How  around  the  roughness  elements  than  described  in 
the  current  roughness  models. 
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VIII.  SHOCK.  WAVE/TURBULENT  BOUNDARY  LAYER 
INTERACTIONS  IN  HIGH  SPEED  FLOWS 

Two-Dimensional  Shock-Wave/ 

Turbulent  Boundary  Layer  Interactions 

The  development  of  two-dimensional  turbulent 
separated  regions,  induced  in  a  compression  corner,  and 
at  the  base  of  an  incident  shock,  with  increase  in 
interaction  strength  in  hypersonic  flow  as  obtained  by 
Holden71*,  are  illustrated  in  Figures  49  and  50,  Separation 
Is  first  observed  in  the  laminar  sublayer  and  a  well- 
defined  separation  bubble  is  clearly  visible  in  Figure  49. 
The  initial  development  of  the  separation  region  takes 
place  by  an  elongation  into  the  laminar  sublayer,  with 
the  separation  and  reattachment  shocks  combining  within 
the  boundary  layer  to  form  a  single  shock.  Only  when 
the  separation  point  has  fed  well  forward  of  the  junction 
is  a  well-defined  plateau  region  formed.  Then,  in  contrast 
to  laminar  Interaction  regions,  the  separation  shock 
originates  at  the  bottom  of  the  boundary  layer  and  is 
contained  within  the  boundary  layer  until  it  coalesces 
with  the  reattachment  compression  process.  In  separated 
regions  induced  by  an  externally  generated  shock, 
separation  first  takes  place  in  the  region  where  the 
incident  shock  strikes  the  laminar  sublayer,  see  Figure 
50b.  The  separation  point  moves  forward  with  Increasing 
strength  cf  the  incident  shock  until  the  separation  shock 
becomes  visible  in  the  invlscid  flow  downstream  of  the 
incident  shock!  as  yet,  separation  is  still  downstream  of 
the  point  where  the  incident  shock  passes  through  the 
edge  of  the  boundary  layer.  For  large  incident  shocks, 
boundary  thickening  occurs  ahead  of  the  incident  shock 
in  an  analogous  fasion  to  laminar  flow  separation. 
However,  as  in  wedge-induced  separated  regions,  viscous- 
lnvlscid  interaction  takes  place  almost  entirely  within 
the  original  boundary  layer.  The  structure  wedge  and 
shock-induced  turbulent  interaction  regions  at  Mach  !3 
are  very  similar  to  those  at  Mach  8;  however,  as  we 
might  anticipate,  the  viscous  interaction  region  and  the 
associated  shocks  are  even  more  firmly  embedded  within 
the  original  boundary  layer. 

Surface  measurements  with  high  frequency 
instrumentation  indicated  that  turbulent  separated 
regions  were  highly  unsteady,  and  typically  tire  separation 
point  would  oscillate  in  a  streamwlse  direction  with  an 
amplitude  of  approximately  one-quarter  to  one-third  of 
the  local  boundary  layer  thickness,  at  frequencies  in  the 
range  from  1  to  10  kHz  (see  Figure  51),  The  unsteady 
character  of  the  records  from  transducers  in  the 
recirculation  region  indicated  it  could  be  unrealistic  to 
assume  that  a  laminar  sublayer  model,  in  the  conventional 
sense,  could  be  used  to  describe  the  lower  part  of  the 
recirculating  region  as  is  done  in  some  of  the  triple  deck 
calculations.  The  mean  distribution  of  skin  friction,  heat 
transfer  and  pressure  to  the  wall's  bounding,  both  shock- 
and  wedge-induced  interaction  region  were  similar  for 
well-separated  flows  with  identical  total  pressure  rises, 
both  the  pressures  and  heat  transfer  distributions  are 
characterized  by  well-defined  plateaus  in  the 
recirculation  region  and  large  gradients  in  the  separation 
and  reattachment  regions.  The  maximum  heat  transfer 
rates  generated  in.  the  reattachment  regions  of  these 
flows  is,  of  course,  of  considerable  importance.  Holden75 
found  that  for  separated  interaction  regions  the  maximum 
pressure  and  heat  transfer  measurements  over  the  Mach 
number  range  from  2  to  13  could  be  correlated  in  the 
form  shown  in  Figure  52, 

The  influence  of  Reynolds  number  on  the  size  of 
two  dimensional  separated  regions  remains  unresolved. 


Most  of  the  early  studies  of  shock  wave-turbulent 
boundary  layer  interaction  were  made  in  the  turbulent 
boundary  layer  over  a  tunnel  wall.  Major  discrepancies 
were  found  between  experimental  facilities;  for  example, 
the  measurements  of  Bogdonoff  and  Kepler76  differed 
considerably  from  those  of  Gadd77  for  identical 
freestream  Mach  numbers  and  interaction  strength,  when 
both  experimenters  had  indicated  that  there  was  little 
effect  of  freestream  Reynolds  number  on  the  length  of 
the  separated  region.  However,  the  measurements  of 
Green7**,  Roshko  and  Thoinke7^,  Law®61,  Settles, 
Bogdonoff  and  Vas°*j  and  Appels*7,  all  made  under 
adiabatic  wall  conditions,  indicate  that  increasing 
Reynolds  number  decreases  the  size  of  a  turbulent 
separated  region.  In  contrast,  the  studies  of  Chapman, 
Kushn  and  Larson**5,  Kudin**1*,  and  Holden71*,  Elfstrom*  ’ 
and  Appels**6  at  hypersonic  speeds,  all  conducted  on 
"highly  coded"  models  mounted  in  the  test  section,  have 
shown  the  opposite  trend.  As  In  the  case  of  incipient 
separation,  the  answer  may  lie  in  changes  in  the 
equilibrium  structure  of  a  turbulent  boundary  layer  with 
Reynolds  number.  This  hypothesis  Is  supported  by  present 
measurements  of  shock-  and  wedge-induced  separated 
flows  at  Reynolds  numbers  of  up  to  3  x  I07  by  Holden75. 

In  the  experimental  studies  by  Holden,  where 
dynamic  measurements  of  the  skin  friction  In  turbulent 
interaction  regions  were  made,  the  separation  conditions 
were  defined  when  the  time  average  of  surface  shear  at 
one  point  only  on  the  surface  was  zero  (see  Figure  51). 

Comparisons  with  Navier-Stokes  Solutions 

The  complexity  of  the  flow  field  in  regions  of  shock 
wavc/turbulont  boundary  layer  interaction  Is  such  that  it 
is  unrealistic  to  expect  to  describe  such  regions  In  any 
detail  within  the  framework  of  the  boundary  layer 
equations.  Indeed  their  are  some  who  would  question 
whether  the  time  or  mass  averaged  NavHr-Stokcs 
equations  capture  the  basic  fluid  mechanics  associated 
with  the  intrinsically  unsteady  nature  of  separated 
regions.  In  hypersonic  flows  the  effects  of 
compressibility  on  the  structure  and  development  of 
turbulence  must  also  be  considered. 

While  there  have  been  strenuous  efforts  to  obtain 
predictions  of  2D  and  3D  turbulent  interaction  regions, 
it  is  currently  recognized  that  "successes"  with  "Navier- 
Stokes"  code  in  describing  some  3D  turbulent  interactions 
regions  ii  a  result  of  the  dominance  of  the  pressure  and 
inertial  terms  in  these  flows.  In  these  latter  comparisons 
(References  88  &  89)  St  was  found  that  the  modeling  of 
turbulence  could  be  changed  without  significantly 
changing  the  numerical  solution.  For  2D  interactions  it 
appears  the  modeling  of  turbulence  is  more  critical.  To 
obtain  good  agreement  for  these  latter  flows  some  very 
gross  assumptions  must  be  made  in  the  turbulence  model. 
Shang  and  Hankey55,  for  example  chose  to  apply  an 
empirical  relationship  (selected  by  matching  the  length 
of  the  separated  region)  to  rapidly  decrease  the  turbulent 
scale  size  through  the  interaction  region  as  shown  in 
Figure  53.  Horstmann5*,  however,  found  the  best 
agreement  with  Settles52  measurements  in  wedge-induced 
separated  regions  using  a  two  equation  model  (or 
turbulence  scale  size  and  vorticity,  as  shown  in  Figure 
54.  Working  with  this  same  turbulence  model,  however, 
Horstmanr.52  was  unable  to  predict  the  occurrence  of 
separation  on  two  incident  shock/turbulent  boundary  layer 
configurations  studied  by  Holden55  at  Mach  11.2.  As 
shown  in  Figures  55  and  56  botli  these  flow  fields  are 
clearly  separated  and  yet  the  numerical  solution  fails  to 
predict  the  characteristic  plateaus  in  either  the  heat 
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transfer  or  pressure  distributions.  The  modeling  of 
turbulence  in  separated  interaction  regions  at  hypersonic 
Mach  numbers  .should  account  for  the  effects  of 
compressibility,  and  the  .generation  of  Turbulence  by  the 
unsteady  movement  of  the  incident  and  induced  shocks 
as  they  traverse- and  interact  with  a  major  region  of  the 
turbulent  boundary  layer.  Clearly  further  detailed 
experimental  work  and  insightful  theoretical  modeling 
are  required  to  develop  numerical  prediction  techniques 
which  are  capable  of  describing  turbulent  interaction 
regions  in  detail 


Three-Dimensional  Shock  Wave/Turbulent  Boundary 
Layer  Interaction 


Many  of  the  conceptual  problems  associated  with 
the  use  of  the  boundary  layer  equations  to  describe 
separated  regions  induced  by  shock  wavc/turbulent 
boundary  layer  interaction  are  circumvented  by  the  direct 
solution  of  the  Navier-Stokes.  equations.  However,  in 
their  place  we  find  tiie  equally  thorny  problem  of 
specifying  a  detailed  model  of  turbulence  for  flows  with 
exceedingly  large  strcamwise  pressure  gradients.  Despite 
the  lack  of  success  In  developing  credible  turbulence 
models  for  two-dimensional  Interaction  regions,  or 
perhaps  because  of  it,  three-dimensional  turbulent 
interaction  regions  have  become  the  focus  of  attention 
of  the  Navier-Stokes  solvers. 


Comer  Interaction  Studies 


The  axial  corner  flow  or  swept-shock  Interaction  has 
been  one  of  the  principal  configurations  selected  to 
Investigate  three-dimensional  regions  of  shock 
wave/boundary  layer  interaction.  The  swept-shock,  which 
Is  generated  by  a  wedge  or  fin  mounted  perpendicular  to 
a  flat  plate,  Impinges  normally  onto  the  flat  plate 
boundary  layer.  The  initial  studies  in  this  area  by 
Stalker9**  and  Stanbrook95  were  followed  by  the  more 
detailed  Investigations  of  McCabe96,  Peake  and 
Rainblrd97,  Oskam  et  al98,  Cousteix  and  Houdeville99, 
Dolling  and  Bcgdonoff  100»  10  *,  Dolling  and  Murphy10’, 
and  Dolling*03.  The  latter  extensive  series  of  studies 
was  conducted  at  Mach  3  under  adiabatic  wall  conditions. 
While  Incipient  separation  is  relatively  easy  to  define  for 
two-dimensional  turbulent  interactions,  this  concept  has 
generated  considerable  controversy  In  three-dimensional 
flows.  While  McCabe96  suggests  that  separation  should 
be  defined  on  the  basis  of  converging  steamlines, 
Stanbrook93  and  others  have  used  criteria  based  on  the 
inflection  points  in  the  pressure  distribution.  The 
occurrence  of  separation  was  correlated  in  simple  terms 
by  Korkegi109,  who  found  that  in  low  Macli  number  flow, 
deflection  angle  g  w;  for  incipient  separation  varies  as 
the  Inverse  of  the  upstream.  Mach  number,  i.e.,  = 

0.3Mq,  while  for  2  <  M  <  3.4  Korkegi  suggests  that  pi/p 
is  Independent  of  Mach  number.  Goldberg's103  anq 
Holden's106  measurements  at  Mach  6  and  11  respectively 
do  not  agree  with  the  Korkegi  correlation. 


Studies  with  the  emphasis  on  the  heating  in  swept- 
shock  interaction  regions  have  been  conducted  by 
Neumann  and  Burke*07.  Law108,  Token109,  and 
Scuderi110  and  Hqlden111.  Figure  57  shows  typical 
distributions  of  heat  transfer  and  pressure  along  a 
streamwise  cut  through  the  interaction  region  together 
with  nomenclature  which  is  in  conventional  use.  While 
the  heat  transfer  and  pressure  distributions  exhibit  a 
uniform  and  monotonic  increase  through  attached 
interaction  regions,  when  the  flow  separates,  distinctive 
plateau  regions  are  formed  in  the  heat  transfer  and 
pressure  distributions,  as  depicted  in  Figure  57.  As  noted 
above,  at  low  Mach  numbers  (M  =  2-?4)  and  for  adiabatic 


surfaces}  a  large  body  of  data  exists  on  tile  mean 
characteristics  of  swept-shock  interactions.  Strangely, 
this  body  of  3D  o'ata  has  been  found  to  be  in  better 
overall  agreement  with  the -Hung  arid  MacCormack1 17, 
Horstmann91,  Shang  and  Hankey90,  Settles  and 
Horstmarin97  solutions  to  the  Navier-Stokes  equations 
than  the  relatively  less  complex  two-dimensional  flow- 
separation  over  a  flat  plate/ wedge.  Thc.ie  results  are 
not  as  sensitive- to  the  turbulence  model  and  suggest  that 
the  gross  features  of  the  flows  are.  controlled  principally 
by  invsicid- effects. 


In  a  corner  flow,  tire  swept-shock  generated  by  tire 
inclined  fin  impinges  on  the  turbulent  boundary  layer  in 
a  plane  perpendicular  to  the  fiat  plate.  The  basic 
mechanism  of  pressure  rise  through  die  interaction  is 
therefore  controlled  principally  by  the  component  of 
freestream  Mach  number  normal  to  this  shock  (M_  Sin£). 
A  highly  simplified  visualization  of  the  viscous/inviscid 
interaction  with  flow  separation  is  sketched  in  Figure 
58.  Here,  we  consider  the  flow  in  tire  plane  norma!  to 
the  plane  of  tire  shock  to,  be  similar  to  that  in  transonic 
flow.  When  flow  separation  occurs,  a  three-dimensional 
vortex  is  formed,  the  pressure  in  which  is  relatively 
constant  at  the  "two-dimensional"  plateau  level,  as  we 
will  show  later.  The  stroamwisc  distribution  of  heat 
transfer  in  this  region  Is  also  found  to  be  relatively 
constant,  and  Indeed  we  and  others  using  skin  friction 
and  oil  flow  measurements  have  correlated  the  first 
appearance  of  a  plateau  region  in  the  heat  transfer  with 
a  significant  change  in  the  flow  structure  which  is  linked 
with  flow  separation,  in  fact,  Token109  has  shown  that 
the  McCabe96  criteria,  based  or,  an  examination  of 
surface  oil  streaks  in  the  neighborhood  of  flow  separation, 
arc  less  sensitive  methods  for  detecting  flow  sepa,  ation 
than  observations  based  on  changes  in  the  heat  transfer 
distribution  with  increased  interaction  strength.  In  a 
recent  experimental  study  at  Mach  II,  Holden111  used 
the  Incipient  formation  of  a  plateau  in  the  heat  transfer 
distribution,  together  with  a  marked  increase  in  the 
fluctuation  levels,  in  the  output  of  the  thin  film 
instrumentation,  as  marking  the  onset  of  flow  separation. 
As  shown  in  Figure  59  Holden's  measurements  indicated 
that  iri  hypersonic  flow  over  highly-cooled  surfaces  the 
turbulent  boundary  is  more  tenacious  in  resisting  boundary 
layer  separation  than  predicted  by  the  methods  derived 
by  McCabe96  and  Korkegi109.  Holden's  measurements 
of  the  peak  pressure  ratio  through  the  interaction  and 
the  plateau  pressure  rise  are  in  better  agreement  with 
calculations  based  on  an  inviscid  flow  mode!  in  the  2D 
theory  of  Ueshotko  and  Tucker113  than  the  correlations 
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of  Scudcri MU  as  shown  in  Figure  60.  Holden  found 
that,  as  in  tire  studies  of  two-dimensional  separated 
interaction  regions,  the  peak  heating  can  be  related  to 
tire  overall  pressure  rise  by  a  simple  power  law 
relationship  as  shown  in  Figure  61.  Figure  62  shows  that 
the  maximum  pressure  rise  through  the  interaction  region 
can  be  calculated  with  good  accuracy  from  inviscid  flow 
relationships.  While  there  appears  to  be  merit  for  the 
development  of  simple  prediction  methods  in  describing 
the  flow  in  terms  of  tire  normal  flow  Mach  number,  this 
is  clearly  a  gross  oversimplification  and  it  should  be 
noted  that  tire  plateau  pressure  measurements  obtained 
in  the  current  study  were  relatively  independent  of  M0 
Sin  9  . 


Skewed  Shock/Boundary  Layer  Interaction 


Another  approach  to  exploring  flow  separation  in 
regions  of  three-dimensional  shock  wave/boundary  layer 
interaction  is  to  begin  with  a  two-dimensional  or 
axisymmettic  interaction  and  sweet)  this  interaction  (or 
introduce  angle  of  attack  for  the  axisyminetric  case)  to 
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progressively  introduce  crossflow  into  the  interaction, 
region.  Experimental  studies  of  this  type  have  been 
conducted  by  Ericsson,  Reding  and  Guenther**^,  Settles 
arid  Perkins**^,  and  Settles  drid  Teng**°.  Settles,  who 
studied  the  interaction  region  over  swept  and  unswept 
fiat  plate/wedge  configurations  in  an  adiabatic  Mach  3 
airflow,  found  that  introducing  crossflow  increased- the 
scale  of  the  separated  interaction  region.'  Considerable 
effort  was  expended  in  this  latter  study  to  determine 
the  Reynolds  number  scaling,  and  the  length  from  the 
upstream-tip  of  the  wedge  for  the  flow  to  become  quasi- 
two-dimensional.  However,  the  effect  of  changing  the 
overall  spanwise  scale  of  the  model  oh  the  scale  of  the 
interaction  was  not  examined  explicitly.  The 
measurements  of  surface  and  pitot  pressure  through  the 
Interaction  regions  were  in  good  agreement  with  solutions 
to  the  Navier-Stokes  equations  obtained  by  Horstman9*; 
however,  some  key  features  of  tire  flow  were  poorly 
predicted.  It  is  known  that  agreement  with  pressure 
data  is  not  the  most  definitive  of  tests.  More  recently 
Holden*06  performed  studies  of  crossflow  effects  on  the 
size  and  properties  of  the  interaction  region  induced  by 
a  swept-oblique-shock  incident  on  a  turbulent  boundary 
layer  over  a  flat  plate  at  Mach  11  and  Re,  =  30  x  10°. 
Experiments  were  conducted  for  two  strengths  of  incident 
shock,  the  first  (  9  SG  =  12.5°)  to  generate  a  separated 
condition  close  to  incipient  separation,  and  the  second 
(  0  SG  =  *5°)  t0  generate  a  well-separated  flow. 
Distributions  of  heat  traasfer  and  pressure  as  well  as 
schileren  photographs  of  the  unswept  or  two-dimensional 
flow  condition  and  the  30  degree  swept  condition  are 
shown  In  Figures  63  and  6 4.  It  is  clear  from  the  well 
defined  plateau  regions  in  the  distributions  of  pressure 
and  beat  transfer,  as  weli  as  the  well  defined  separation 
shock  in  the  schlieren  photograph,  that  a  well  separated 
region,  extending  two  inches  in  length,  is  induced  beneath 
the  stronger  incident  shock.  The  measurements  made  of 
the  distribuiton  of  heat  transfer  and  pressure  beneath 
the  well  separated  flow  induced  by  both  the  1 2.3  degree 
and  the  15  degree  shock  generators  swept  at  angles  of 
0  and  30  (shown  in  Figures  63  and  64)  indicated  that  the 
induced  crossflow  has  little  effect  on  the  size  and 
characteristics  of  the  interaction  regions,  if  there  is  a 
perceptible  effect,  It  is  a  decrease  in  the  length  of  the 
separated  region  with  increased  crossflow.  The 
significant  differences  between  Holden's  and  Settles92 
measurements  of  the  variation  of  interaction  length  with 
sweep  angle  and  those  obtained  in  these  studies  are  shown 
in  Figure  65.  While  Settles  finds  an  almost  threeiold 
increase  in-  separation  length  at  sweep  angles  ot  40 
degrees,  Holden  found  10  percent  reduction  in  this  length. 
Further  studies  are  required  to  resolve  this  issue. 

Three-Dimensional  Shock-Shock  Boundary  Layer 

Interaction 

The  heating  rates  generated  in  three-dimensional 
interaction  regions  by  shock-shock  interaction  can  pose 
serious  problems  for  the  designer  of  TP5.  Heating  levels 
up  to  two  orders  of  magnitude  larger  than  the  stagnation 
point  value  can  be  generated  at  hypersonic  speeds  by 
shock-shock  interaction  over  the  leading  edge  of  fins, 
inlets  and  structures  connecting  two  vehicles.  These 
regions  o!  sharply  peaked  heating  levels  are  accompanied 
by  high  pressures,  and  unlike  the  stagnation  point  their 
position  cannot  be  defined  with  ease.  Perhaps  the  most 
dramatic  example  of  damage  reuslting  from  shock 
impingement  heating  was  that  recorded  on  the  ventral 
fin  supporting  a  ramjet  model  on  the  X-I5.  At  Mach 
6.7,  the  heating  resulting  from  shock-shock  interaction 
caused  a  bum-through  in  the  fin,  ns  shown  in  Figure  66, 
which  also  graphically  demonstrates  the  high  gradients 
generated  in  such  flows. 


The:  first  detailed  investigation  of  the  fluid 
mechanics  of  the  shock-shock  interaction  was  performed 
by  Edney**2,  who  bn  the  basis  of  a  detailed  schlieren 
study  defined  six  distinct  flow  regimes  which  resulted 
from  the  orientation  of  the  shuck  systems,  three  classes 
of  interactions  are  shown  schematically  in  Figure  67 
along  with  the  typical  levels  of  enhanced  heating.  Type 
TV  interference,  in-  which  the  shock-shock  interaction 
generates  a  jet  which  impinges  on  the  surface  gives  the 
greatest  heating  enhancement.  Edney's  Schlieren 
photographs  of  this  interaction  are  shown  in  Figure  67 
for  the  shock-shock  interaction  on  a  biunt  fint.  The 
Type  IV  flow  structure  for  the  unswept  sin  configuration, 
which  is  shown  in  Figure  67,  exhibits  a  supersonic  jet, 
bounded  by  shear  layers,  that  impinges  on  the  fin  leading 
edge  resulting  in  the  highest  heating  levels.  Because 
the  shock  layer  flow  is  subsonic,  the  structure  of 
interacting  flow  is  not  readily  predicted.  Similarly  the 
shear  layer  flows  are  highly  sensitive  to  the  Reynolds 
number.  The  loading  edge  sweep  has  clearly  a  marked 
effect  on  the  interaction  heating  as  the  studies  of  Hicrs 
and  Loubsky**^  antj  Gulbran'*9  have  demonstrated 
experimentally.  The  measurements  of  heating  distri¬ 
bution  obtained  by  Hlers  and  Loubsky****  dearly  show 
the  magnitude  of  the  impingement  beating.  Hains  and 
Keyes*20  and  Keyes  and  Morris12*  expanded  on  the  Edney 
studies  to  provide  correlations  of  the  heating  levels 
encountered  in  the  various  classes  of  interaction.  These 
studies  demonstrated  a  strong  sensitivity  to  Reynolds 
numbers  and  were  followed  by  studies  by  Birch  and 
Keyes*22  who  attempted  to  correlate  conditions  for 
transition  in  these  flows,  in  general,  the  studies  that 
have  been  conducted  of  shock-shock  interaction  were 
conducted  at  relatively  low  Reynolds  numbers  for  limited 
ranges  of  Mach  number.  There  remains  a  need  for  high 
Reynolds  number  studies  In  which  detailed  measurements 
are  made  of  the  distribution  oi  heat  transier  in  fully 
turbulent  interaction  regions. 

IX.  COMPRESSIBILITY  EFFECTS  ON  TURBULENT 
MIXING 

While  in  supersonic  and  hypersonic  flows  constant 
pressure,  turbulent  boundary  layers  do  not  appear  to  be 
strongly  influenced  by  compressibility  effects,  studies  of 
the  supersonic  mixing  between  two  coaxial  streams 
suggest  that  there  is  a  significant  decrease  in  spreading 
rate  with  increased  Mach  number  *23,  Similar  effects 
have  beep  ibserved  in  the  turbulent  near  wakes  of  slender 
bodies  travelling  at  hypersonic  spccds*24,125,l26i 
Whethci  such  an  effect  is  associated  with  fundamental 
changes  in  the  mechanism  of  turbulence  production  or 
dissipation  (a  turbulence  compressibility  effect)  or  is 
associated  with  the  relative  Mach  number  or  density 
between  the  two  streams  remains  to  be  determined. 
However  such  an  effect  Itas  important  implications  for 

the  perf . ance  of  hypersonic  rain  jets  where  fuel 

(hydrogen)-air  mixing  will  occur  at  supersonic  and  perhaps 
hypersonic  speeds,  it  is  clear  from  both  the  fundamental 
and  practical  viewpoint  that  an  understanding  of  this 
phernonena  is  of  key  importance  to  me  development  of 
rational  prediction  techniques. 

In  the  ir  review  of  free  shear  layer  mixing 
measurements,  Birch  and  Eggers*2^  concluded  that  there 
is  a  significant  reduction  in  the  spreading  parameter  with 
increased  Mach  number  presenting  their  correlation  of 
tire  measurements  in  a  form  reproduced  in  Figure  68. 
More  recently  Padova,  et  al,*22  have  examined  tire 
coaxial  mixing  of  non-reacting  gas  flows  for  relative 
Mach  numbers  up  to  3  for  density  ratios  up  to  2.7;  their 
measurements  fall  between  the  Birch  and  Eggers  data 
and  those  of  Maydew  and  Rced*2S  Vvho  found  little 
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variation  of  spreading  parameter  with  Mach  number. 
Currently  there  are  a  number  of  unresolved  issues 
associated  with  the  downstream  effects  of  initial 
conditions;  however  a-  decrease  in  turbulent  spreading 
rate  has  also  been  found  in  high  Mach  number  turbulent 
wakes.  It 'has  been  observed  that  in  hypersonic  flow  the 
turbulent  wake  of  slender  bodies  over  which  the  boundary 
layer  is  fully  turbulent  grows  at  a  significantly  slower 
rate  than  a  wake  where  transition  occurs  downstream  of 
the  rear  stagnation  point.  Schlieren  photographs  suggest 
the  existence  of  small  scale  turbulence  in  the  wake 
behind  the  turbulent  body,  while  when  transition  occurs 
in  the  wake  larger  eddies  are  observed  as  sketched  in 
Figure  69..  In  his  analysis  of  "low  altitude"  wake  flows, 
Flnson  suggests  that  the  weak  growth,  of  the  large  scale 
turbulence  in  the  near  wake  of  the  turbulent  body  results 
from  a  damping  o l  turbulence  production  by  the  small 
eddies. transported  from  the  boundary  layer  over  the  body. 
Others  suggest  that  remnants  of  the  turbulent  boundary 
layer  serve  only  to  mask  observation  of  the  inner  wake 
that  they  surround,  and  the  intrinsic  difference  lies  in 
the  significantly  different  Mach  number  profile  that 
exists  across  the  wake  in  the  two  cases.  Orenbergcr 
has  suggested  that  the  occurrence  of  slow  wake  growth 
is  associated  with  a  supersonic  velocity  defect  across 
the  wake.  He  suggests  that  when  the  largest  eddy 
structures  flow  supersonically  relative  to  the  surrounding 
inviscid  wake  the  compression  waves  which  are  generated 
act  to  dampen  the  lateral  motion  of  the  eddies.  Recently 
Bushnell,  et  al.*29  have  conducted  numerical  studies 
which  demonstrate  that  the  interaction  between  an  eddy 
and  a  locaily  supersonic  flow  leads  to  edJy  breakup 
through  biflication  as  sketched  in  Figure  70.  One  may 
conjecture  thrt  a  similar  analysis  would  show  that  the 
growth  of  large  eddy  structures  would  be  inhibited  by 
analogous  mechanisms.  In  investigating  such  effects  it 
is  clear  that  we  are  attempting  to  understand  the  effects 
of  compressibility  on  turbulent  structure  and  dlffusivity. 
It  is  also  clear  that  if  such  effects  are  manifest  in  free 
shear  layer  mixing  then  they  may  also  be  important  in 
separated  flows  induced  by  shock/boundary  interaction 
and  in  flows  with  transpiration  cooling  at  hypersonic 
speeds.  This  is  an  area  where  fundamental  theoretical 
and  experimental  research  is  urgently  needed. 

X.  CONCLUDING  REMARKS 

With  the  renewed  interest  in  hypersonic  vehicles 
whose  design  requirements  are  more  complex  and 
sophisticated  than  the  essentially  ballistic  systems  which 
have  been  developed  to  date,  there  has  coine  the 
realiz?tion  that  the  current  predictive  capabilities  arc 
inadequate  in  a  number  of  important  areas.  While  there 
has  been  significant  increase  in  computing  power  during 
the  past  decade,  relatively  little  progress  has  been  made 
in  the  understanding  and  modeling  of  the  fundamental 
mechanisms  which  must  be  incorporated  into  the  machine 
codes.  During  the  past  decade  fundamental  research  in 
hypersonic  flow  has  almost  stopped,  and,  in  addition  to 
the  absence  of  work  to  develop  new  hypervelocity 
facilities,  a  number  of  hypersonic  facilities  have  been 
destroyed.  It  is  in  the  area  of  high  altitude  hypervelocity 
flight  where  real  gas,  non-equilibrium,  viscous 
interaction,  slip  and  catalytic  wall  effects  are  important, 
that  the  lack  of  ground  test  facilities  is  the  most  evident. 
While  in  the  short  term  flight  tests  will  help  in  obtaining 
measurements  in  the  hypervelocitiy  transition  regime, 
the  development  of  ground  test  facilities  in  which  clean 
reacting  airflows  at  velocities  of  up  to  20,000  ft/sec  are 
generated,  should  receive  high  priority.  Facilities  are 
also  required  iu  which  turbulent  flows  can  be  studied  on 
models  in  airflows  at  Mach  numbers  up  to  20. 


In  the  flow  regime  between  free  molecular  and 

continuum  flow  where  viscous  interaction,  non- 

equilibrium,  surface  slip  and  catalysis  are  important, 
there  is  a  real  need  to  obtain  detailed  flow  field 

measurement  in  addition  to  those  which  can  be  obtained 
on  the  surface.  These  measurements,  taken  in  an 

environment  where  there  is  -hon-equilibrium  chemistry, 
would  provide  important  information  with  which  to  model 
the  conditions  close  to  the  wall  in  modifications  to  the 
Navier/Stokes  equations  to  describe  flows  at  the 
continuum  end  of  the  transition  regime,  in  the  laminar 
continuum  flow  at  the  lower  altitudes,  where  non- 
equilibrium  air  chemistry  can  be  embodied  with  Navter- 
Stokes  codes,  questions  associated  with  catalytic  wall 
effects  still  remain.  However,  Navier-Stokes  solutions 
to  the  laminar  flow  over  the  Space  Shuttle  incorporating 
both  real  gas  and  viscous  interaction  effects  would 
provide  important  insight  into  the  modeling  of  these  flows 
with  less  demanding  techniques. 

The  prediction  of  boundary  layer  transition  and 
characteristics  of  transitional  flows  in  complex 
interaction  regions  is  one  of  the  most  difficult  and  yet 
important  problems  In  the  design  oi  hypersonic  vehicles. 
In  addition  to  the  sudden  changes  which  can  occur  in 
the  aerodynamic  performance  of  intakes  and  control 
devices  when  transition  occurs,  the  thermal  loads  on  the 
leading  edges  of  swept  wings  and  intakes  when  the 
boundary  layer  becomes  turbulent  may  be  a  key  feature 
controlling  the  acrothermal  design.  As  discussed  in  the 
review,  the  dynamic  loads  generated  in  transitional 
regions  of  shock  wave/boundary  layer  interaction  may 
also  lead  to  significant  design  problems.  The  tripping 
of  the  boundary  layer  at  the  lower  altitudes  by  surface 
roughness,  or  blowing  from  a  transpiration  cooling  system, 
is  another  important  phenomenon.  In  this  case,  however, 
there  is  some  reason  to  believe  that  it  is  possible  to  use 
transition  measurements  from  ground  test  facilities  to 
predict  transition  in  flight.  When  transition  occurs  in 
or  near  regions  of  shock  wave/boundary  layer  interaction 
as  in  the  flow  over  indented  nososhnpes,  these  flows  are 
not  only  difficult  to  describe  theoretically  but  also 
sensitive  to  Mach  number,  Reynolds  number  and  surface 
blowing  and  roughness.  The  use  of  wind  tunnel 
measurements  made  under  such  conditions  to  evaluate  a 
prediction  technique  must  be  handled  with  great  care. 

The  design  of  ablative-  or  transpiration-cooling 
systems  for  components  of  a  hypersonic  vehicle  is  an 
area  where  a  significant  research  effort  is  needed  to 
develop  new,  more  accurate  prediction  techniques.  In 
addition  to  problems  associated  with  the  inherent  three- 
dimensional  structure  of  the  transition  process  on  ablative 
heat  shields,  defining  the  roughness  height  and  surface 
characteristics  of  an  ablating  surface  lias  yet  to  be 
successfully  handled.  The  theoretical  techniques  used  to 
describe  the  ablation  of  a  heat  shield  in  the  presence  of 
blowing  and  rougiiness  are  at  best  highly  simplified,  and 
tiie  almost  complete  absence  of  experimental 
measurements  on  the  effects  of  combined  rougiiness  and 
blowing  in  high  Reynolds  number  flows  over  highly-cooled 
surfaces  does  not  help  this  situation.  Because  of  the 
rapid  titinning  of  tiie  boundary  layer  which  occurs  in 
regions  of  shock  wave/boundary  layer  interaction,  the 
effects  of  surface  roughness  can  exert  a  strong  influence 
on  the  local  iieating  rate  in  these  regions.  Also  tiie 
momentum  deficit  induced  at  tiie  base  of  a  rougli  wall 
boundary  layer  makes  sucli  a  flow  separate  more  readily 
when  subjected  to  strong  adverse  pressure  gradients.  The 
influence  of  transpiration  cooling  on  tiie  development  of 
a  turbulent  boundary  layer  in  a  siiock  interaction  region 
is  another  problem  which  could  be  of  importance  in  tiie 
design  of  ram  jets  for  hypersonic  propulsion. 
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In  hypersonic  flows,  regions  of- shock  wave-turbulent 
interaction  are  one  of  the  most'  stressing  aerothermal 
problems  as  well  as  one  of  the  most  difficult  to  compute 
accurately;  At  the  heart  of  the  problems  associated 
with  theoretical  description  of  these  flows  is  the 
description -of  the  development  of  the  turbulence  at  the 
base  of  the  boundary  layer  as  separation  occurs  and  as 
the  boundary  layer  undergoes  a  radical  thinning  in  a 
recompression  process  where  the  pressure  can-rise  two 
orders  of  magnitude  over  several  boundary  layer 
thicknesses.  There  are  questions  on  compressibility 
effects  in  these  flows  and  the  effects  of  the  intrinsic 
flow  unsteadiness  of  separated  regions.  Do  the  unsteady 
distortions  of  the  shock  waves  which  traverse  the 
boundary  layer  in  the  separation  and  reattachment  regions 
contribute  significantly  to.  -the  turbulence  generated  in 
these  flows?  In  fact,  can  these  flows  be  described  within 
the  framework  of  the  time  or  mass  averaged  Navier- 
Stokes  equations?  Are  there  Mach  number  effects  on 
the  separated  shear  layer  mixing  similar  to  those  which 
are  believed  to  exist  in  free  shear  layer  mixing  and 
turbulent  wakes?  It  is  clear  that  only  through  carefully 
directed  and  conducted  experimental  studies,  involving 
detailed  mean  and  fluctuation  measurements  in  the  flow 
field  as  well  as  on  the  surface,  can  the  insight  and 
information  be  generated  to  resolve  these  issues. 
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Figure  1  Performance  Map  of  Current  Hypersonic 
Facilities 
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Figure  2  Velocity-Altitude  Comparison  oi  Flight  Vehicles 
With  Shock  Tunne!  Capability 
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Figure  3  Problems  Resulted  From  Laminar  Viscoiis/lnviscid 
Interaction 


-  FAIRING  OF  COMPUTtD 

GROUND-TEST  DATA 

•  Sl!*.| }  APP*0XIMATI 

A  SWEPT  CYLINDER 

J _ I 

.1  1.0 


Figure  4  Comparison  of  Transition  Criterion  At  Shuttle 
Flight  and  Ground-Test  Conditions  (Ref.  6) 
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Figure  5  Batt’s  Transition  Correlation  (BMD-TR-81-58) 
(Ref.  53) 
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Figure  8  Gun  Tunnel  and  Operating  Cycle 
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Figure  9  Comparison  of  Performance  of  Shock  Tunnels  and  Free  Piston  Facilities 
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Figure  12  Comparison  of  Space  Shuttle  Orbiter  Data  From 
Calspan  Hypersonic  Shock  Tunnel 


Figure  10  Shock  Tunnel 


Figure  11  Comparison  of  Measured  and  Predicted  Heating 
Rates  for  STS-2:  (Zoby  JSR  Vol.  20  No.  3) 
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Figure  14  Skin  Friction,  Heat  Transfer,  and  Pressure 

Distributions  on  the  Flat  Plate-Wedge  Models 
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Figure  15  Skin  Friction,  Heat  Transfer  and  Pressure 

Distributions  on  the  Flat  Plate  <  20"R  Cylindrical 
Arc-Wedge  Models  (Ref.  11) 

(Mgo  *14.0,  Re/FT  *  7.2  x  104) 
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Figure  16  Comparison  Between  the  Integral  Solution  and 
Experimental  Measurements  in  a  Well  Separated 
Wedge-Induced  Interaction  Region  (Ref.  11) 
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Figure  18  Typical  Flow  and  Heat  Transfer  Laminar  Corner 
Flows 


Figure  21  Correlation  of  Transition  on  Slender  Cones 
Showing  Effects  of  Tunnel  Sice 
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Figure  19  Schematic  Representative  of  Model  With  Cross- 
Flow  in  Swepi-Shock  Laminar  Boundary-Layer 
Interaction 
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Figure  20  Shock-Wave  Structure  in  Axial  Corner  Comprised 
of  Two  Intersecting  Wedges 
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Figure  22  Correlation  of  Transition  Measurements  in  Cnlspan 
Shock  Tunnels  with  Ballistic  and  Downrange 
Measurements  (Ref.  57) 
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Figure  23  Transition  Reynolds  Number  vs.  Unit  Reynolds 

Number;  Present  Data  vs.  Potter  and  Sheet*  (Ref.  67) 
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Figure  24 


ASchlieren  Photography  Showing  Development  of  Transition  of  a  Conical  Boundary  Layer 
in  Hypersonic  Flow  (Ref.  57) 


Figure  25  Schlieren  Photographs  of  "Wave-Like''  instabilities  (Ref.  57) 
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Dependence  of  R«  Upon  £096  Mach  Number  For 
Rough  Surfeces,  After  Poll  (Ref.  44) 
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Transition  Correlation  on  Rough  Leading  Edge, 
After  Poll  (Ref.  44) 
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Figure  33 


21 


i«F 


20 


1.2  0 
□ 

O 

°*<  r 

o 


o 


Figure  34 


Flow  Over  the  Rough  MRV  Configuration  With 
a  30°  Flap  Deflection 
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Influence  of  Mach  Number  on  the  Distribution 
of  Heating  to  the  Smooth  Model 


(b!  MACH  13 

Figure  35  Effects  of  Mach  Number  on  Flow  Pattern  Over 
the  Indented  Nosetips  (Ref.  61) 
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Figure  44 


Heat  Transfer  Measurements  on  Scant  and  Smooth 
Hemispherical  Noietip  Showing  How  Small 
Blowing  Brings  Down  Heating  Levels  To  Smooth- 
Wall  Values 


(b>  BLOWING  SURFACE 


Figure  46 


Schematic  Diagrams  of  the  Flow  Structure 
Between  and  Above  the  Roughness  Elements  on 
Blowing  and  Non-Blowing  Surfaces 
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Figure  45 


Comparison  Between  the  Turbulent  Theories  of 
Lin  &  Crowell  for  Smooth-Wall  and  Rough-Wall 
Measurements  of  Holden  on  12"  Diameter  Hemis¬ 
phere  (M  =  11.2,  Re0  -  11  x  106,  K  =  12.5) 


Figure  47 


Schlieren  Photographs  of  the  Flow  Over  the  Sharp 
Biconic  Nosetip  Showing  the  Individual  Shocks 
for  the  Roughness  Elements  (Ref.  61) 
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Figure  49  The  Development  of  a  Wedge-Induced  Separated 
Flow  (M  oe  =  8.6  ReL  =  22.5  x  106)  (Ref.  75) 
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Figure  50  The  Development  of  a  Shock-Induced  . 

Separated  Flow  (Moo=  8.6  Re,  =  22.5  x  10b) 
(Ref.  75)  L 


Figure  51  Output  From  Skin  Friction  Gages  in  the 

Separation  Region  ELFSTROM 


Turbulent  Separated  Flows  (Ref.  75) 
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Figure  57  Typical  Surface  Heating  and  Surface  Piessuie 
Distributions  Through  Swept  Shock/Turbulent 
Boundary  Layer  Interaction 


(b)  SEPARATED  FLOW 


Figure  58  Schematic  Representation  of  Attached  ami 
Separated  Regions  in  Swept-Shock/Boundary 
Layer  Interactions 


Figure  59  Variation  of  Shock  Generator  Angle  to  Induce 
Incipient  Separation  with  Mach  Number 
(Ref.  56) 
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COMPONENT  OP  FREESTREAM  MACH  NUMBER 
NORMAL  TO  SHOCK  WAVE  (M  iln  ©  ) 

:igure  60  Correlation  of  Plataau  Pressure  Measurement  From 
Swept-Shock  Interaction  Studies  (Ref.  106) 


NORMAL  TO  SHOCK  WAVE  (M„  sin  ©  ) 

Figure  62  Correlation  of  Maximum  Pressures  Recorded  in 
Swept-Shock  Interaction  Regions 


Figure  61  Correlation  of  Peak  Heating  Rates  in  Skewed-  and 
Swept-Shock  Interaction  Regions 
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Figure  63  Streamwise  Distributions  of  Heat  T ransfer  and 

Pressure  Through  Skowed-Oblique-Shock  Boundary 
Layer  Interaction  (  Q  =  15°  V  =0°) 
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Figure  66  Damage  Resulting  From  Shock  Impingement  on  Vertical  Support 


